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ABSTRACT 


The analysis of this region of the Milky Way is based on star counts to magnitude 15, color excesses 
of stars determined here and elsewhere, and published counts of extragalactic nebulae. The obscuration 
is ascribed to two cloud systems: the first, at the distance of 500 parsecs, is chiefly responsible for the variety 
in the appearance of this part of the Milky Way; the second, at 800 parsecs, is nearly uniform in density 
throughout the low latitudes. However, an equally good representation of the counts and colors is ob- 
tained by employing constant coefficients of absorption. In all divisions of the region the star-density 
functions decrease with increasing distance from the sun; the gradient in the plane of the galactic equator 


is somewhat less steep than in the previously studied region of Auriga. 


An earlier study of this region of the Milky Way by one of us! was based on plates 
taken with the 4-inch Ross-Lundin camera of Harvard Observatory. To facilitate the 
analysis it was assumed that the southern part of the region is mostly free from dust 
clouds to within 6° of the galactic equator. Magnitude measures of stars then available 
in this area showed little or no color excess, though there was no assurance of the absence 
of dust beyond those stars. The star-density gradients for the apparently unobscured 
zones were then extended arbitrarily to the equator. It was pointed out that the absorp- 
tion introduced in the analysis is, therefore, a minimum throughout and is subject to up- 
ward revision when further data become available. 

The present study is confined generally to the southern galactic latitudes. The material 
supplied by the Harvard plates on 17 centers is supplemented by plates on 20 centers, 
which extend the original region to higher latitudes and somewhat to the east. The re- 
gion now under discussion lies generally between galactic latitudes +1° and —34°, and 
longitudes 75° and 105°. The supplementary plates were taken at the University of 
Illinois Observatory with the 4-inch Ross-Fecker camera, which is similar to the Harvard 
Observatory camera. The scale of the plates is 290’’ per millimeter. 


THE STAR COUNTS 


The method of counting is the same as that described in the earlier paper. Magnitude 
standards for calibrating the reference scale are available on plates containing Selected 
Areas 8, 9, 19, 20, 21, 22, 42, and 43 and are carried by overlapping plates to the re- 


1 Harvard Circ., No. 424, 1939. 
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mainder of the region. Standards are taken from Bergedorfer Spektral-Durchmusterung 
for the brighter stars and from the Mount Wilson Catalogue? for the fainter ones. No cor- 
rections are made to our counts for distance from center of the plate or for any other 
reason. 

In preparation for the analysis of the counts, the region is divided into areas arbitrarily 
lettered, in each of which the surface distribution of the stars seems fairly uniform and 


95 90° 
OBSCURATION BY CLOUD AT 500 PARSECS 


-20 


2 95° 90° 8s° 80° 75° 
OBSCURATION BY CLOUD AT 800 PARSECS 


Fic. 1.—Divisions of Cassiopeia region, showing obscuration at 500 and 800 parsecs. The —22.5° to 
—32.5° zones, which are unobscured, are not shown. 


different from that in the neighboring areas. The region is further divided by parallels of 
galactic latitude into zones 2.5° wide, in accordance with the usual practice. The arrange- 
ment of the divisions is shown in Figure 1, in which the zones of higher latitude are omit- 
ted to save space. The following changes in the designations of divisions for the present 
study may be noticed by comparing Figure 1 with the corresponding one in the earlier 


2 Mt. W. Contr., Carnegie Inst. Washington, No. 402, 1930. 
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paper. Area T is added. Area K is extended to the — 10° zone, and areas G and R are ex- 
tended to the —7.5° zone. The original areas J, Q, and R in the —2.5° zone prove to be 
considerably different from the similarly lettered areas in the higher latitudes and are 
now distinguished from them by the double letters /7, QQ, and RR. We have also omitted 
in the new diagram small patches in the —7.5° and — 10° zones which are darker than 


TABLE 1 
MEAN LOG N(m) FOR THE DIFFERENT DIVISIONS 

pen MAGNITUDE Oxsscu- | Appi- | Const. 

RATION TIONAL CoEFF. 

tone AT 500 | at 800 | pER1000 

zer* | 9.0 | 10.0 | 11.0 | 12.0 | 13.0 | 14.0 | 15.0 | PARSECS| Pamsucs | Panszcs 
ee) Aa + 2.5°| 82 | 0.52 | 0.99 | 1.37 | 1.73 | 2.07 | 2.38 | 2.57 | Om8 1=0 1™5 
Be. Me 0 97 .41 | 0.81 | 1.23 | 1.59 | 1.88 | 2.18 | 2.44] 1.5 1.1 2.6 
) 0 83 .69 | 1.16 | 1.56 | 1.96 | 2.36 | 2.67 | 2.88} 0.0 1.0 0.7 
| eae 0 90 .55 | 1.00 | 1.41 | 1.78 | 2.10 | 2.41 | 2.66] 0.5 2 1 
: ee 0 99 .55 | 0.98 | 1.41 | 1.76 | 2.07 | 2.38 | 2.65 | 0.5 1.2 1.7 
aa —2.5| 91 .72 | 1.19 | 1.55 | 1.91 | 2.27 | 2.58 | 2.90 | 0.0 1.0 0.7 
aes: — 2.5} 81 .52 | 0.90 | 1.30 | 1.66 | 2.03 | 2.34 | 2.60 1.0 0.9 1.5 
on — 2.5] 87 .39 | 0.94 | 1.36 | 1.76 | 2.14 | 2.45 | 2.70] 0.5 1.0 12 
OD ee ain —2.5| 94 .46 | 0.99 | 1.441 1.78 | 2.11 | 2.41 | 2.73 | 0.5 1.1 1.3 
) eer — 2.5} 95 .70 | 1.04 | 1.51 | 1.93 | 2.28 | 2.56 | 2.87 | 0.0 1.0 0.7 
RR... — 2.5 | 100 .68 | 1.08 | 1.48 | 1.85 | 2.20 | 2.50 | 2.78 | 0.2 1.0 1.0 
lit -— § 90 .78 | 1.18 | 1.51 | 1.90 | 2.28 | 2.62 | 2.89 | 0.0 0.8 0.5 
Pe 2 ee — § 84 .69 | 1.12 | 1.50 | 1.87 | 2.22 | 2.56 | 2.82 | 0.0 1.0 0.7 
ee -— 5 80 .60 | 0.98 | 1.42 | 1.78 | 2.19 | 2.46 | 2.75 | 0.3 1.0 1.0 
: er -— § 78 .55 | 1.02 | 1.53 | 1.90 | 2.32 | 2.62 | 2.91 | 0.0 0.8 0.5 
(2 ee -— § 95 .50 | 0.90 | 1.35 | 1.72 | 2.06 | 2.36 | 2.62 | 0.5 1.1 1.3 
5 94 .54 1.00 |} 1.48 | 1.84 | 2.19 | 2.51 | 2.79 | 0.0 1.0 0.7 
| —§ 100 .62 | 1.09 | 1.56 |.1.90 | 2.26 | 2.54 | 2.85 | 0.0 1.0 0.7 
— 7.5)| 91 | 1.15 | 1.53 | 1.90 | 2.31 | 2.62 | 2.90 | 0.0 0.5 0.3 
| ee — 7.5 | 85 .74 | 1.10 | 1.48 | 1.87 | 2.24] 2.55 | 2.85 | 0.0 0.9 0.6 
Ki. — 7.5 77 62 | 1.12 | 1.54 | 1.92 | 2.35 | 2.66 | 2.94| 0.0 0.5 0.3 
"| ao —7.5)| 96 36 | 0.90 | 1.36 | 1.67 | 1.98 | 2.33 | 2.63 | 0.5 1.2 1.3 
Ri: : — 7.5 | 100 65 | 1.03 | 1.46 | 1.83 | 2.21 | 2.53 |} 2.85 | 0.0 0.9 0.6 
a — 7.5 | 101 .53 | 0.87 | 1.38 | 1.72 | 2.04 | 2.36 | 2.72 | 0.5 1.0 1.2 
—10 76 |} .72 | 1.17 | 1.50 | 1.89 | 2.29 | 2.62 | 2.90 | 0.0 0.3 0.2 
P......| —10 95 .49 | 0.88 | 1.28 | 1.62 | 1.98 | 2.32 | 2.61 | 0.4 1.0 1.1 
—10 97 .71 | 1.05 | 1.40 | 1.76 | 2.12 | 2.48 | 2.77:| 0.0 1.0 0.6 
ae Ne —12.5| 93 .56 | 0.97 | 1.35 | 1.71 | 2.07 | 2.42 | 2.68 |. 0.0 0.8 0.5 
—15 93 .53 | 0.92 | 1.32 | 1.66 | 2.01 | 2.35 | 2.62 | 0.0 0.5 0.3 
wa aad |; —17.5} 92 48 | 0.87 | 1.27 | 1.62 | 1.95 | 2.30 | 2.57 | 0.0 0.3 0.2 
—20 90 .37 | 0.76 | 1.19 | 1.57 | 1.90 | 2.25 | 2.51 | 0.0 0.2 0.1 
—22.5| 90 .23 | 0.61 | 1.01 | 1.45 | 1.81 | 2.18 | 2.45 | 0.0 0.0 0.0 

BN prince —25 89 .17 | 0.60 | 1.01 | 1.44 | 1.80 | 2.16 | 2.42; 0.0 0.0 0.0 
ah OR —27.5| 89 .20 | 0.63 | 0.98 | 1.39 | 1.72 | 2.05 | 2.34} 0.0 0.0 0.0 
eee —30 94 .18.} 0.60 | 0.96 | 1.36 | 1.67 | 1.99 | 2.29} 0.0 0.0 0.0 
—32.5| 95 | 0.10 | 0.56 | 0.90 | 1.30 | 1.64 | 1.94 | 2.24!) 0.0 0.0 0.0 


* All galactic co-ordinates in this paper are referred to the conventional pole at right ascension 12540™, declina- 
tion +28°, 


their surroundings and are probably wisps detached from the main body of the nearer 
dust cloud to the north. Owing to their small areas, these patches are omitted from the 
analysis as before. 

The values of the mean log N(m) representing the counts in the different divisions ap- 
pear in Table 1, where N(m) is the number of stars per square degree brighter than ap- 
parent photographic magnitude m. Before the analysis of the counts is made, these are 
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converted to log A(m) at magnitude intervals of from m = 9.5 to m = 14.5, where A(m) 
is the number of stars per square degree between m — } and m + 3. 


ANALYSIS OF THE STAR COUNTS 


The analysis is made with the (m, log 1) tables whose construction and use have been 
described by Bok.* The amount of absorption introduced in the analysis is determined 
for the higher latitudes chiefly by the nebula counts and for the lower latitudes by the 
star colors that are available; these independent data will be discussed later. The absorp- 
tion over the entire region can be ascribed to two layers of dust clouds at the distances of 
500 and 800 parsecs, whose effectiveness in the various divisions is shown in Table 1 and 
also by the different degrees of shading in Figure 1. No considerable amount of obscura- 
tion is indicated south of the — 20° zone. In that zone the 800-parsec cloud makes its ap- 
pearance and increases in density until its maximum absorption of around 1.0 mag. is 
attained in the — 10° zone. From there to the equator this cloud remains practically uni- 
form in density except in areas G and K, where it is appreciably thinner. On the other 
hand, the 500-parsec cloud is irregular in density and is therefore chiefly responsible for 
the variety in the appearance of the Milky Way in Cassiopeia. Fragments of the nearer 
dust cloud appear as far south as the — 10° zone, in area P and probably in a few other 
places which are too small to analyze. Area H, a conspicuous dark spot in the photo- 
graphs, is the densest patch of the 500-parsec cloud south of the galactic equator. Im- 
mediately north of the equator the dark area A marks the beginning of the well-known 
“Cepheus flare” of heavy obscuration. 

The earlier study of the Milky Way in Cassiopeia was considered as provisional be- 
cause it introduced minimum absorption and did not arrive at a physical picture of the 
absorption. While satisfactory agreement between observed and computed log A(m) in 
the different divisions of the region was obtained by employing constant coefficients of 
absorption and superposing an absorbing cloud at 800 parsecs, both coefficient and ab- 
sorption in the cloud were required to increase in the more heavily obscured divisions. 
Fairly satisfactory agreement was obtained by introducing two clouds at 200 and 800 
parsecs, respectively; but here again the increase of the obscuration in the darker areas 
had to be divided proportionately between the two clouds. 

The present analysis is well served by two dust clouds which operate independently. 
It will be noticed in Figure 1 that each bright or dark area is an irregularity in one or the 
other of the two cloud layers and never in the two together. However, we find that the 
hypothesis of constant coefficients of absorption can produce an equally good agree- 
ment bet ween the observed and the computed log A(m). When the values of the unob- 
scured log A(m) are kept about the same as in the two-cloud analysis, the appropriate 
coefficients per kiloparsec of distance from the sun are shown in the last column of Table 
1. Moreover, we have not found any evidence from the available data independent of the 
star counts that seems to favor one construction over the other. 


COMPARISONS WITH RELATED DATA 


Color excesses of stars of spectral classes B8—A7 and chiefly between red magnitudes 
9.0 and 12.0 have been determined in Selected Areas 8, 19, and 43 by Dr. Lois Kiefer and 
Dr. J. C. Bell. The red magnitudes were measured on Eastman 103E plates, taken with 
the Ross-Fecker camera, the plates given 60 minutes’ exposure behind a ciné-red filter. 
Three or four such plates were employed for each Area. Magnitude standards were ob- 
tained by pairing each plate with a similarly exposed and developed plate centered at the 
north pole.4 The blue magnitudes and spectral classes of these stars were taken from the 
Bergedorf catalogue. The color excess for each star is the amount by which the observed 


3 Distribution of the Stars in Space, University of Chicago Press, 1937. 
* Harvard Ann., 89, No. 5, 1935. 
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color index (blue minus red magnitude) exceeds the normal color index to be expected for 
a star of the particular spectral class. 

In Table 2 the color excesses for the individual stars are assembled into means with 
respect to the spectral classes and red magnitudes of these stars. When the mean color 


TABLE 2 
COLOR EXCESSES 
Mean Color Distance Mean Color Distance 
a tral ye of Red Excess (Par- — of Red _ Excess (Par- 
Mag. (Mag.) secs) Mag. (Mag.) secs) 
Area Cc at +2.5° Area D at 0° 
10 11.5 +0.24 1150 
| + 830 34, | 11.5 | 40.37 
eee 10.5 33 630 27 11.0 + .25 870 
11 + 
10 95 + 19 10.5 + .28 690 
16 10.0 + .28 550 
7 11.5 + .37 600 
A1-A2 9.5 + .13 520 
7 11.0 + .33 500 6 90 + (10 440 
A3-A7 11.5 + .20 480 
14 10.7 +0.14 360 28 11.5 + .28 690 
_ 30 | 11.0 | + .22 580 
AI-A2...) Jig | 405 | + | 460 
Areas E and F at 0° 8 9.1 + .11 290 
20 11.5 .20 480 
B8-B9 9 11.0 +0.32 1260 || A3-A7...| 410 11.0 Fi 17 400 
10 | 11.6 | + .24 | 1200 11 10.5 | +0.13 330 
AO-Al... 7 11.1 + .il 1150 
5 10.6 + .09 910 
4 11.9 + .10 1600 
A2 | 48 11.4 + 23 690 Area G at —2.5 
6 | 12.0 | — .02 960 
U2 11.1 | -+0.06 | 480 7 | 11.6 | 40.30 | 1700 
9 11.0 + .43 1100 
B8-B9...| 17 | 105 | + .30 | 660 
me Beek 6 | 10.0 | + .25 550 
19 | 11.5 | 40.57 | 720 13 | 12.0 | + .17 | 1600 
11 11.5 + .27 1100 
20 11.0 + .53 600 AO-Al1... 
12 11.0 + .27 870 
9 10.0 + .61 350 - 
10 | 11.5 | + .60 | 440 6 | 11.8 .07 | 1600 
A1-A2. 8 | 11.0 | + .48 420 || A2....... {11 113 "20 690 
6 10.5 + .52 » 320 
9 11.5 + .46 330 12 12.0 — .14 960 
A3-A7... 11 11.0 + .60 220 A3-A7... 11 11.5 .00 790 
7 10.5 +0.49 200 t 11.0 +0.21 500 


excess is multiplied by the factor 3, the result is supposed to represent on the Interna- 
tional Scale the photographic absorption, K, by the intervening dust. The corresponding 
mean distance of each group of stars is computed in the usual way by comparing the ap- 
parent red magnitude with the absolute magnitude adopted for the particular spectral 
class and correcting for absorption. 

When the values of the photographic absorption derived from Table 2 are plotted for 
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each division against the distances, it is noticed that the scatter of the points is in some 
cases greater than in most of the previous color determinations made here. The only de- 
parture from former procedure is that the blue magnitudes are taken from the Bergedorf 
catalogue instead of being measured on our plates. Whatever the cause may be, the scat- 
ter is so pronounced for areas E, F, and G that the absorption in these areas cannot be 
regarded as very reliably determined. The results in areas D and H are satisfactorily 
consistent. In the latter area, however, the full amount of the absorption which we ex- 
pect behind the 800-parsec cloud makes its appearance at 200 parsecs from the sun and 
remains constant as far out as the points go, an effect that is equally difficult to reconcile 
with the constant-coefficient hypothesis. 

The values of the absorption at 800 parsecs, which we adopt from the evidence of the 
Illinois colors, are included in Table 3 for comparison with the absorptions employed in 
the analysis of the star counts. The agreement would be considered quite good if we 
could assume that the stars whose colors were determined in areas C, E, and F are mostly 
in front of the 800-parsec dust cloud and that those in areas D, G, and H are mostly be- 
hind that cloud. The stars observed for color in Selected Area 43 are somewhat too few to 
provide a definitive value of the absorption in the —17.5° zone; yet they suggest absorp- 
tion of the order of 0.5 mag., as compared with the value 0.3 mag., which we have 
adopted in the analysis. 

Many photoelectric color excesses of class-B stars by Stebbins, Huffer, and Whitford® 
are available in this region. We have converted them to photographic absorptions on the 
scale we employ here by multiplying by the factor 7, and we list the means for the various 
divisions in Table 3. The corresponding mean distances are calculated from the distance 
moduli given by these investigators. The comparison of these results with the absorptions 
introduced in our analysis is conveniently made under three heads. 

1. In six divisions within 5° of the galactic equator these B stars have a mean distance 
slightly exceeding 800 parsecs. We therefore compare their mean absorptions with those 
in our analysis at 800 parsecs. The differences (SHW minus Illinois analysis) are quite 
small, with the exception of area F; the average difference is +0™1. 

2. In four divisions from the —7.5° to the —15° zone the B stars considered have a 
mean distance somewhat less than 500 parsecs, and the majority are in front of the 800- 
parsec cloud. Thus it is not surprising that the average absorption for all these stars is 
0™5 less than is employed in our analysis at 800 parsecs. 

3. There are six of these B stars scattered in the —22.5° to —30° zones. Since their 
average distance is only 300 parsecs, the agreement of theif absorption with that em- 
ployed in the analysis is fortuitous. 

Table 3 also contains the values of the obscurations derived from Hubble’s® nebula 
counts, where we take his log N(m) = 1.90 as representing no obscuration. In the higher 
latitudes, where no considerable amount of dust is likely to lie beyond the range of the 
star counts, these values are employed to guide the amounts of absorption introduced in 
the analysis. No such agreement is expected, of course, near the galactic equator. 

As an example of the agreement between the absorptions by dust clouds introduced in 
our analysis and the values derived independently, we consider area D. This bright area 
along the galactic equator is an opening in the nearer dust-cloud system. No obscuration 
enters until the 800-parsec cloud is encountered; here K = 170. The Illinois color ex- 
cesses of 226 stars between 300 and 900 parsecs suggest the value of K = 0™8 produced 
by a cloud at the distance of 600-800 parsecs; or they may be represented by a coefficient 
of 1"2 per kiloparsec. The photoelectric color excesses of 15 stars mostly between 400 and 
1600 parsecs give the average K = 0™9. W. Becker,’ from color excesses of stars in the 
northern part of Selected Area 19, finds no reddening of the stars as far as 1000 parsecs 


’ Ap. J., 91, 20, 1940; Mt. W. Contr., No. 621, 1940. 
6 Ap. J., 79, 8, 1934; Mt. W. Contr., No. 485, 1934. 7 Zs. f. Ap. 17, 285, 1939. 
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from the sun, where absorption sets in and increases to K = 1™0 at 2000 parsecs. Bok’s*® . 


color excesses of 81 stars show no absorption as far as 800 parsecs and, from there on, 
absorption whose rate of increase is represented by the coefficient 1™2/kpc. Zug® has 


TABLE 3 


OBSCURATIONS FROM STAR COLORS AND NEBULA COUNTS 
COMPARED WITH ADOPTED OBSCURATIONS 
From SHW Cotor Excesses ADOPTED IN ANALYSIS 
ILLINOIS NEBULA 
AREA ZONE CoLors Counts 
(At 800 Mean Mean (ToTaL Additional ; 
ParsEcs) Obscu- | Distance | Osscura- (At 800 
ration (Parsecs) TION) Parsecs) 
+ 2.5% 0.9 18 1.9 0.8 1.0 
0.0 0.5 7 1.9 630 >3.0 0.5 1.2 
—10.0j......... 4 0.4 500 0.9 0.0 1.0 
—30.0 }........ 0.0 0.0 0.0 


measured the color excess of one open cluster, NGC 7654, in area D within the range of 
our counts; the distance is 1400 parsecs. Multiplying his value by 1.9, we obtain K = 170. 
With allowance for possible distance errors, we find reasonable agreement betweeen all 
these results and the value K = 1™0 at 800 parsecs, which we have employed in the 


analysis. 


8 Ap. J., 90, 249, 1939. 


® Lick Obs. Bull., No. 454, 1933. 
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THE GRADIENTS OF STAR DENSITY 


Table 4 shows the actual values of log A(m) for the different zones, provided that the 
correct amounts of absorption by the intervening dust have been adopted in the analysis, 
The same provision applies to the density factors of Table 5. These factors are fractions 

TABLE 4 
ADOPTED VALUES OF UNOBSCURED LOG A(m) 


MAGNITUDE 
ZONE 
9.5 10.5 11.5 12.5 13.5 14.5 
1.16 1.59 2.00 2.35 2.62 2.86 
ay See 1.06 1.53 1.94 2.30 2.58 2.80 
Seas 1.00 1.48 1.89 2.25 2.55 2.78 
0.97 1.44 1.86 2:22 2.52 2.77 
0.88 1.34 1-77 2.14 2.39 2.73 
0.82 1.26 1.75 2.02 2.58 
0.71 1.14 1.48 1.87 2.12 2.39 
0.65 1.06 1.45 1.76 2.06 
4 eee 0.59 1.00 1.38 172 2.00 2.25 
0.41 0.79 1.20 1.57 1.89 2.17 
0.42 0.81 1.18 1.54 1.89 2.12 
0.41 0.76 1.08 1.45 1.77 2.00 
0.32 0.68 1.08 1.36 1.67 1.91 
TABLE 5 
DENSITY FACTORS, D(r), FOR NO OBSCURATION 
DISTANCE IN PARSECS 
ZONE 
125 200 300 500 800 1250 2000 3000 
- SS ane 1.00 1.00 0.90 0.90 0.75 0.45 0.20 0.03 
at ke 1.00 1.00 90 90 as 45 20 03 
| 1.00 1.00 90 85 45 20 03 
1.00 1.00 90 85 45 20 03 
1.00 0.90 85 70 .60 .40 20 03 
1.00 0.90 80 65 .50 30 10 03 
1 1.00 0.80 70 45 .30 20 10 02 
See 1.00 0.80 65 40 .30 15 06 0.02 
1.00 0.60 30 18 .2 15 
1.00 0.60 30 20 16 
27.5 1.00 0.60 30 20 .10 06 
1.00 0.60 30 20 .10 .06 
—32.5..... 1.00 0.50 0.25 0.15 0.10 0.05 


of the star density in the immediate neighborhood of the sun, which is taken as unity. 
The factors at 3000 parsecs are not considered very reliable; they are included because 
they were used in the analysis. The negative gradients with increasing distance from the 
sun show that the sun is situated in a region of high star density with respect to the direc- 
tion of Cassiopeia. 
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Since its inauguration by Bok" the analysis of the Milky Way based on extensive star 
counts over the various regions has proceeded far enough to suggest certain features of 
the galactic system within a radius of some 2000 parsecs around the sun. We mention a 
few of these which can be noticed in the results of analyses at this observatory. 

1. The gradients of star density are generally negative with increasing distance from 
the sun. Thus the sun seems to be situated amid a high congestion of stars as compared 
with its less immediate surroundings. 


TABLE 6 
COMPARISON OF DENSITY FACTORS IN 0° ZONE 
DISTANCE IN PARSECS 

GALACTIC 

REGION LONGITUDE 
125 | 200 | 300 | 500 | 800 | 1250 | 2000 3000 
ARG, Oe. Ss 130° 1.00 | 1.00 | 0.90°} 0.90 | 0.60 | 0.30 | 0.15 | 0.05 
Cassiopeia........ 90 1.00 | 1.00 | 0.90 .90 75 45 .20 03 
Ophiuchus........ 340 1.00 | 1.00 | 1.00 | 0.90 | 0.70 | 0.50 | 0.50 | 0.40 


2. The gradients appear to be steeper in the direction of the galactic anticenter. At 
least this is true for the regions of Auriga," Cassiopeia, and Ophiuchus,"” whose density 
factors in the galactic plane are shown in Table 6. 

3. Exceptions to the rule of 1 are found in the Scutum cloud and the little Sagittarius 
cloud, where the density factors exceed the unit value assigned to the sun’s vicinity. 
Other probable exceptions are found in Aquila!’ and in Selected Area 9, in Perseus. 

4. Features often associated with the galactic system as a whole, including the dust 
clouds as well as the star clouds, are prominent in the “local system”’ within 2000 parsecs 
of the sun. 


10 Harvard Circ., No. 371, 1931. 12 Ap. J., 96, 224, 1942. 
u Ap. J., 94, 482, 1941. 13 Ap. J., 94, 493, 1941. 
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ABSTRACT 


From 341 radial-velocity observations of the spectroscopic binary 8 Coronae Borealis made in the in- 
terval 1931-1942 with the three-prism Mills spectrograph, the writer has derived the orbital elements of 
the system, whose period is 10.496 years. Some evidence is found for the presence also of a shorter period 


of 320 days. 


The star 8 Coronae Borealis (1900. a = 15237, 6 = +29°27’, ptm 3.72, class FOp) 
was announced as a spectroscopic binary by Dr. J. H. Moore! in 1907. In 1912 Mr. J. B. 
Cannon? published orbital elements for the star, basing his investigation upon 153 spec- 
trograms taken with one-prism dispersion at the Dominion Observatory in Ottawa in 
the interval 1910-1912. On account of the large accidental errors to which his velocity 
determinations were subject, he had considerable difficulty in finding a period for the 
velocity oscillation, but by combining a period of 40.9 days with a secondary one of 
490.8 days he was able to represent the observations fairly well by an approximate set 
of elements. 

In 1931 a request came to the Lick Observatory from the parallax observers at the 
Allegheny Observatory to check, if possible, Mr. Cannon’s two periods; and Dr. Moore 
turned over to me the investigation of this star. Observations were begun in May, 1931, 
and continued until October of that year, when the star could no longer be observed, the 
interval covered being almost five times the period of 40.9 days. The 58 velocities ob- 
tained that season hovered fairly constantly around the value of — 16 km/sec and could 
not be represented by Cannon’s orbital elements. One striking fact, however, was re- 
vealed, namely, that there was a progressive run in the individual residuals from the 
season’s mean velocity, the range being about 2.5 km/sec, indicating the presence of a 
true variation with a period greater than six months. 

Since observations throughout a single night showed no variation in velocity, a very 
short period could be left out of consideration. In February, 1932, another set of observa- 
tions was begun and was continued until October of that year, 55 additional plates being 
obtained in the interval. Again the same characteristics as in the previous year were evi- 
dent, the range of the residuals derived from the annual mean velocity being of the same 
magnitude and exhibiting the same progressive tendency as did those of the year before. 
Neither series of observations could be represented by either of Cannon’s periods, and it 
was at once evident that the period would have to be of the order of ten or eleven years 
in order to represent the velocities of the present series, together with those obtained 
here in the years 1900-1902 with the original Mills spectrograph and those since 1905 
with the new Mills spectrograph. Observations were therefore continued through the fol- 
lowing ten years, thus affording a continuous series extending over one cycle, and these 
form the basis of the present investigation. 

In the interval from May, 1931, to September, 1942, 341 spectrograms were taken by 


* Contributions from the Lick Observatory, Ser. II, No. 8. 
1 Lick Obs. Bull., 4, 162, 1907. 


* Pub. Dom. Obs., Ottawa, 1, 375, 1912. 
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the writer with occasional help from Dr. Moore and Dr. G. F. Paddock. They were made 
with the three-prism Mills spectrograph, which has a linear dispersion of about 10 A/mm 
at \ 4500. Prior to the spring of 1941, Eastman 40 plates were employed, for which the 
normal exposure time for this star was 56 minutes. At that time the prisms were coated 
with a film of lithium fluoride, and the collimator and camera lenses with one of mag- 
nesium fluoride, which effected a 60 per cent decrease in the required exposure time, 
thus permitting the use of Eastman 33 with an exposure time of only 30 minutes. The 
slit-width throughout the program was 0.025 mm. 

The spectrum of 8 Coronae Borealis contains numerous excellent lines suitable for 
radial-velocity measurement and is classified by Harvard as FOp. It is peculiar in “the 
relative intensities of the lines and resembles the spectrum of a Persei, except that the 
lines appear to be less narrow.” Mount Wilson assigns the spectrum to an earlier class, 
A7s. No secondary spectrum is visible, nor has any evidence of change in width of the 
lines been found. The measurable region of the spectrograms lies between Ad 4390 and 
4620, the line \ 4500 being central. In this region the spectrum of 8 Coronae Borealis has 
about 200 measurable lines, all of which may readily be identified by means of the data 
given in Dunham’s description of the spectrum of a Persei.* Of these, 60 lines that were 
free from blends were selected for the measurement on a plate of this star which was 
taken for use as a standard on the Hartmann spectrocomparator. This spectrogram, ob- 
tained on April 7, 1934, was measured by me several times on a Toepfer measuring en- 
gine, yielding the adopted velocity of — 22.10 km/sec, corresponding to the displacement 
of the stellar lines with reference to those of the comparison spectrum. The other 340 
plates were measured against this standard, and data for these are given in Table 1. 

A velocity-curve was obtained by plotting the data for the series of observations 1931- 
1942, given in Table 1, and has a double amplitude of 18 km/sec with a period of approxi- 
mately 10.5 years. Upon this there appears to be superimposed a variation having a 
double amplitude of nearly 3 km/sec with a period of a little less than a year. 

Preliminary elements for the long-period system were first derived from the velocity- 
curve drawn through normal points formed by taking the annual means of the velocities, 
the eight early observations of 1900-1907 being utilized in determining the period. The 
period being regarded as known, differential corrections to the other five elements were 
computed by the method of least squares, which yielded the following values: 


PRELIMINARY ELEMENTS 


P= 10.504 years 
T = 2428964.3 + 49.5 JD 


K= 9.09 + 0.42 km/sec 

e= 0.51 + 0.03 
o= 176.2 + 5.5 

—19.11 km/sec 


The 341 observations were then regrouped into 57 normal places, consisting of month- 
ly means of the velocities. On the basis of the preliminary elements derived above from 
annual means, the residuals given for the 57 normal places were formed, and a second 
set of elements was derived by a least-squares solution, in which the period was in- 
cluded. These elements, which resulted in a 50 per cent decrease in the sums of the 


5’ Contr. Princeton Obs., No. 9, 1929. 
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136 F. J. NEUBAUER 
TABLE 1 
RADIAL-VELOCITY OBSERVATIONS OF 8 CORONAE BOREALIS, 1931-1942 
Vel. Vel. 
Date JD (Km/Sec) Date JD (Km/Sec) 
2420000-+ 2420000+ 

1931 May 5...... 6466.724 —16.30 Pets. 6616.617 —14.16 
.884 17.32 19.632 14.92 
.923 16.43 26.622 13.86 
67.729 16.90 29 .606 13.47 
.886 16.68 33.606 13.97 

68.746 16.66 
71.823 16.24 1932 6756 .972 15.21 

77.742 16.74 
.788 16.99 Mar. 282. 94.923 16.17 

.850 17.46 
88.822 17.10 6805 .826 15.50 
.862 17.76 871 15.90 
.900 17.25 6.843 15.93 
908 16.48 
883 17.56 15.55 

.878 16.26 May 10...... 37.882 16.74 
Te? 96.872 14.96 .923 16.32 
6503. 871 15.20 38.919 16.31 
a 4.805 16.06 .960 16.08 
15.93 39.741 16.39 
9.712 16.68 16.10 
.766 15.48 .858 45:73 
14.758 15.16 .906 16.16 
.803 15.24 .956 16.04 
40.931 16.27 
23.888 15.78 43.853 16.18 
24.815 16.67 44.738 16.16 
34.713 16.16 .787 15.40 
37.742 15 53 15.65 
.790 15.60 (Te 46.785 15.84 
45.682 15.29 53.837 45.52 
49 686 14.97 54.748 15.18 

.721 14.90 
54.717 15.39 15.55 
(fe 60.671 15.29 .832 15.48 
66.683 15.29 67.785 15.74 
72.691 14.98 .834 16.10 
« 81.650 14.55 68.705 16.15 
.746 16.33 
98 . 654 15.26 78.805 15.34 
6600 .641 15.16 15.85 
1.634 16.00 82.739 15.93 
a 2.646 15.39 .790 15.66 
5.639 14.97 84.799 15.44 
7.643 14.71 87.790 15.25 
8.610 15.08 .840 14.71 

9.628 14.06 
10.625 14.64 91.840 14.26 
12.629 —14.87 .895 —14.31 
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Vel. Vel. 
Date JD (Km/Sec) Date JD (Km/Sec) 
2420000+ 2420000+ 
1032 July: 6895 .705 —15.17 1964 Ans. 7534 .897 —15.20 
751 15.35 .974 15.46 
.736 14.91 .888 15.20 
.782 15.06 20. 53.852 16.15 
.913 15.54 
.677 13235 May 18...... 75.852 15.69 
50.654 14.69 .904 14.26 
.699 14.30 
7655 .709 13.44 
1083 27. 7130.940 15.25 atte 15.44 
.990 15:17 
19033) Mare tt: 7872 .871 15.70 
40.914 14.77 .922 15.31 
.967 15.43 75.883 16.62 
50.885 14.69 93.024 16.46 
.948 15.19, 
May 3...... 7925.913 16.66 
64.875 15-77 .962 16.85 
.923 15.43 30. 884 16.75 
971 15.59 .927 16.75 
75.886 16.04 36.904 16.11 
.950 15.22 50.816 16.43 
May 19. ....... 7211.748 15.56 So 58.866 16.73 
.798 15.32 68.857 16.36 
fee 16.730 15.38 .914 16.98 
.773 15.80 75.756 16.97 
Fane’ 44.756 15.69 93.798 17.27 
.801 15.89 843 16.96 
14.49 65.656 17.67 
.828 15.00 
58.715 15.71 1936 Feb. 4...... 8202 .929 18.58 
.758 14.66 18.88 
60.757 14.12 25.947 19.26 
.799 15.13 .989 19.23 
65.806 14.69 
67.805 15.16 Mat. 32.923 20.92 
68 .796 13.68 33.007 20.41 
84.714 14.84 42.912 19.44 
.955 19.03 
7313 .661 13.84 48 .945 20.00 
14.654 14.80 
Pope. 72.820 20.34 
22.960 13.69 863 19.79 
43.612 13.83 83.751 19.74 
63.600 14.45 june 8328 .852 19.96 
36.826 19.52 
7469 .006 16.59 43.828 19.25 
.055 17.19 
1937 Apr. 16...... 8639 .938 25.84 
Bar. 98.938 16.56 .982 25.51 
.987 —15.84 49.949 — 26.08 
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138 F. J. NEUBAUER 
TABLE 1—Continued 
Vel. Vel. 
Date JD (Km/Sec) Date JD (Km/Sec) 
2420000+- 2420000+ 
19a easy 8658 .938 —26.76 1959 2 9385 .890 — 24.50 
.992 26.11 .929 24.75 
70.697 25.54 .973 25.83 
26.24 
801 26.36 9499 697 22.31 
.850 26.70 .740 23.05 
.892 26.31 9506 .666 22.98 
.935 26.51 
.977 26.75 13.648 22.19 
26. 8740 .727 30.28 1940 Feb. 9...... 9669 .052 19.53 
.770 31.70 
Mare 91.925 18.87 
dept. 86.654 32.84 96.989 19.43 
.697 32.22 97.048 19.63 
9701 .958 19.00 
1938 Feb. 20...... 8950 .023 30.58 2.023 19.27 
53.939 32.16 12.958 18.27 
90.015 32.74 32.971 19.00 
9006 .019 33.04 33.017 19.27 
40.972 19.49 
998 34.19 May 21...... 70.847 17.90 
Lae 33.918 34.82 .890 17.95 
36.938 33.47 
42.827 33.22 June 13... 93.784 17 
49.833 34.43 
July 16...... 9826.755 17.47 
58.836 33.03 .795 17.04 
1939 Mar. 14...... 9336 .906 25.33 .778 17.03 
.956 25.69 cis 29.785 17.15 
38.910 25.10 
.955 25.30 63.670 16.59 
51.848 25.81 
53.885 26.11 Sept. 10.3... 82.651 16.80 
.928 25.47 
.979 26.30 2430000+- 
1941 Apr. 0097 .989 16.41 
57.793 25.29 104.973 16.47 
.863 25.58 5.013 15.94 
.906 25.00 
60.921 26.27 21.912 16.09 
971 25.40 
61.780 25.85 June 54.828 16.09 
.826 25.67 are 58.815 15.97 
67 .966 25.85 .848 1577 
ie 81.901 24.95 July 10...... 85.732 15.28 
82.892 25.80 
.935 25.65 1942 Mar. 1...... 0420 .020 14.80 
21.954 14.85 
85.725 pA 23.954 14.74 
.769 .978 15:35 
.810 26.979 15.52 
.847 —24.17 35.902 15.41 
.921 —14.65 
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Vel. Vel. 
Date JD (Km/Sec) Date JD (Km/Sec) 
2430000-++ 2430000+ 
1942 Apr, £9 0464 .861 —15.07 1942 June 19...... 0529.705 —13.42 
.885 14.23 31.715 13.35 
74.936 14.53 .750 13.89 
75.803 13.99 40.679 14.06 
.839 14.47 hae 13.67 
May 3...... 82.779 13.59 July 5...... 45.780 13.73 
815 14.25 .809 14.02 
85.776 14.27 47 .698 13.75 
.809 14.26 429 13.34 
.840 14.16 49 693 13.58 
.882 14.13 .723 13.99 
915 13.73 54.716 13.40 
96.813 14.33 .742 13.78 
.851 13.51 59.670 13.87 
98 .698 13.57 .692 13.30 
.738 13.95 BGS 66.734 14.02 
.772 13.71" 13.73 
.802 13.07 
.833 13.48 73.685 13.95 
.861 12.88 ° 711 14.39 
.891 12.86 75.663 14.49 
14.15 
0513 .726 13.93 96.673 15.17 
.759 13.75 .699 14.52 
.798 13.03 
.835 13.03 Sept, 0612 .640 14.40 
23.687 14.00 .664 — 14.65 
721 —13.87 


squares of the residuals, are adopted as the final ones for the long-period system and are 


as follows: 


P= 


lI 


asin? 


f(m) = 


FINAL ELEMENTS 
10.496 + 0.001 years 


2428971.3 


9.19 


+ 18.09 JD 


+ 0.23 km/sec 


+ 0.025 
+. 


—18.03 km/sec 
442,000,000 km 


0.24 


© 


The probable error of a normal place of unit weight is +0.67 km/sec. Data for the 57 
normal places are entered in Table 2, the first two columns of which contain, respective- 
ly, the mean date of the observation in Julian Days and the corresponding mean radial 
velocity. In the third column are given the residuals O — C for the normal places com- 
puted on the basis of the final elements, the number of observations combined to form 
the normal place or its assigned weight being listed in the final column. The fourth col- 
umn will be explained later. 
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TABLE 2 
NORMAL PLACES 
Corrected 
Vel. Orbit for 
JD (Km/Sec) (O—C) Secondary Weight 
(0’—C’) 
2420006+- 

—16.91 —2.00 —0.60 13 
15.95 —1.12 +0.18 12 
ee 15.61 —0.85 +0.33 8 
15.15 —0.49 +0.36 
.. 14.96 —0.37 +0.15 13 
15.21 —0.83 —1.22 1 
Sa 16.17 —1.72 —0.32 1 
16.01 —1.69 —0.29 6 
16.00 —1.68 —0.48 18 
15.65 —1.34 —0 18 
i! 14.86 —0.56 +0.26 7 
14.28 +0.07 +0.29 4 
15.08 —0.66 +0.50 6 
15.61 —1.17 +0.13 5 
15:52 —1.00 —0.10 4 
15.79 —1.23 —0.75 2 
Ue 14.74 —0.15 —0.03 11 
14.11 +0.58 —0.30 
14.14 +0. 64 —0.64 2 
7460-969 16.54 —1.43 —0.14 4 
15.45 —0.14 +0.65 6 
14.98 +0.45 +0.74 2 
es 14.44 +1.31 +0.19 2 
16.02 +0.92 +1.47 4 
16.59 +0.72 +0.40 6 
68.098....... 16.76 +0.81 —0.19 4 
17.12 +0.64 —0.68 
18.38 —0.07 +0.73 
8214.460....... 18.99 +0.79 +1.09 4 
40.148....... 19.96 +0.27 +0.17 5 
19.95 +0.75 +0.15 3 
$3530 19.55 +1.97 +0.69 3 
2647956... 25.81 +1.67 +0.27 3 
68 26.45 +1.72 +0.94 9 
8740.748....... 30.99 —1.20 +0.16 2 
BE 31.90 —1.76 —0.65 2 
8981-961. 32.43 +0.35 —0.93 2 
9056 24S. . 34.12 —2.11 —0.71 6 
33.03 —1.32 +0.06 1 
9345.671....... 25.63 —0.88 +0.22 8 
25.46 —1.32 +0.08 20 
| —22.65 —1.16 —0.90 4 
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Corrected 

Vel. Orbit for é 
JD (Km/Sec) (O-C) Secondary | 
(0’—C’) 
2420000-+- 

9669 .052....... —19.53 —0.45 +0.73 1 
9710.048....... 19.08 - — .38 +0.93 6 
19.25 — .93 +0.33 3 
17.74 + .17 +1.01 3 
9828.098....... 17.34 + .10 +0.22 6 
73.166... 16.70 + .37 —0.33 2 

2430000+- 
0105.764....... 16.24 — .59 +0.07 5 
164.556....... 15.64 — .22 —0.36 4 
0426.958....... 15.83 — .40 +0. 26 7 
487 .448....... 13.91 + .60 +0.44 21 
190s 14.46 + .05 —0.87 11 
554.062....... 13.71 + .70 —0.68 12 
SOO: GIS —14.46 —0.08 —0.86 8 


The velocity-curve of the long-period system of 8 Coronae Borealis in the interval 
1931-1942 is given in Figure 1. It was computed from the final elements; and upon the 
diagram the 57 normal points representing the 341 observations have been indicated by 
small circles. The small secondary variation in velocity, to which reference has already 
been made, is apparent on the diagram, as well as in the run of the residuals in the third 
column of Table 2, especially for the first four years of observation, 1931-1934. Since the 
amplitude of the variation is small, considerable difficulty was encountered in deriving a 
curve which would represent the residual velocities from the long-period velocity-curve. 
After numerous trials the best representation was obtained by the curve for the second- 
ary variation reproduced in Figure 2, in which the 57 normal places are represented by 
circles and the 9 normal points into which they were combined by dots. From this curve 
the following approximate elements were derived: 


P’ = 320.13 days w’ = 102° 

T’ = 2426475 JD a’ sin i’ = 4,400,000 km 
K’ = 1.4 km/sec f'(m) = 0.000034 © 
= 0.7 


When the 57 normal places were corrected for the effect of the secondary variation 
having the elements derived above, the residuals O’ — C’ with reference to the long- 
period velocity-curve, given in the fourth column of Table 3, were obtained. By this pro- 
cedure the probable error of a normal place of unit weight was reduced to +0.4 km/sec, 
while the sum of the squares of the residuals was reduced from (O — C) 336.89 to 
(O’ — C’) 97.35. Although the 70 per cent reduction in residuals thus effected would seem 
to offer evidence of considerable weight in favor of the reality of the short-period varia- 
tion, it cannot be considered conclusive, for the secondary curve was obviously so de- 
rived as to effect a reduction in the original residuals of the normal places. After the re- 
moval of the secondary variation, the 57 normal places were again plotted on the long- 
period velocity-curve, and the diagram thus obtained is reproduced in Figure 3. While 
the representation of the observations has obviously been somewhat improved over that 
shown in Figure 1, especially in the portion of the curve for the years 1931-1934, it must 
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be admitted that little or no improvement appears to have been effected in the repre- 
sentation of the observations for other parts of the curve obtained in later years. 

Since the period of secondary variation indicated is nearly a year, consideration has 
been given to the question of its possible origin in instrumental or other errors of a sea- 
sonal nature. Thus it is known‘ that small errors in the velocity determinations made 
with the Mills spectrograph are produced by changes in collimation due to flexure in the 
tube of the 36-inch refractor. Such errors depend upon the declination and hour angle of 
the star and are of opposite signs for spectrograms obtained with the telescope on oppo- 


TABLE 3 


RADIAL-VELOCITY OBSERVATIONS OF 8 CORONAE BOREALIS 
PRIOR TO 1931 AND IN 1943 


Date JD (O-—C) (0’—C’) Remarks 
1900 Apr. 18....| 2415128.89 —17.8 Pub. Lick Obs. 
Vol. 16, 1928 
May 13.... 53.79 18.1 
01 May 6.... 5511.83 18.1 
02 Feb. 13.... 5794.06 15.1 
05 Apr. 11.... 6947 .03 23:2 
07 Mar. 29.... 7664 .06 29.4 
Apr. 8.... 74.96 31.0 
July 17... 7774.71 26.8 
13 Feb. 13.... 9798 . 64 15.7 Cape* 
Mar. 19.... 9826.57 22.3 Cape 
16 June 16....} 2421031.30. 26.3 Cape 
30 Feb. 16.... 6023 .975 17.68 —0.37 —0.37 | Lick 
Mar. 11.... 46.924 17.27 —0.22 + .89 
43 July 17....| 2430922.708 13.93 +0.43 + .25 
Sept. 4.... 971.649 15.86 +1.45 + .07 
.674 15.53 +1.08 — .30 
981.641 —15.98 +1.55 +0.23 


* Observations made at the Cape of Good Hope, Cape Ann., 10, Part VIII, 1928. 


site sides of the pier, positions of the instrument required in observing the same star at 
different seasons of the year. The effects of such errors, however, have been carefully de- 
termined for all positions of the telescope, and corrections for them are applied to all ob- 
servations made with this spectrograph, including those of the present investigation. 
Constants used in computing the “reduction to sun,” errors which would introduce 
seasonal effects in the velocity, also have been checked. Furthermore, an examination of. 
the radial-velocity measures of other stars made with this instrument in the interval 
1931-1942 shows no evidence of the presence of variations of the character here consid- 
ered. The evidence, therefore, appears to be quite definite that secondary variations of 
several kilometers per second do occur in the radial velocity of 8 Coronae Borealis, al- 
though the available data are not sufficient to determine satisfactorily whether these are 
strictly periodic. In order to throw further light upon the nature of these small varia- 
tions, further observations will be required, and to this end it is planned to obtain several 
series of spectrograms within the next few years. 


4 Pub. Lick Obs., 16, xxv, 1916. 
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Observations made with the three-prism dispersion prior to and later than those listed 
in Table 1, which were used in the above discussion only for a determination of the 
period of the long-period system are given in Table 3. Data for the first four, taken with 
the original three-prism Mills spectrograph, are published in Publications of the Lick Ob- 
servatory, Volume 16. The radial velocities entered for the Lick observations 1905-1930 
are my measures of these spectrograms made on the Hartmann spectrocomparator 
against the same standard plate that was used in the measurement of the plates listed in 
Table 1. After completion of the present investigation, four spectrograms were obtained 
in 1943, and these are included in Table 3. The residuals (O — C) with reference to the 
long-period curve, and the (O’ — C’) after correction for the secondary variation are 
entered in the fourth and fifth columns. On account of the uncertainty in the short-period 
O’ — C’ residuals have not been given for observations earlier than 1930. 

Mr. Cannon’s 153 one-prism observations of this star were made in the interval Febru- 
ary 23, 1910—July 23, 1912, during which time the velocity computed on the basis of the 
elements here derived for the long-period system varies from —15.2 to —14.3 km/sec. 
His observations for this interval range from —25.2 to —18.9 km/sec and, when com- 
pared with the computed values, yield a mean residual of (O — C) = —7.03 km/sec. 

The trigonometric parallax of this star given in Schlesinger’s General Catalogue of 
Stellar Parallaxes (1935) is 07033, which yields a value for a sin i of 0709. There is thus a 
possibility that the two components could be separated near apastron with a large tele- 
scope, the favorable time for such observations being in 1944-1946. In this connection 
it should be noted that the Allegheny observers have found from their parallax plates 
evidence of the orbital motion of the brighter star in a period which is of the order of that 
obtained for the long-period system in the present investigation. 
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RADIAL VELOCITIES OF 283 STARS OF SPECTRAL 
CLASSES R AND N* 


RoscoE F. SANFORD 
Mount Wilson Observatory 
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ABSTRACT 


Radial velocities of 283 stars of spectral classes R and N are given in Table 1. By intercomparison 
with velocities obtained from spectrograms of high dispersion, these are found to be free from systematic 
errors. Thirty-four class-N stars (28 known to be light-variables) with emission lines of hydrogen have 
a mean difference of —20.4 + 0.9km/sec between the velocity from the emission lines (mostly Ha 
alone) and that from the absorption lines. 

Galactic concentration is marked for class-N stars, slight for class-R. 

Mean absolute magnitudes of —2.3 + 0.2 for class-N stars and of —0.4 + 0.4 for class-R stars have 
been derived from differential galactic rotation. 

Interstellar sodium lines —The radial velocities of 2 class-N stars and 4 class-R stars are sufficiently 
large to separate the stellar from the interstellar sodium lines. Velocities and approximate values of the 
equivalent widths have been derived from the interstellar lines. The distances and absolute magnitudes 
have been obtained from the equivalent widths. The mean magnitude for the 2 class-N stars is —1.8 and 
for the 4 class-R stars, 0.0. ; 

Stars of high velocity.—Table 4 contains 1 class-N star and 7 class-R stars, whose velocities, when freed 
from solar motion, are numerically larger than 100 km/sec. All 7 class-R stars have spectra with strong 
g-bands, but with few atomic lines. 


A catalogue of radial velocities, mostly of stars of spectral classes N and R, was pub- 
lished by the writer eight years ago in Mount Wilson Contributions, No. 525.! This cata- 
logue contains 191 stars, of which 187 were observed at Mount Wilson and 4 at Santiago, 
Chile, by the D. O. Mills Expedition. The distribution in type is: class N, 147 stars; 
class R, 40; class S, 3; peculiar, 1. 

Some of these velocities were measured with the Hartmann spectrocomparator, but 
most of them with a micrometer measuring engine. The velocities for the standards used 
on the Hartmann spectrocomparator were based upon velocities obtained by J. H. 
Moore? for class-N stars and upon those obtained by W. C. Rufus? for class-R stars. For 
the micrometer, lines were chosen and their wave lengths determined on the basis of 
Moore’s velocities for class N and Rufus’ velocities for class R. 

In Mt. W. Contr. No. 525 it was concluded that Moore’s velocities for class-N stars 
and those derived both by the Hartmann comparator and by micrometric measurement 
at Mount Wilson needed a systematic correction of —7 km/sec, and this correction was 
applied. The velocities of class-R stars appear as measured originally. 

Observations of class-N and class-R stars were resumed at Mount Wilson in 1935 and 
have been continued when conditions permitted during the next eight years. Improved 
photographic plates, especially the Eastman 103 F-emulsion for the yellow-red region, 
have greatly facilitated observing the fainter stars and have made it possible to observe 
many of the brighter stars with high dispersion. 

Spectrograms of the faint stars of class N having a dispersion of 65 A/mm were ob- 
tained in the yellow-red region with the Cassegrain plane-grating spectrograph.* For 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washington, No. 689. 
1Ap. J., 82, 202, 1935. 
2 Lick Obs. Bull., 10, 79, 1923. 
3 Pub. Astr. Obs. U. Michigan, 2, 45, 1916. 
4 Mt. W. Contr., No. 432; Ap. J., 74, 188, 1931. 
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class R, spectrograms with a dispersion of about 70 A/mm at Hy for the brighter stars 
and 100 A/mm for the fainter ones were obtained in the blue-violet region with one- 
prism spectrographs. Spectrograms (dispersion either 20, 10, or 5.9 A/mm) of 30 bright 
class-N and 10 bright class-R stars have been obtained with the coudé plane-grating spec- 
trograph. 

Additional velocities are now available for 106 of the stars in Mt. W. Contr. No. 525 
and for 45 class-N and 48 class-R stars not previously measured. These 93 stars, observed 
for the first time, were largely selected from the Henry Draper Catalogue and from the 
lists of discoveries made by objective prism at Mount Wilson® and at the Dearborn Ob- 
servatory.® The improved values for velocities in Contribution No. 525, together with the 
velocities for the 93 other stars, make it desirable to issue a completely new catalogue. 

The catalogue.—All stars of classes N and R with measured radial velocities are listed 
in Table 1. Of the 283 stars included, all have been observed at Mount Wilson except 
4 class-N stars observed at Santiago, Chile. The designation MSB indicates the Mount 
Wilson lists of stars discovered on objective-prism plates;> Lee indicates stars similarly 
discovered at the Dearborn Observatory.® 

The magnitudes have been taken from various sources. Visual magnitudes for non- 
variable stars are from the Henry Draper and Bonner Durchmusterung catalogues, from 
the Dearborn Observatory Annals, and, lastly, if no value is available, from estimates by 
the writer at the finder or at the guiding telescope. These estimates are in parentheses. 
The magnitudes for the maxima and minima of variable stars are, for the most part, from 
Schneller’s “Catalogue and Ephemieris of Variable Stars.’’? Photographic magnitudes are 
underscored. 

All the stars are placed in one or the other of the Harvard spectral classes N or R. 
Those in class N have not, however, been subdivided with the letters a, b, and cc, asin © 
the Harvard classification, or with the numbers 1-9, as in C. D. Shane’s® classification. 
Only the broad classification N has been used, with e appended when the spectrum has 
bright lines and p when the spectrum is peculiar. 

Keenan and Morgan® have published a new system of classification for stars with car- 
bon bands. Two of the data which they use as criteria are given in Table 1; the eleventh 
column contains estimates of the strength of the D lines on an arbitrary scale of 1-9, and 
the twelfth column, estimates of the strength of the Swan band-head at \ 6192 on a scale 
of 1-5. By means of these two columns at least a rough conversion of class-N stars to the 
classification of Keenan and Morgan can be made, if desired. 

The subdivisions of class R in the tenth column of Table 1 resemble closely those of 
the Henry Draper Catalogue, except that the progression in the strength of the isotope 
bands of CC" has here relatively more influence. These class-R stars may be put into 
Keenan and Morgan’s classification by using the nearly linear relation which it bears to 
the Henry Draper classification.’ 

Although their classification may provide a better temperature sequence than the 
Harvard classification, it is still unsatisfactory in some respects; it seems not to be con- 
sistent with the observed relatively low intensity of the D lines in the spectrum of some 
extremely red class-N stars, like VX Andromedae, for example. The contents of the tenth, 
eleventh, and twelfth columns may therefore serve as a compromise classification until 
such difficulties are removed. 

The mean absorption-line velocities in Table 1 have been formed with appropriate 


5 Pub. A.S.P., 45, 306, 1933; 54, 107, 1942. 

6 Dearborn Obs. Annals, 4, Part 16, 1940. 

7 Kleinere Veriff. Sternwarte Berlin-Babelsberg, No. 21, 1939. 
8 Lick Obs. Bull., 13, 123, 1928. 

° Ap. J., 94, 501, 1941. 10 [bid., Fig. 1. 
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weight for the velocities of class-R stars by Rufus and of class-N stars by Moore. A sys- 
tematic correction of —7 km/sec was first applied to Moore’s values. 

The letters accompanying the velocities express the writer’s judgment of the probable 
accuracy. The uncertainties are about as follows: a, 1 km/sec; b, 2-3 km/sec; c, 4-5 km/ 
sec; d, 6-8 km/sec. Velocities with the letter e usually depend upon a single low-dis- 
persion spectrogram taken primarily to check the spectral classification. It was not prac- 
ticable to obtain more spectrograms. The underscored velocities depend largely or 
wholly upon those obtained from coudé spectrograms. The emission-line velocities in the 
last column of Table 1 are derived for the most part from Ha. 

The radial velocities from absorption lines—The velocities of class-N and class-R stars 
in Table 1 are on the same basis as those in Mt. W. Contr. No. 525. Before using these 
velocities it had to be ascertained whether they were free from systematic errors, which 
was done by comparing the velocities obtained with the coudé spectrograph with those 
for the same stars with the smaller dispersions. 

The velocities from coudé spectrograms could be so used because experience has shown 
that this spectrograph is practically free from instrumental errors and because the follow- 
ing considerations assured accurate wave lengths for the absorption lines used for veloc- 
ity measurement. Many atomic lines with well-known wave lengths are well isolated for 
measurement in the blue-violet region of classes N and R and in the yellow-red region of 
early class R. Lines of comparable quality in the yellow-red region of the spectrum of 
late class-R and class-N stars belong to the band spectrum of cyanogen. Wave lengths for 
these lines were obtained by referring them to the lines in the neon comparison spectrum 
on a coudé spectrogram of the carbon arc. Comparison showed that the weighted mean 
difference, coudé minus low dispersion, for 26 suitably observed class-N stars is —0.4 + 
0.6 km/sec and for 7 class-R stars +0.0 + 0.8 km/sec. There appears, therefore, no 
reason to apply systematic corrections to the velocities of class-R or class-N stars obtained 
with lower dispersion. 

Velocities from bright hydrogen lines.—Thirty-four stars of class N with hydrogen emis- 
sion lines appear in Table 1. All are known variable stars with the exception of six of the 
Mount Wilson objective-prism discoveries and one of Lee’s stars. These also may turn out 
to be variable when adequately observed. At any rate, the evidence is strong that emis- 
sion lines of hydrogen are favored by variability. Sufficient spectroscopic observations 
have been made at all phases of a number of these variables to show that the emission 
lines are at their maximum strength near maximum light and practically vanish for a 
considerable interval before and after light-minimum. 

For the most part, the emission-line velocity has been derived from Ha alone. This 
velocity is invariably smaller algebraically than the absorption-line velocity. The weight- 
ed mean difference (emission minus absorption) for these 34 stars is — 20.4 + 0.9 km/sec. 
The correlation of this difference with the length of the period of light-variation is weak, 
there being, perhaps, a slight tendency to larger negative values for the longer periods. 

Galactic distribution.—Stars of class N are plotted by galactic co-ordinates in Figure 1, 
which shows marked galactic concentration; only about 10 per cent of the stars are more 
than 30° from the galactic plane, and these are predominantly bright stars whose real dis- 
tances from the plane cannot be great. In contrast, stars of magnitude 9.0 or fainter lie 
for the most part within 10° of the galactic plane. The galactic concentration of stars of 
class R (Fig. 2) is not marked. 

The mean absolute magnitudes.—Studies of the motions of distant celestial objects have 
amply demonstrated the differential effects of galactic rotation. This rotation produces 
a double sine-curve of velocity variation for the 360° of galactic longitude. Previous val- 
ues of the semi-amplitude of this curve for a distance of 1000 parsecs range from 15 to 
about 20 km/sec. A recent discussion by R. E. Wilson" of the best available data gives 
17.7 km/sec/kiloparsec. 


1 Mt, W. Contr., No. 631; Ap. J., 92, 170, 1939. 


; 
) 
) 


TG+ N | | 98 TH+ | 6 GSW | ded AV 
96+ N ous S$ -| 62T | Bf 6e+ | 6°22 | 
: bu eee" | 9T-| Bet | 9 get | B°TT B0TLE 
Pt- ON $°6 | os+ | 6°8 9 SH | ded 
It- T+] | & TS+ | O°F F GSH | 026,0S+ 
P 9 N L-| | 8% | GSW | ded OV 
99 + 9 N ose? | Tt | 6°¢e © GSW | 
s N 0°6 | 6°9 60T 6T | sree wed 
B6- 6-| | OS | 6°02 osete ted K 
ag+ L an 8-| PIT] 42 | L°9 TSs6T 
Q 9 N O° | 8T- | | Be+ | Tree puy 
9T- ou vee et 90T | Sst | GSN | Jed ZA 
eee OT+ Pu eee ost es 6 e 2°oe STI9T 
a+ Q ON e°et| t°s 99- | T8T | Se | 
agp. | ¢ 9 N tte @-| SOT] se TSt+ | 
BOLT- | fe) ou o°6 | | Set | LH | 9°6 TAV A 
eee o¢c- ou eee eeenee 6°6 8 e 66 8s L°3¢ 9S90T eeeree 
6S- 9 N 96 PT LG+ | T | SBD ma 
9 an 9Le 8°s8 | 3°9 | | | 6L88 TOS 
Ott L N ts | 9c- | 86 | 9°OT T9SL osd Z 
| ¢ ON 1(g*0T)| - | Tt es+ | O'S T | Gh ASH 
8 N eee o°s | 8°4 | 98 get | 2°22 puy dv 
eee 2 OF- ¢ 9 N eee eeee (3°6) e 48 asw 
N s°6 | 8I- | 38 6 | O° FT 9bST puy XA 
2 N | g° et- | 38 bP | SST LS 
P 8°ST | OST | ol + | | | 049 yO | xn 
* *sqy | | ex (sfeq) cut | | Tt | *toeq 
*qgoedg 
(oes /uy)*TeA | £47 potied Tus ey O06T aH 


N GNV SASSVIO dO SUVLS dO SHILIOOTSA 


++ 


a 
+++ 


‘+ 


+ 


ON 


OFNOM NODONH CONDON 


NNO OW 


ajo 


DARrO 


2 Tajo 
YHODW 
eee ee @ 


amy nn 
any AY 
0$ 


LLTT 
£862 
wep av 
t40 Ta 
ZOPT 


wep ZZ 
Nd 
any A 

LT 


any dv 
TS6 oS T+ 
neL 
neL OL 


LY 
any WO 
any 


8T 9eT 
aay An 
SSLolTt 
tata AS 


any XL 
M 


L 


eee err eee er ewme eee of emeee eeme eeeee eevee i 
eee on “meee oe eeeee eee ee 
J 
| 
| 
| 
ea 
‘ 
DOONM HONWONO 000 
Qn AYONH FAN FOMNDR Ag ] 
onor oD ANN WM oe *HAQ> 4 
faq : 13338 
+ 
i 
i] 
| 
i 


eee © ° 86 LT+ | 28T | 06 + THO 
© 96+ N Tet| | TOT | OT Set | 6°FT O9TLS wep Nd 
oe | | \ \(S°OT) | | LT- | 6°0T | GSR 
| O°6 | 6°L 62+ | SIT | 69+ | 6°OT L9T9S weg ny 
© 9 N 86T | ST LT- | GSH | 998ToLT- 
eee 99+ 2 9 N eee eeece (0°6) eo tog 2e 6T- 3°6 2s asn 
eee qa 8 + ou eee eeee (9° TT) L + Zet 8s re) + 
P 9 N | 6°9 T-| | 9F | M 
Q OT+ 9 N see | 94T | 819 + | OZ9Ts AY 
eee © LS+ g N | 6829+ | 82 GSN | TSPT.9+ 
eee ¢ N ove at 9+] SLT | STL + | 6°0S GSW | UOW OG 
oL- dou oes get | 96 | | HT | 6%60S amy MO 
P set N v°OT| 2°6 S8Tt | 42 - | 
PO | N eee 4°6 @-| 2 L- | | 
Lit N ot T+| 64T | + | | O9 
© N eee 1(0°6) 9-| O6T | 6% 2T- | 
N +OTT 8°s | 9°8 | IST | Te+ | wep MA 
*sav | wn (sfeq) | xen T toed 

ao 
(ces/my)*tea | £4 opng Tus BD 006T qu 


eee 
eee ° 8 + 2 
eee ° 9oT- e 
| eget-| + 
eee Pp 
| 
eee B St+ 
| 
ee- | Per! 
eee . 
eee 
eee as + ° 
eee Qa 
eee ° St+ 
eee ss an 
| 
eee 
eee Q 26- 
| 
eee ¢ + 
eee 6S+ 
eee 8 . 
OT- 
eee Bs + 
eee es + 
eee Got ° 


od 

> 

+ 
e 


eee Ti+ 

eee te+ 

eee ¢ 
a9 + 
eee q 
eee qa 
LT+ 
eee 
eee 
eee ° 
eee + . 


Oro ee 


eee ee ee 


*o we 


orwno wo oO 


ze 


OOADO 


200 


. 


al 


pe 


ao 


° e 

~ 


| 
et 


ornn 


wo 
a 


aA 


Tt 


ot 


| 
| 
| 
| 


= 
= 
= 
= 
= 
4 
= 


|b | 9 N at | ec] POBTLT | 39S XU 
N OT | _6 2 b+ | Teso4t | Ab 
20 cu O°2T | O°TT b- | | BS | S6POLT | 438 SS 
eee q . Ou eee eeee L + ¢ 9 + Ltt 
° 6T- : ou 8°6 | 0°6 | 62 ST- | 9°ST | 
eee q 8 + Py eee eevee ot t’s 
eee ° S L N eee eeee (0° OT) ss 2°R @eeeee asn 
| by | 92 6 + | O°F 46099T | 
cot-| ° cu 60T TT]}| 0° OT OT- | Tee | 08 | ¥°8S 
eee oot- Pu eee (OT) OT+ 9 Oot t+ Stt ee, 
: eu | os+ | 09 6 9+ | | 
9¢- N | 9°6 9-| gee 62 Se- | 8*°0F TTST9T | xs 
6T+ ¢ 9 N | | LE BT- | | 23S ZS 
98+ 9 N 2° TT]| 0°6 4L-| 9T2 | SE TH | SOZO9T | IL 
9 N ot | 9TS | SS | 4°02 ASH | 
P6°- 9 N €-| | 9T | B°ST 69 ASH | 
eee ° ¢ N eee 0°6 S°L 8 + LT? Tt 2e- 2°Re as 
0 toe S°OT | O°4 ce+ | tee | St et- | estert | udo A 
eee L + N eee eeee 2°6 + 90¢ Le T° 
eee 2 . Pu eee e@eee L°OT 43 9°L LLLSPT 
9 den | 8°4 | OF 9) OS+ | S°T 9T | SLSPPT | 
BSti-| 2 9 oN O°ST | os+ | o¢ 2S | | A 
wg sqv o| BN (sfeq) | vu 

*qgo0eds 
(ces /uy)*toq | £9 poyaed epng dH 


| 
| | 
| | 


son | 8 N | | 6 +] | Go got | | 
eee 9° eee (6) 82 re Lt+ 6°02 T ee] 
| i on 4S} | | T°9 | 9+] | GE | oses6t | 349 
6+ | | 8 on | grat | 6 Let | | XM 
N tte | | 9 OSS | Te- | LELZ6T | LY 
| | 8 on sey | 4°68 | | B+! SH | 98 Bet | B°6 | 340 su 
| 1 a T-| 9¢ | | 9°99 | 4S GSH 
eee ¢ ¢ N eee eeee 2°6 4 o 82 
eee N eee L°8 Lt PT 65 + 
eee 2 ou eee 8°OT 6c L oe so96st 
eee oe eee eeee 2°6 ss 6 get 8° Cr 
oe | Bot | | ose| ve et- | 9°08 go9gst | 23s mn 
| | N ost | | | GE | Se see | | 340 LL 
N | | ve | ec | 89 asm | 349 AH 
eee a ¢ 9 N eee grat 6- tt of 9 + toby 
N | BOT | O'S | Stt+ | 9% | OG | oc asm | 349 
| | 4 on | Tet | | t-/| OST FLSA 
| | N | O84, | | | SB out | esseet | eid xn 
eee Lb- e e 2u eee vse es L°at 
| | | 9 on | | | LE | Le “e+ | asm | 0 
eee St 9 N eee 2°h St- Ss 9T- e3608T 
eee eee (0°6) + 22 St To+ 8°2t 99 asn 
@ee Z 2 cu eee (0°6) 9 + Le+ 8°6 asn 
su cee | | | O8- | & | Tes SSS6LT | 
| pest | | | BOL) 9-| 2 s¢ tT - | SST64t 
eee ° e ou eee 9°OT Lee 92 Lt- 61 
| | 9 | | | ase] os g - | Tes | TOV A 
| iss | | 9°6 | H+! HT | HT | O° ooge4t | thy an 
18 | | | HE] BT B- | ALESLT | 99S 
|» | | | BSS] T B- | | 398 S 
ou | 09 - | | | LySSyBT | SSTSLT | 398 AU 


| ] 
| 
| 
i 
j 
| 


eee eee eeee (s°6) oe 64 92 T°6t oe asn 
ou | 89- It | te- | ¢*9 susete | 
ane Q cree 1(5°6) | | 98 Sb+ | 2°9 | | 
LT- N L°s | - | Te | 6°TS | AL 
eee q cu eee eeee ( LS- Ze oc L e 
9u ove Tat | ¢- | Sf | CL9STS | XL 
Q LT+ on (6) | oS | Of LI+ | O°0F 
eee ° 8t- N eee (¢°8) 6 8S+ SL asn 
09- RPuz- | 2 6 On OFF | 0°6 8T- 9g | 3 | S*t | Bea zy 
° cu Le- Lb | bf OS+ | L°6S | 
9 N errs 1(9°6) | 8 - £9 |S Sb+ | ges6oe | 
Le- oN Gut 28. | Lb | | L°TS gzssoe | 3eq xu 
2 cu vids S°IT | S°OT | | OS+ | 9°TS etssoe | 
eee OT+ N O°L 0°9 | set | 04g902 | Zed sa 
q ST- N 69 | 82 09+ | LLLSOS | °°°°°° 
OT- 9 Lee L- | 8S TH+ | | 2% 
| at-| a4 oN st]| | te- | at | 9¢0- | te] | aby 
ope TSt | | o |S | tTs9s6t | B49 sa 
eee ° wu eee eeee 8°6 L e L°6e POOLET 
2 8 en | 8°9 e+ oS | | 62 GSW | A 
*sqy | | wn (sfheq) | | T |*to0ed ash 

*g0edg zo 
(ces | potszedg Tus eH an 


| | 

| 


*SeuTI 09 UMOCUH ST ON 
pus uvuesy TI LOL9Y UMoU SUH ZA CH 


‘out “Bt ‘*Qnd HOTT GH 
* TT 
ZM UT se IT TI LOLOY Xm GH 
=e 
“9861 “LES “Ting “ar y6T 
‘Lee *TIng “08S AS ‘wo"s 
°ON UT Suotm pey eA *Ss980T dH 
‘OOT ‘8ST ‘*and HOTT s 
dH 
“TOT ‘Spteapy pue “anv ZX °m6°9S 
pus og ‘I SW ‘[ III 0] Jo ‘any AN GH 
Area fem 00Te, “doy *966TS dH 
‘ce ‘°qnd °8Q0 HOTT dH 


T OL SHLON 


ee ou eee 6°6 ez ¢ oe o°Ls eserves 
8L- Q 9F- 4 On ne (2°OT)| - | SZ 9S+ | O° tS | SL 
eee 3B e cu 8°s ese 99 os Ss + O° zecces 
2 ecu sec Se- SL | St | | PUY LS 


| 

| 


JO Xopur Aq [enstA 07 


156 


157 


158 ROSCOE F. SANFORD 


Within 30° of the galactic plane there are 171 stars of class N whose velocities can be 
studied for evidence of galactic rotation. Only HD 189711 was omitted on account of its 
outstandingly large velocity. These velocities, freed from solar motion (apex a 270°, 
6 +28°, and V» = 20 km/sec), have been gathered into thirteen velocity means (Ta- 
ble 2). The last column shows the mean radial velocity freed from solar motion and re- 
duced to mean apparent visual magnitude 8.3. The effect of differential galactic rotation 
is evident from the plot of the data of Table 2 in Figure 3. 


TABLE 2 
NORMAL VELOCITIES—CLASS N 
No. of Stars ms i b 
8.9 32 0 +17.8 
8.3 49 +12.0 
9.0 72 +1 5.7 
9.2 132 +4 + 4.0 
8.5 150 0 + 5.4 
8.9 170 0 +13.0 
Seer 8.5 190 0 +19.1 
7.9 218 +5 + 8.4 
7.8 328 +2 + 3.8 
KM. 
SEC ° 
+20 
240 200 160 120 80 40 0 320 280 


GALACTIC LONGITUDE 


Fic. 3.—The mean velocities of class-N stars in Table 3 are plotted as dots and connected with 
straight lines. The mean velocities of class-R stars in Table 4 are plotted as circles; they have been cor- 
rected for K and dV (p. 159) and reduced to visual magnitude 9.8 (general absorption included). 


The mean distance, 7, for class-N stars of magnitude 8.3 has been derived from the 
data of Table 2 by a least-squares solution of the conditional equation for differential 
galactic rotation: 


7A sin 2(1 — Io) cos? b = p, 


wherein A = 17.7 km/sec/kiloparsec, J) = 325°, and / and fare from Table 2. The values 


of 6 in the fourth column differ so little from zero that cos? b was assumed to be unity. 
The mean distance found for class-N stars of apparent visual magnitude 8.3 is 1.10 + 
0.08 kiloparsecs. The evaluation of the mean absolute magnitude depends to some extent 
upon what value is adopted for general absorption per kiloparsec. If 0.35 visual magni- 
tude is assumed, the mean absolute visual magnitude becomes —2.3 + 0.2. This may be 
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compared with —1.88 derived by R. E. Wilson” from proper motions and the radial 
velocities of the earlier catalogue. 

Table 3 gives thirteen mean velocities derived from 62 stars of class R in the same way 
as the data for class N in Table 2. HD 187216 and HD 209621 were omitted because of 
their large velocities. p is for mean apparent visual magnitude 9.0. Inspection shows evi- 
dence of differential galactic rotation, but a lack of symmetry with respect to the axis of 
zero velocity suggested including the well-known K-term and a correction to Vo, the 


TABLE 3 
NORMAL VELOCITIES—CLASS R 

No. of Stars m i b ry 
yee 8.8 30 + 6 + 2.8 
9.4 52 +10 + 7.4 
8.8 73 —10 — 7.7 
9.5 90 — 3 —14.0 
8.5 127 + 6 —14.3 
9.6 139 0 —11.0 
9.3 163 +18 + 2.6 
9.2 185 +15 + 4.3 
214 +10 + 2.6 
8.4 315 +12 + 3.9 
8.4 348 —12 + 2.9 


adopted solar motion, in making the solution for differential galactic rotation. A condi- 
tional equation of the form 


K + dV cos (i — 20°) + 7A sin 2(/ — 1) = p 


was used (6 assumed to be 0). Least-squares solution gave 
K = —3.5 +1.3 km/sec, ~ 
dV = +1.9 +1.6 km/sec, 
r= 0.67 + 0.10 kiloparsecs. 


This mean distance is for a mean apparent visual magnitude 9.0. Absolute visual mag- 
nitude —0.4 + 0.4 for class-R stars follows from the above mean distance and mean 
magnitude, if we use a general absorption of 0.35 mag. per kiloparsec. Wilson’s!* value 
is M = —0.50 + 0.20. 

The value of K might be interpreted as evidence of a systematic error in the velocities 
if it were not that intercomparison of velocities from spectrograms of high and low dis- 
persion has shown no such error. No obvious explanation is at hand. 

The value dV is only slightly larger than its probable error; hence, the use of a value 
differing significantly from the +20 km/sec actually applied for solar motion receives 
scant justification. 

The values of p in Table 3 have been corrected for K and dV, derived above, and then 
made to correspond to a mean visual magnitude of 9.8 (corrected for general absorption). 


12 Mt. W. Contr., No. 631; Ap. J., 92, 188, 1940. 
13 Mt. W. Contr., No. 618; Ap. J., 90, 492, 1939. 
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The resultant values for p, plotted as circles in Figure 3, are in fair accord with the values 
for class-N stars of mean visual magnitude 7.9 (8.3 uncorrected). This should be true, 
since class-R stars, which are fainter than class-N stars by the difference of their absolute 
magnitudes (1.9), should be at the same distance and therefore give similar curves for 


differential galactic rotation. 
Table 4 contains those stars whose velocities after correction for solar motion are nu- 


merically larger than 100 km/sec. 


TABLE 4 
STARS OF HIGH VELOCITY 

8.3-9.0 RO — 163 122 —44 

Lee 11.3 R5 —137 229 +71 

9.2-9.6 R6p —135 74 +79 

9.6 R3 —129 85 +27 

8.7 Nl —167 17 —12 

8.8 R3 —379 47 —27 

9.9 RO —132 63 —64 

TABLE 5* 
INTERSTELLAR D LINES 
Equiv. 
Star m Spect. —— Width M 
0.5(D2+D1) 

8.8 R3 0 km/sec O.15A 2.2 
Var. RO — 3 .49 —1.7 
Var. N —14 .28 —1.9 
| 8.7 N —10 .56 —-1.7 
8.8 R3 —10 0.27 +0.4 


* Part of the data of this table appears in Pub. A.S.P., 54, 257, 1942. 
t V Ari: mag 8.3-9.0. 
tV Aql: mag 6.7-8.2. 


Keenan and Morgan" have noted that the four stars marked with asterisks in Table 4 
are characterized by very strong g-bands and weak atomic lines in the blue region. Lee 
107, HD 187216, and HD 224959 are similarly characterized. All seven have large nega- 
tive velocities. It may be worth noting that one of these stars (HD 5223) has the largest 
proper motion known for any star of class R or N. This suggests that this group of stars 
may be of lower absolute magnitude than the average class-R star. 

Interstellar sodium lines.—Four class-R and two class-N stars have velocities large 
enough to permit the interstellar D lines to be separated from the stellar D lines. It has 
been possible to derive both the interstellar velocities and the total absorption of the inter- 
stellar D lines from the coudé spectrograms (Table 5). The values for 0.5 (D2 + D1) 


M Ap. J., 94, 509, 1941. 
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have been used to read off distances for the stars of Table 5 by the curve of equivalent 
width versus distance given by P. W. Merrill for class-B stars. 

The absolute magnitude of HD 5223 from the interstellar D lines is +2.2, noteworthy 
as much fainter than the magnitudes obtained for the other stars of the table. This is 


- TABLE 6 
ABSOLUTE MAGNITUDES 


Method Authority Class R Class N 


Proper motions and radial velocities 


Radial velocities (Table 1).......... R. F. Sanford —0.4 —2.3 
Total absorption of interstellar Dlines| R. F. Sanford —0.8 (0.0) —1.8 


quite consistent with this star’s abnormally large proper motion, already alluded to. The 
mean of the absolute magnitudes of the four class-R stars is 0.0, or —0.8 if HD 5223 
is omitted. The two class-N stars give a mean absolute magnitude of —1.8. 

Absolute magnitudes for stars of classes R and N by three independent methods are 
brought together in Table 6. The class-N stars are unquestionably more than a magni- 
tude brighter than class-R stars. 


1% Mt. W. Contr., No. 569; Ap. J., 86, 35, 1937. 
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ABSTRACT 


The results of wave-length measurements on eight spectrograms of Comet Whipple-Fedtké-Tevzadze 
(1942g) are presented. The spectrograms, obtained with relatively high spectral purity (projected mono- 
chromatic image of the slit at the plate ~ 0.5 angstroms for \ 3880), cover the ordinary photographic re- 
gion. In addition to a general consideration of the spectrum of the comet, several points of particular in- 
terest are discussed. These include, first, the resolution of the R branch of the A 38830,0 CN band into 
its individual lines and, second, results of increased resolution in the cases of several features of the \ 4050 
group of CH; bands, in particular the duplicity of the main maximum at A 4051. Third, reasons are given 
for believing that the emission at \ 4752, which for other comets has been attributed to the CC" isotopic 
carbon molecule, is not in the present instance due to this molecule. Fourth, the interesting but possibly 
fortuitous wave-length coincidences between certain weak cometary emissions and bands, ascribed tenta- 
tively to the molecule SiO2, are noted. 


INTRODUCTION 


Comet Whipple-Fedtké-Tevzadze (1942g) was a relatively bright object, being be- 
tween the third and fourth magnitude for a few days.! At the time of its maximum bril- 
liance it was unusually well situated for observation from northern latitudes. During late 
February and early March, 1943, some time after perihelion, which occurred on February 
6, 1943, the comet unexpectedly increased in brightness; and during this period twelve 
spectrograms of its head were obtained at Victoria. In securing these observations two 
purposes were kept in mind. The first was to obtain spectrograms of the heretofore neg- 
lected visible spectral region to the red of \ 5000. The second was to photograph the spec- 
trum of the comet in the ordinary photographic region, below \ 5000, under conditions 
of the highest practicable spectral purity. Some of the twelve plates were obtained with 
panchromatic emulsions and were useful in both these spectral regions. 

The study of the visual region of the spectrum of Comet 1942g has been completed and 
was published? in a joint paper with P. Swings and R. Minkowski, who had secured data 
on the same region of cometary spectra at the McDonald and the Mount Wilson ob- 
servatories. The main point of interest in that study was the tentative identification of 
the molecule NH2 as the source of several of the strongest visible emission bands. The 
present article presents the results of measurements on eight spectrograms covering the 
ordinary photographic region. While most of the emission features of which wave-length 
measurements are given have been measured as such in the spectra of previous comets, 
our increased spectral purity has made possible the resolution of others among them into 
components not detected up to the present time. 


THE OBSERVATIONAL MATERIAL 


The plates were obtained with the single-prism form of the universal spectrograph. 
Two different cameras were employed, one having a focal ratio {/3 and the other f/5. 
All the data concerning the spectrograms thought to be of interest or value are assem- 
bled in Table 1, and they require little comment. In addition to the high spectral purity 
mentioned earlier, one further point may be noted, namely, that during the period over 


1 See, e.g., J. Ashbrook, Pop. Astr., 51, 362, 1943; N. T. Bobrovnikoff, Pop. Astr., 51, 481, 1943. 
2 P. Swings, A. McKellar, and R. Minkowski, Ap. J., 98, 142, 1943. 
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which the plates were taken the comet was at the relatively large heliocentric distance, 
r = 1.4astronomical units. Most slit spectrograms of comets have been secured at values 
of r of 1 astronomical unit or less. Further reference will be made to this later. The dis- 
tance of Comet 1942g from the earth in late February and early March, 1943, was ~ 0.5 
astronomical unit. 

THE MEASUREMENTS 


The results of the measurements on the eight plates are given in detail in Table 2. 
After a brief description of the spectral feature measured, the wave length as found from 


TABLE 1 
DATA ON SPECTROGRAMS OF COMET WHIPPLE-FEDTKE-TEVZADZE (1942g) 
RaDIAL SPECTRAL 
VELOCITY SPECTROGRAPHIC Purity, 
ere OF DISPERSION PROJECTED Pu QuALiy 
D ComMET Surt-Wiptr* os 
ATE (U.T.) HELIOCENTRIC GRAPHIC 
ESPECT GRAM 
TO At At At At 
EarTH 3880 4740 3880 | 4740 
Astronomical 
units km/sec 
1943, Feb. 24.535.... 1.38 +7.9 |57 A/mm)130 A/mm|0.85 A/1.94 A|Eastman |Good 
103a0 
25-4385. 1.38 8.0 |32 72 0.57 {1.30 |Eastman |Good 
103a0 
27 .486.... 1.39 8.3 130 0.85 |1.94 |Eastman /Fair 
103F 
28 .481.... 1.39 8.5 |32 72 1.14 |2.60 |Eastman /Good 
103F 
Mar. 1.468.... 1.40 8.6 |57 130 0.50 |1.14 |Cramer Excellent 
Hi-speed 
3.467 1.40 8.9 |57 130 0.41 (0.94 |Cramer Fair 
Hi-speed 
7.469.... 1.42 9.5 157 130 0.57 |1.30 |Eastman |Weak 
103a0 
9.460.... 1.43 +9.8 |57 130 0.50 |1.14 |Eastman /|Good 
103a0 


* In expressing ‘‘spectral purity’’ in terms of the width in angstrom units of a monochromatic image of the slit at the plate, 
as we have done, it is assumed that the spectral lines have rectangular profiles. Therefore we are not adhering to the strict defi- 
nition of ‘‘spectral purity’’ as given originally by Schuster. The latter definition makes it a function, among other things, of the 
aperture of the collimator. However, our usage indicates the maximum resolution to be expected and would seem to provide 
important information on resolving-power not shown by presenting the dispersion alone. 


each plate is listed and is succeeded by the mean value from all the plates. Then follows 
the identification of the emitting molecule where this is known, and, finally, the labora- 
tory wave length of the feature is given. Directly below the date in the second to the 
ninth columns, the quality of the spectrogram is noted. 

Results of the individual wave-length measurements, obtained from the micrometer 
settings by the Hartmann dispersion formula, were calculated to the hundredth of an 
angstrom unit. Many of the emissions were so diffuse that the last figure has little sig- 
nificance. Thus, in obtaining the mean cometary wave length, in cases where either the 
number of individual measurements was only one or two or the feature was diffuse, the 
wave length has generally been given only to the tenth of an angstrom. In deriving the 
mean value, the wave lengths found from the plate of March 1.468, in view of its excel- 
lent quality, were given double the weight of each of the others. 
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Three other points should be mentioned. First, the measured wave lengths were cor- 
rected to take account of the radial velocity of the comet with respect to the earth. Sec- 
ond, wave lengths of all band edges or band heads measured were further corrected for 
the effect of the lack of infinitely high spectral purity. That is, from the wave length of 
each red edge of an emission band measured, there was subtracted an amount equal to 
half the projected spectrographic slit width at that wave length, while this quantity was 
added to the wave lengths of the violet edges of the bands measured. Third, in the few 
cases in which an intensity minimum is identified as an actual Fraunhofer absorption line 
in the scattered solar spectrum (e.g., \ 4045) its wave length has been corrected also for 
the effect of the radial velocity of the comet with respect to the sun. 


DISCUSSION 


A reproduction of the spectrum of Comet 1942g is shown in Plate V. The most note- 
worthy feature of general interest regarding the spectrum of this comet, in the ordinary 
photographic region, is the relatively great strength of the \ 4050 group of bands due to 
CHz. This group and the omnipresent and always strong \ 3883 0, 0 sequence of CN 
bands are the most outstanding of the emissions and are of approximately the same total 
intensity. The usually very strong \ 4737 1, 0 sequence of the Swan system of C; is less 
intense than either the CN or the CH2 bands. The A 4315 band of CH is faint but definite- 
ly measurable. These exceptional intensity ratios are doubtless a consequence of the 
considerable heliocentric distance of the comet, r = 1.4 astronomical units. At such a 
distance the breakdown of “‘parent”’ polyatomic molecules to form CN, CH, C2, etc., by 
photodissociation under the influence of sunlight, has not proceeded to such an extent as 
at closer approach of the comet to the sun. Therefore, as pointed out previously,® the 
emission bands due to polyatomic molecules should be, and apparently are actually found 
to be, more intense with respect to those of diatomic molecules at the greater heliocentric 
distances. Much the same general observations as above and similar conclusions are con- 
tained in a brief note by Swings and Struve‘ on the spectrum of Comet 1942g. 


THE BANDS OF CN 


On every plate the \ 3883 0, 0 CN band exhibited the complex structure which has 
been the subject of study over the last few years.** This is evident from the measure- 
ments given in Table 2. One spectrogram, that of March 1, taken with projected slit 
width at \ 3880 of only 0.5 angstrom, showed unusually fine definition. The appearance 
of the band is shown in Plate V. On this plate the rotational structure of the R branch 
(nonhead-forming branch) of the \ 3883 band is just barely resolved. Reference to Table 
2 will show that nearly all the individual lines from the short-wave-length extremity at 
the line R(12) up to the band origin have been measured. The intervals between succes- 
sive lines are about 0.7-0.8 angstrom unit. In so far as the writer is aware, such resolution 
has not been obtained up to the present on the CN bands in the spectrum of any comet. 
Use will be made of this plate in a spectrophotometric study to determine the rotational 
distribution of the cometary CN molecules. 

It will be noted in Table 2 that for the sharp lines in the CN band from R(12) to R(1), 
while the agreement between the accurate laboratory and the cometary wave lengths, 
determined mainly from one low-dispersion plate, is fairly good, there appears to be a 
small but definitely systematic difference between them. The mean difference for ten 
lines is \ (comet) — A (lab.) = +0.09 angstrom unit. Unfortunately, there are no other 
spectral features with well-determined laboratory wave lengths and sharp enough in the 


3 See, e.g., P. Swings, Rev. Mod. Phys., 14, 190, 1942; G. Herzberg, Rev. Mod. Phys., 14, 195, 1942. 
4 P. Swings and O. Struve, Pub. A.S.P., 55, 150, 1943. 

5 P. Swings, Lick Obs. Bull., 19, 131, 1941. 

6A. McKellar, Rev. Mod. Phys., 14, 179, 1942. 
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cometary spectrum to test the reality of this difference, except, perhaps, the \ 4300 and 
d 4304 lines of CH. The measurements of these are not definitive enough to permit a de- 
cision on this point. Under the circumstances of observation the existence of a systematic 
error of this amount cannot be excluded. If, however, the shift of wave length can be at- 
tributed to a physical cause which may be referred to the comet, the only obvious ex- 
planation would appear to be an average motion of recession (toward the tail of the 
comet), the line-of-sight component amounting to some 7 km/sec. While no definite con- 
clusion can be drawn from the present observations, the matter is of sufficient interest to 
warrant further attention in the future when greater spectrographic power is available for 
studies of this kind. 

The matters of the over-all wave-length range of the CN band and its general intensity 
distribution are worthy of note. In the present case the short-wave-length limit of the 
band is at \ 3866.2, the R(12) line, and the long-wave-length limit at \ 3881.5, the P(15) 
line. Thus the band appears to be developed to such an extent that lines originating from 
the rotational levels up to K’ = 14 appear. The maximum of intensity in the P (head- 
forming) branch occurs at \ 3880.16, between the lines P(10) and P(11). For comet Cun- 
ningham (1940c), it was found® that, for the heliocentric distances r = 0.92 and r = 0.54, 
lines from levels up to about K’ = 25 were present and the P-branch maximum occurred 
at P(13) for r = 0.92 and P(20) for r = 0.54. Dufay’ had earlier found from examination 
of spectra of several comets that the K” corresponding to the CN band maximum ad- 
vanced from K” = 6 to K” = 14 as r decreased from 1.8 to 0.5 astronomical units. The 
general intensity profile of the CN band in the spectrum of Comet 1942g is thus in con- 
formity with these recent suggestions regarding the increase in ‘‘rotational temperature” 
(the term ‘“‘temperature”’ being used only in the sense referring to the relative populations 
of the molecular rotational states) with decreasing heliocentric distance. Also, it yields an- 
other point in establishing the quantitative relationship between 7 and the K”’ corre- 
sponding to maximum intensity. The role in the determination of this intensity maximum 
played by the radial velocity of the comet with respect to the sun and discussed in refer- 
ences 5 and 6, must not be ignored. 

The 0,1 sequence of the violet CNV system at \ 4216 appeared on the plates with mod- 
erate intensity. In view of the smaller spectrographic dispersion compared with d 3880 
and the lower intensity of this sequence, nothing worth while is gained by a detailed dis- 
cussion of it. 


THE A 4050 GROUP OF CH» BANDS 


On the spectrogram of March 1, 1943, the resolution in the region of the \ 4050 group 
of bands was also very good. Plate V includes a reproduction of the group from this spec- 
trogram. The projected slit width at \ 4050 was 0.6Q,angstrom. These bands have recent- 
ly been shown by Herzberg* to be due to the triatomic molecule CH2. A comprehensive 
résumé, including previous summaries of the wave lengths of this group of bands in the 
spectra of many comets, is given by Swings, Elvey, and Babcock.° 

Measurements on the present plates have revealed certain structure not previously 
noted. Of perhaps most interest is the duplicity of the strongest emission feature at 
d 4051. On four of the eight plates it has been found to be double, with components at 
d 4050.43 and \ 4051.75, the component of longer wave length being perhaps a little 
over twice as intense as the other. On the plate of March 1, the region from A 4013 to 
d 4020 appears to contain four intensity maxima. Also, the duplicity of the emission at 
d 4068, heretofore suspected, is confirmed, components being measured at \ 4066.17 and 


d 4069.99. 


7J. Dufay, C.R., 206, 1948, 1938. 
8 Loc. cit.; and Ap. J., 96, 314, 1942. 
®P. Swings, C. T. Elvey, and H. W. Babcock, Ap. J., 94, 342, 1941. 
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According to current ideas on the resonance or fluorescence origin of cometary emis- 
sion bands, the observed intensity distributions over these bands would be expected to 
show certain definite differences from those encountered in laboratory studies. These dif- 
ferences would be a consequence of the intensity distribution over the region of the solar 
spectrum absorbed in the course of the production of the bands. In the case of the CH2 
bands, the three very strong Fe absorption lines at AA 4045, 4063, and 4072 in the solar 
spectrum would doubtless have considerable effect in this way. They would cause intensi- 
ty minima in the CH2 bands at these three wave lengths and possibly other wave lengths, 
depending on the as yet unknown detailed band structure. The CH2 bands shown in 
Plate V appear to exhibit this effect, particularly for \ 4072, where the emission centered 
at 4074 seems to be cut off sharply on its violet edge. The effect mentioned has, of 
course, to do only with the cometary emissions, not with that part of the cometary spec- 
trum which is simply reflected or scattered solar radiation and which exhibits the normal 
Fraunhofer lines. 

Unfortunately, the laboratory investigation and analysis of the \ 4050 group of CH2 
bands, from high-dispersion spectrograms, at present being conducted by Dr. Herzberg, 
is not yet completed. While it could be suggested that the main sharp maxima are prob- 
ably Q branches, yet the detailed interpretation of the structural features of the band 
discussed above and of the intensity profile to be expected under cometary conditions of 
excitation cannot profitably be undertaken now. The band structure should be re-ex- 
amined, however, when the complete laboratory data on the system become available. 


THE SWAN BANDS OF C2 


As already noted, the Swan bands of C; appeared with somewhat lower intensity than 
usual. On the ordinary plates the 1, 0 sequence, \ 4737, was readily measured, and a trace 
of the 2, 0 sequence (the 4, 2 band) was found to be present. On the panchromatic plates 
the 0, 0 sequence, 5165, and the 0, 1 sequence, A 5635, occurred. 

One point of interest may be noted, namely, that on the two best plates a sharp emis- 
sion feature was measured near \ 4752. This is the wave length of the head of the 1,0 
Swan band due to CC", the diatomic carbon molecule made up of two atoms of the car- 
bon isotope of atomic weight 13. This band is the analogue of the \ 4737 Swan band due 
to ordinary CPC”, 

The relative abundance of C” to C' has been reliably measured" in terrestrial sam- 
ples of carbon as about 100 to 1. Also, a spectrographic analysis of some meteoritic sam- 
ples" has indicated that the relative abundance of C™* in these is not sensibly greater than 
this value. If the abundance ratio C” to C'* is a, the relative abundances of the isotopic 
molecular species CC”: C?C: CC are as a?: 2a:1, that is, for terrestrial or meteoritic 
carbon, approximately 10,000: 200: 1. 

Bobrovnikoff” and later Swings," on the basis of their own measurements on the spec- 
tra of Comets Halley (1910) and Brooks (1911) and also Wright’s measurements on the 
spectrum of Comet Brooks, concluded that the isotopic bands of CC and C¥C at 
\ 4745 and \ 4752 were present, hence that the carbon isotope, C'’, was a constituent of 
the heads of these comets. While this result is not specifically questioned, yet in the case 
of Comet 1942g the relatively low intensity of the main, C°C, Swan system and also the 
absence of the C’C' band head at \ 4745 lead us to conclude that the feature measured 
at \ 4752 is not due to the isotopic carbon molecule C'*C'’. Perhaps one would think that 
the fact that CC" is not a strictly homonuclear molecule and so probably has a different 


10 F, A. Jenkins and L. S. Ornstein, Kon. Akad. v. Wetensch., Amsterdam, 35, 1212, 1932; A. O. Nier 
and E. A. Gulbransen, J. Amer. Chem. Soc., 61, 697, 1939. 

1 F, A. Jenkins and A. S. King, Pub. A.S.P., 48, 323, 1936. 

?2,.N. T. Bobrovnikoff, Pub. A.S.P., 42, 117, 1930;and Pub. Lick Obs., 17, 441, 1930. 

13 P. Swings, Ap. J., 95, 270, 1941; and M.N., 103, 86, 1943. 


| 
| 
) 
) 
) 
> 


170 ANDREW McKELLAR 


rotational distribution of molecules than C?C” and CC’ could be invoked to explain, 
partially at least, the absence of \ 4745. But the origin of this band is at about 4 4738, the 
band head occurs at about K”’ = 17, and the rotational constant B” = 1.6, so this argu- 
ment is not strong. In any case, the matter of relative abundances inclines us to the view 
that in the present instance the emission at \ 4752 is due to some as yet unidentified 
emitter. Swings! has pointed out that for Comet Cunningham (1940c), on spectrograms 
on which the very intense Swan bands were strongly overexposed, there was no trace of 
the isotopic bands. He advances the suggestion that the relative abundance of C” and 
C's may be different in various comets. If our present view should be erroneous and the 
\ 4752 emission in the spectrum of Comet 1942g should actually be due to C'%C", this 
fact, together with the data on Comet Cunningham, would provide very strong evidence 
in favor of Swings’s suggestion. 

An additional point regarding the measured wave lengths of the heads of the 1,0, 2,1, 
and 3,2 bands of C2 is evident from Table 2. These wave lengths are consistently 
about 1 angstrom unit greater than the laboratory values. Since, as was earlier stated, 
care was taken to correct the measured wave lengths by an amount equal to one-half the 
projected slit width, this cannot be the cause. No reason for this difference is apparent. 


THE CH BANDS 


While relatively weak, the \ 4315 0, 0 band of the A, 7IT transition was definitely pres- 
ent. The Q-branch maximum at \ 4313 and the sharp lines of the R branches at \ 4304 
and \ 4300 were measured. No trace of the less intense band due to the *2, *II transition, 
with heads at A 3872 and \ 3889, was detected. 


OTHER EMISSIONS 


Certain emission features occurring between the 0,1 CN band at \ 4216 and the CH 
band near \ 4300 may now be dealt with. Swings!’ has recently identified three of these 
at AX 4231, 4238.5, and 4254.4 as the R(O) + R(1), the Q(1) + Q(2), and the P(3) lines 
of the 0, 0 band of the violet CH* system. Incidentally, the first of these, suspected by 
Swings to be double on one of his plates, has been found to be so for Comet 1942g. 

While examining wave-length material on the spectra of polyatomic molecules in con- 
nection with the spectrum of Comet 1942g, it was noted that the strongest bands of a sys- 
tem tentatively ascribed to SiO,"* occurred in this same region. Furthermore, it appeared 
that very good wave-length coincidences occurred between the most intense SiO2(?) 
bands and certain of the cometary emissions. The state of affairs in this respect is pre- 
sented in Table 3, which gives five columns containing data on the spectra of comets, fol- 
lowed by columns giving the laboratory wave lengths of CH* band lines and of SiO2(?) 
band heads and intensity maxima. The letters “V,” “R,” and “M” following the wave 
length in the last column indicate the heads of bands degraded to the violet, those de- 
graded to the red, and headless bands, respectively. All the numbers in parentheses give 
relative intensity estimates except in the CH* column, where, of course, they refer to 
the rotational quantum numbers used to designate the band line. When considering the 
wave-length region covered in Table 3, it should be called to mind that the two mole- 
cules responsible for all the emission bands known in the spectra of cometary tails have 
bands in this region. There are CO* band heads at AX 4248.9, 4252.4, 4272.0, and 4274.3; 
and there are V,*+ band heads at \ 4236.5 and \ 4278.1. 

Strangely, the wave lengths of the CH* band lines and the SiO2(?) bands are very 
similar, and the general agreement between cometary and laboratory values for the two 
molecular spectra is about the same. Favoring the conclusion that some of the cometary 


4R. C. Pankhurst, Proc. Phys. Soc. (London), 52, 707, 1940; R. W. B. Pearse and A. G. Gaydon, © 


The Identification of Molecular Spectra, pp. 179-80, London: Chapman & Hall, Ltd., 1941; L. H. Woods, 
Phys. Rev., 63, 426, 1943. 
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emissions may be due to SiO,(?) are (1) the close wave-length coincidences and (2) the 
fact that the bands occur in the comet spectrum at the heliocentric distance r = 1.4. 
Opposing this conclusion is the following greater number of points: (1) The lines at \ 4228 
and \ 4232 are definitely sharp and so are more likely to be the single lines of the CH* 
band. (2) SiO. (quartz) in solid form is a most refractory compound and, if present in 
molecular form, must have been formed by dissociation of some parent polyatomic mole- 
cule. No obvious parent comes to mind. In this connection carbon is also very nonvolatile, 
and yet C2 bands appear in the coolest hydrocarbon flames and in cometary spectra, 


TABLE 3 
COMETARY EMISSION BANDS IN THE AA 4230-4280 REGION 
Comet Unidentified 
Whipple- Comet Comet Comet Nuclear oe Wave Lengths 
Fedtké- Cunningham Brooks Halley | Emissions, of Bands 
Tevzadze (1940c)* (1911¢)t (1910)t | Summary by “all Due to 
(1942g) Baldet Si02(2) 
r=14 AU. r~0.6 r~0.5 r~0.7 (1926)§ 
1 
4239 .4(0) 4240 .4(1—2n)} 4238.5(0)(long |....... 4238(1) 4239.37 @Q(2) 4240 M(10) 
line, tail) 
4254. 3(1n) 4254.5(0+)(long}.......| 4255.8(2) | 4254.38 P(3) 4254.4R (7) 
line, tail) 
4265(0) 4264.5} 4264.2(2) | 4262.13 P(4) 4262 .8R (6) 
CO*?) 


*P. Swings, Ap. J., 95, 270, 1942; Contr. Lick Obs., Ser. II, No. 3. 


t W. H. Wright, Lick Obs. Bull., 


7, 8, 1912. 


tN. T. Bobrovnikoff, Lick Obs. Pub., 17, 445, 1931. 
§ F. Baldet, Ann. obs. Astr. phys. Paris, 7, 58, 1926. 
|| A. E. Douglas and G. Herzberg, Canad. J. Research, 20, 71, 1942. 
© R.C. Pankhurst, Proc. Phys 
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where, of course, the molecular carbon results from dissociation. Also, Bobrovnikoff 
states’ that for Comet 1882 I at r = 0.87 before perihelion and for Comet 1914 V at 
r = 1.21, the sodium D lines were observed in emission, despite the fact that the boiling- 
point of sodium is 877° C and, at r = 1, the temperature of a black body is only around 
0° C. (3) If SiOz were present, one might expect to find bands due to SiH and SiN and 
perhaps CO, in the spectrum of the same comet. A careful search in the present case has 
not revealed these bands. (4) We do not know whether or not the bands produced in the 
laboratory and ascribed to Si02(?) arise from the ground state of the molecule, a condi- 


1° N. T. Bobrovnikoff, Rev. Mod. Phys., 14, 171, 1942. 
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tion necessary for their appearance by resonance and virtually necessary for appearance 
by fluorescence under cometary conditions. 

It must be remembered that the emissions we have been considering in this section 
are among the weaker ones occurring in cometary spectra. In view of this and also the evi- 
dence against the significance of the possible S702(?) origin of the bands, as opposed to 
the less weighty favorable evidence, the writer is inclined to the belief that the SiO2(?) 
bands are probably not present. Nevertheless, it was thought worth while to discuss the 
matter as we have done, for it cannot be said to be settled with certainty. Two future de- 
velopments—first, better observational data on bright comets and, second, further labo- 
ratory work on the Si0,(?) bands, clarifying their molecular origin and including their 
detailed vibrational and rotational analysis—should enable a definite conclusion to be 
reached on the matter. 

There still remain in Table 2 about a dozen unidentified emission features. Study of all 
available laboratory data on molecular spectra has failed to reveal any reasonable identi- 
fication for them. The recent studies resulting in evidence for the occurrence of C Hz bands 
and probably N Hz bands in cometary spectra would seem to enhance the possibility that 
an appreciable portion of these unidentified features may subsequently be found to arise 
from polyatomic rather than from diatomic molecules. 

In conclusion, the need for and the value of good future spectrographic observations of 
comets should be emphasized. In the dicussion of the measurements which are the basis 
of the present paper there have arisen several problems, including possibly significant 
wave-length shifts of the sharp CN band lines, the identification of the weak \ 4752 emis- 
sion with the C¥8C'’ band, and the identification of emissions in the \A 4230-4280 region, 
which cannot be settled with the data now at hand. Only by the study of better cometary 
spectrograms—and there is no substitute for these—can conclusive answers to the above 
questions be advanced. When the occurrence of a bright comet makes it possible, slit 
spectrograms secured with the highest practicable combination of dispersion and spec- 
tral purity should be obtained. Relatively few observatories have the telescopic and 
spectrographic equipment necessary for obtaining such spectra, and, when such factors 
as the transitory nature of comets and the diffuseness of even their nuclei are considered, 
it is understandable why so few high-purity cometary spectrograms have been secured. 
In urging that steps be taken to remedy this deficiency, by those who may be in a position 
to do so when a new bright comet should appear, we may state with certainty that such 
cometary spectrograms, when obtained, will be of great value. Not only will they help us 
to arrive at decisions on the above-mentioned matters, but also they will doubtless bring 
forward, and aid in the solution of, other and important problems in cometary physics. 
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ON THE SPECTRUM OF COMET 1940c 
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ABSTRACT 


The spectrum of the comet was measured on two spectrograms obtained at the Perkins Observatory 
on December 23 and 24, 1940. The structure of the main cyanogen and carbon bands is discussed. The 
fine structure of the cyanogen band (0, 0) can be identified with certain rotational lines in this band. The 
presence of sharp lines in the carbon band (1, 0) is established, but their connection with the band is 
doubtful. A list of other lines and bands measured on the plates is given. 


OBSERVATIONAL MATERIAL 


Owing to unfavorable weather, only two measurable spectrograms of Comet 1940c 
(Cunningham) could be obtained between December 1, 1940, and January 4, 1941. They 
were made with the two-prism spectrograph and the 350-mm camera, attached to the 
69-inch reflector. The plates were Eastman 103-0. The width of the slit was 0.35 mm, and 
the length 10 mm. The projected slit width was 1.46 A for \ 3883 and 3.80 A for 4737. 
Because of the shortness of exposure, most cometary bands appear as sharp lines with a 
width much smaller than the above figures would indicate. 

In view of a thorough investigation of the spectrum of this comet made elsewhere,! it 
might seem superfluous to give the results of my measurement. However, since our spec- 
trograms have a much greater dispersion (16.1 A/mm at A 3883 and 41.7 A/mm at 
\ 4737) than ordinarily used in comet spectroscopy, many details either only suspected 
by the former investigators or not noticed at all are quite clear. The difference in the 
intensity of the various bands in cometary spectra is so great that the faintest bands can 
be seen only when the strongest bands are overexposed and their structure concealed. 
While many of the fainter bands cannot be measured on our spectrograms, the structure 
of the stronger ones is striking. 

The details for the spectrograms are as shown in the accompanying tabulation. 


No. Mid.-Exp. U.T. Exp. Time 
1940 Dec. 23.993 60" 


The comet in the Cassegrain focus appeared as a strongly condensed globular body 
of 10” in diameter, surrounded with a faint coma. Visual observers of that period state 
that the nucleus was stellar. Therefore, what appeared to be the nucleus on the slit of the 
spectrograph was in reality a dense spherical envelope around the stellar nucleus. 

The slit was oriented along right ascension, and the nucleus of the comet held in the 
center of the slit. The guiding, in spite of the low altitude of the comet, was satisfactory. 
The measurement of the ee shows that the nucleus was displaced from the 
geometrical center of the slit b ry only 0.11 mm or 2"8 in angular measure. The nuclear 
line on the spectrograms is 11°5 wide, corresponding to the spherical envelope around 
the nucleus mentioned above. 

1 Swings, Elvey, and Babcock, Ap. J., 94, 320, 1941; Contr. McDonald Obs., No. 38. I adhere to the 
designation “‘1940c’’ used by these writers, although the correct designation may be 1940d, depending on 
whether Comet 1940a (Kulin) is a comet or an asteroid (Pop. Astr., 48, 490, 1940). 
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As the position angle of the radius vector on December 25.0 was 206.5°, the slit made 
an angle of 63.5° with the radius vector on the west side. The west side of the comet was 
therefore nearer to the sun than the east side. As a consequence of this fact, the struc- 
ture of some bands shows a curious difference between the west and the east sides. 

The plates were measured in both directions and reduced in four sections by the 
Hartmann formula. The main difficulty of measurement was the curvature of the lines 
introduced by the long slit. The wire of the micrometer was set on the tips of the iron 
comparison lines, and a correction of —0.0124 mm was applied to the measured abscissa 
of the comet lines. No evidence was found of any appreciable systematic error due to this 
difficulty in measurement. From the internal agreement of the measures of the cometary 
lines, as well as from the comparison lines not used in the derivation of the Hartmann 
afr the reliability of measurement is estimated at 0.03 A for \ 3883 and 0.06 A for 
4737. 

All bands except the strongest parts of the CN (0, 0) and (1, 1) bands appear on the 
spectrograms as sharp lines, considerably sharper than comparison lines of comparable 
intensity. Therefore, the wire of the micrometer was set to bisect the cometary lines 
rather than on their red edge. 

The measured wave lengths in the spectrum of the comet were corrected for the Dop- 
pler effect resulting from the relative motion of the comet and the earth. For December 
23 this was — 27.7 km/sec and for December 24, — 28.7 km/sec. Although the hour angle 
of the comet was over five hours, the Doppler effect due to the rotation of the earth was 
only —0.3 km/sec. 

WAVE LENGTHS 


The finally adopted wave lengths are given in Table 1. Column 7 contains estimated 
intensities on an arbitrary scale, Column SEB the wave lengths reported by Swings, El- 
vey, and Babcock for the same comet, and Column H the wave lengths found by me in 
the spectrum of Halley’s Comet.’ 


THE STRUCTURE OF CN BANDS 


The fine structure of the (0, 0) band of CN at \ 3883 appears to be identical on both 
plates, the average difference being only +0.01 A. It is evident that no modification in 
the CN bands took place between these two dates. 

The main portion of the (0, 0) band is overexposed, but the three lines in it can be seen 
without difficulty, especially farther away from the nucleus. The individual lines in the 
main portion of the (1, 1) band are not distinguishable, and we may speak only of the 
maxima of intensity. 

As to the identification of the fine structure of the bands, I do not think the present ob- 
servational material—or, in fact, any material so far published—is adequate to settle the 
question. Taking as an illustration the five lines in the (0, 0) band measured on both 
plates, we have the comparison between the cometary and laboratory wave lengths* for 
December 24 shown in Table 2. 

. Formally, such an identification must be regarded as highly satisfactory; yet the as- 
signment of rotational lines is in every case different from that given by Swings‘ in his 
adopted cometary wave lengths for the (0, 0) sequence. 

It would seem that the measured features are actually individual rotational lines and 
not groups of lines, but I think this conclusion would be erroneous. The lines in the band 
are rather symmetrically placed, and it is possible to select a combination of rotational 
lines resulting in the actually observed wave length of the line. For instance, taking six 


2 Lick Obs. Pub., 17, 433 ff., 1931. 
3 Laboratory wave lengths from Uhler and Patterson, A p.J., 42, 434, 1915. 


4 Lick Obs. Bull., 19, 131, 1941. 
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TABLE 1 
WAVE LENGTHS IN COMET 1940c 
Dec. 24 Dec. 23 i SEB H Identification Description 
1 57.7 55.2 CN R(i,1) Violet edge of a band 
1 62.5 63.1 CN P(1,1) Red edge of a band 
2\ 67.5 CN P(1, 1) Max. int. in a band 
69.0 CN P(1, 1) Max. int. in a band 
2 69.9 69.8 CN R(0,0) Max. int. in a band 
73.69 3 CN R(O,0) Very sharp line 
76.92 4 76.8 76.3 CN P(0,0) Very sharp line 
79.89 7\ 80.2 79.5 CN P(0,0) Line in the main band 
80.96 10/ CN P(0,0) Line in the main band 
RG a aa 82.32 10 82.0 81.8 CN P(0,0) Line in the main band 
1 43.6 42.5 Diffuse line 
1 97.0 97.4 CN R(O,1) Faint and diffuse line 
$206.15 06.45 1 CN R(0,1) Faint and diffuse line 
1) 10.2 CN P(0,1) Very sharp short line 
11.86 1f 11.4 CN P(0,1) Very sharp short line 
3} 14.7 1325 CN P(0,1) Long line 
14.92 3 15.0 CN P(0,1) Long line 
1 CH R(0,0) Very faint line 
91.67 2 CH R(0,0) Diffuse line 
97.17 2 97.3 96.0 R(0,0) Diffuse line 
99.66 2 CH R(0,0) Sharp line 
03.35 5 ves CH R(0,0) Sharp line 
13.13 8 CH Q(0,0) Very sharp on violet edge 
29.3 29.5 CH P(0,0) Sharp line 
1 44 42.8 CH P(0,0) Sharp line 
1 48 50.8 CH P(0,0) Sharp line 
4961 60.88 59.2 Very hazy line 
63.91 65.0 C, P(4,2) Sharp line 
70.42 71.6 C: P(3,1) Sharp line 
4380.68......... 80.74 80.4 C, P(2,0) Sharp line 
78.26 77.9 C, P(5,4) Sharp line 
84.58 84.7 C, P(4,3) Sharp line 
97.31 97.4 G , Sharp line 
4714.98........ 14.81 15.2 C, P(2,1) Sharp line 
45.8 C2C8 (0,0) Diffuse line 
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rotational lines between \ 3882.58 and d 3882.10, we would get the mean wave length 
3882.35, which is within 0.02 A of the wave length measured in the comet. The width of 
this combination of lines would be only 0.48 A. A much higher dispersion and a nar- 
rower slit are evidently necessary to settle the question of identification. 

One important consequence follows from Table 2, as well as from the data on other 
comets. There are no large radial velocities involved in the comet. Although CN has been 
shown to be an important constituent of halos and jets, constantly forming and expand- 
ing in many comets, the velocities of expansion are always small, seldom exceeding 1 
km/sec. Such a velocity of approach would result in the displacement of the CN lines in 
the (0, 0) band of only 0.013 A to the violet. The sharp character of the CN, as well as of 


TABLE 2 
THE (0,0) BAND OF CN 
Comet i aboratory C-L Ident itt Comet Laboratory C-L Identifi- 
cation cation 
3882.32 A | +0.05A | P(18) || 3876.87A..... 3876.84A | +0.03A | P(3) 
3881.00 — .03 PUTS) 3874.00 —0.05 R(1) 
3879.58 —0.01 P(9) 
TABLE 3 
WIDTH OF CORE OF THE +1 SEQUENCE OF C, 
WipTtH WipTtH 
BanD BanD 
Dec. 23 Dec. 24 Dec. 23 Dec. 24 
2.14 2.02 


other lines, speaks against any great turbulence in the atmosphere of the comet. The ab- 
sence of large velocities of ejection from the nucleus, up to 60 km/sec. required by the 
mechanical theory of cometary forms, is one of the most puzzling facts. 

The distribution of intensity along the CN lines on either side of the continuous spec- 
trum is uniform, the sunward side not being favored. The bands stretch across the whole 
spectrum more than 1’ in width, and undoubtedly extend much farther. 

Plate VI shows the structure ‘of the CN bands (the 0 sequence) on December 24. As 
fainter details are difficult to reproduce, the lower figure (b) shows the lines retouched on 
the original enlargement. 

The (0, 1) CN band at \ 4216 has the structure analogous to that of the main band; 
but, because of the lower dispersion and faintness of the bands, the measures on the two 


dates are not in such good agreement as for the (0, 0) band. 


THE STRUCTURE OF THE SWAN BANDS 
Each plate has a set of five usual bands of the +1 sequence of C:. The cores of these 
bands corresponding to the continuous spectrum are very sharply defined. Their width is 


given in Table 3. 
The projected slit width was 3.80 A for the first band and 3.62 A for the last. There is 
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no significant difference in the width of these bands, except for the last one on December 
23. It is strikingly wider and more diffuse on that date than on the following. It appears 
that some modification in the structure of the Swan bands took place within the one 
day which separated these two plates. 

This conclusion is strengthened by the measurement of the wave lengths of the band. 
Table 4 gives the difference between the measured wave lengths of the P-heads in the 
comet and those in the laboratory.® 

The measurement was repeated but with practically the same result. It appears that 
on December 23 the Swan bands were displaced to the violet by 0.22 A as compared with 
December 24. If interpreted as a Doppler effect, this corresponds to a velocity of ap- 


TABLE 4 
COMET—LABORATORY 
Band Dec. 23 Dec. 24 Band Dec. 23 Dec. 24 
—0.13 A —0.03A |] 4,3......... —0.21A +0.05A 
— .39 — .22 —0.34 —0.01 
Mean....| —0.27A —0.05A 
TABLE 5 
DISTANCE OF SATELLITES FROM MAIN BANDS 
Bands Dec. 23 | Dec. 24 


proach of 14 km/sec. Such an interpretation is in contradiction to the results obtained in 
the measurement of the CN bands; and in all probability we have to deal here with some 
modification in the structure of the P-heads of the Swan bands. The average difference 
between the wave lengths of the +-2 sequence of the Swan bands on the above dates is 
—0.25 A, or 17 km/sec, an almost identical figure; but the measures are scattered, as 
these bands are much fainter than the +1 sequence. 

The cores of the first three bands are accompanied on the violet side by satellites visi- 
ble only in the nuclear line. The wave lengths of these satellites cannot be determined 
precisely, because they are somewhat diffuse; but the measures referred to the centers of 
the satellites are in good agreement for the two dates (Table 5). 

Furthermore, between the bands (2, 1) and (3, 2) there are two very sharp emission 
lines extending across the whole width of the spectrum and especially well developed on 
the western, or sunward, side of the nucleus. Their weighted mean wave lengths are 
4703.52 and 4707.12, corresponding to the “‘subhead”’ measured in Halley’s and other 
comets as one band at \ 4705.7. There is no question that this is a double band; and, from 
the absence of a similar structure between other Swan bands, we must assume that it is 
independent of the Swan spectrum. The same remark applies to the sharp line measured 


at \ 4689.87. 
On the other hand, there is no trace of a band at \ 4727, also measured in many comets 


as a “subhead” of \ 4737 (1, 0). 
5 W. Jevons, Report on Band Spectra, p. 11, Cambridge University Press, 1932. 
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Finally, on the plate of December 24 the isotopic band of the C” — C8 molecule is 
definitely present as a faint diffuse line. The spectrogram of December 23 has no trace of 
it, although, judging from the intensity of \ 4737, it should be present if it was of the 
same intensity as on December 24. The computed distance between the (1, 0) Swan band 
and the isotopic band is 7.55 A, the actually measured distance is 7.22 A. The discrepan- 
cy of 0.33 A is not significant in view of the faintness and diffuseness of the isotopic 
band. 

The McDonald observers could find no evidence of the carbon isotope band on their 
spectrograms; but my results are in agreement with the previously established fact? of 
the considerable variation in intensity of the isotope band with respect to the Swan 
bands. In Halley’s Comet one plate would have the isotope band, while another with the 
Swan bands of the same intensity would have no trace of it. 

The present results do not support Swings’s hypothesis that the diffuse features at 
d 4723 and d 4706 are the (2, 1) and (3, 2) transitions of the isotope band. As has been 


TABLE 6 
RELATIVE INTENSITIES OF SWAN BANDS 


ee 9 10 8 5 5 2500° 
7 10 8 3 2 3000 
ey ae 9 10 10 1 2 2000 


mentioned before, the band at 4706 is double, with a structure entirely different from 
the (1, 0) isotope band at \ 4745. The presumed (2,1) band is missing in Comet 1940c. 
On the other hand, in Comet 1941c, the feature at \ 4723 is very prominent, but neither 
the (1, 0) nor the (3, 2) band is present (unpublished results). 

A few words can be said about the distribution of intensities along the bands in the di- 
rection perpendicular to the continuous spectrum. It has been shown’ that the spectrum 
of various formations around the nucleus may differ considerably from the general spec- 
trum of the head. 

Such observations require, however, the co-operation of several observers unless the 
features are prominent enough to be seen on the slit of the spectrograph. The head of 
Comet 1940c in the 4-inch finder of the 69-inch telescope appeared structureless. A di- 
rect photograph obtained with a 6-inch Brashear lens attached to the tube of the 69-inch 
reflector is overexposed and shows the head as a round body of 4’ diameter surrounded 
by a spherical envelope 6’ in diameter. It is impossible to say exactly which features of the 
comet were on the slit during the exposure. 

However, the Swan spectrum shows some interesting differences in the distribution of 
intensity along the bands, as can be seen in Plate VI (Figs. c and d). The last two bands 
(4, 3) and (5, 4) can be traced much farther on the west side than on the east side of the 
nucleus, but the first three bands are stronger on the east side. The estimated intensities 
of the bands on both plates in the continuous spectrum (core) and on the west and east 
side (about 16” from the nucleus) are shown in Table 6. 

The last column contains the temperature taken from Wurm’s tables® for the distribu- 


6 The theoretical explanation of why C!*? — C'* bands are sharp and C!* — C8 bands are diffuse is given 
by Swings, M.N., 103, 86, 1943. 
7 Observational evidence for this is discussed in my article in Rev. Mod. Phys., 14, 164, 1942. 


8 Zs. f. Ap., 5, 260, 1932. The tables were calculated for the case of absorption, but they will suffice 
to give a general idea of temperatures involved. 
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tion of intensities in the vibrational transitions of this molecule. In fact, we can find the © 
distribution of intensity quite different from that given in Table 6 if we consider various 
distances from the continuous spectrum. At the very edge of the spectrum on the west 
side, for instance, all the bands appear to be of about the same intensity. This shows the 
artificiality of the temperature concept as applied to comets. 

The position of the core of the bands with respect to the continuous spectrum is also 
not quite the same for different bands, as can be easily seen on the reproduction. Table 7 
gives the results of measurement of the spectrogram of December 24. 

At the geocentric distance of the comet, 0.801 astronomical unit, 1’ corresponds to 
570 km. It is seen that the center of the core of the (1, 0) band was displaced by more than 


TABLE 7 
CORES OF SWAN BANDS 
Band Length Displace- Band Length Displace- 
ment ment 
16"1 POW 12"1 0"8 W 


1000 km west, or sunward, with reference to the nucleus, while the center of the (2,1) 
core was displaced by 300 km away from the sun. Of course, the line of the continuous 
spectrum is not very definite, but the positions of the cores were measured in reference to 
the sharp edge of the comparison spectrum. Therefore, the relative displacement of the 
two bands by 1300 km is real. What bearing this has on the physical processes in the 
head of the comet is hard to say, except that the vibrational transitions corresponding 
to these bands are localized in regions not symmetrical in respect to the nucleus. Even 
more striking differences in the distribution of intensity of CH bands with respect to the 
continuous spectrum seem to be present, but the bands are too faint for any precise 
measurement. The CN bands, on the other hand, do not show any difference in this re- 
spect. 

P The rest of the lines and bands given in Table 1 call for no especial comment except 
\ 4313, the strongest band in the CH system A?A — x*x. It is impossible to identify it! 
with the whole complex of eight Q branches stretching from \ 4311.3 to \ 4314.2. The 
band was measured very consistently at \ 4313.21 and \ 4313.13 on the two plates as a 
sharp line of 1.32 A in width. It is especially sharply defined on the violet side, and there is 
no possibility of any emission to the violet of \ 4312.6. 
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ABSTRACT 
In this paper the equation of transfer ; 
dI +1 


is solved by expanding / in terms of spherical harmonics in the form 
I(r, = Ti(r)Pi(u) 


It is shown how successively higher approximations to J can be obtained by retaining more and more 
terms in the foregoing expansion; thus, the first, second, and third approximations are obtained by re- 
taining the first three, five, and seven terms, respectively, in the expansion. The first three approxima- 
tions to /9(7) are explicitly obtained, and the corresponding laws of darkening derived. 

The exact function describing the law of darkening and given by Hopf’s theory is also numerically 
evaluated. A comparison of this exact law with those derived on our second and third approximations 
indicates that the method described in this paper is sufficiently rapidly convergent. 


1. Introduction.—As is well known, the solution of the equation of transfer 


(1) 
dr 
where 
I sin ddd, (2) 


has occupied a central position in the theory of stellar atmospheres and has been the sub- 
ject of investigations by Schwarzschild, Milne, Eddington, Jeans, Hopf, Bronstein, and 
others. An account of these investigations and the successive methods of approximation, 
as developed particularly by Milne and Eddington, will be found in the various treatises 
on the subject by Eddington, Milne, and Unsdld.! However, the methods described by 
these latter writers, while adequate for most purposes, are not sufficiently systematic, in 
the sense that no general procedure is outlined which would enable one to obtain solu- 
tions of increasing accuracy by retaining, for example, more and more terms in an expan- 
sion of some sort for 7. More recently, however, an attempt in this direction has been 
made by L. Gratton? in a paper of considerable interest. Gratton’s basic idea is to seek a 
solution of equation (1) of the form 


I(r, 8) = (1) Pi (cos 8), (3) 


l=1 


1A. S. Eddington, The Internal Constitution of the Stars, chap. xii, Cambridge, England, 1926; E. A. 
Milne, Handb. d. Ap., 3, No. 1, 114-126, Berlin: Springer, 1930; A. Unséld, Physik der Sternatmosphiren, 
pp. 89-104, Berlin: Springer, 1938. 


2 Societa astronomica italiana, 10, 309-325, 1937. 
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in which the P,’s denote the various Legendre polynomials. In other words, we expand J 
in terms of spherical harmonics with the expectation that by retaining a sufficient num- 
ber of terms in the expansion we shall be able to obtain solutions with any desired degree 
of accuracy. This is, of course, an entirely sound procedure; but, unfortunately, Grat- 
ton’s specific treatment of the problem following this idea is vitiated by a number of 
errors and oversights which require corrections. Accordingly, it has been thought worth 
while to re-examine the whole problem de novo and obtain what might properly be de- 
scribed as the first, second, and third approximations to the solution of equation (1). 
We shall, moreover, compare the predictions of these approximate solutions with regard 
to those features (e.g., the law of darkening) for which the Hopf-Bronstein* theory pro- 
vides exact information. 
2. The outline of the method.—Denoting cos 3 by u, we can re-write equation (1) in the 
form 
(4) 


where J, is the first term in the expansion 


Substituting the foregoing expansion for 7 in equation (4) and remembering that the 
P7’s satisfy the recursion formula 


pq Pin (6) 
we find 
— 2/+1 ‘dt <= 
Equating the coefficients of the various Legendre polynomials in equation (7), we ob- 
tain 
and 
1 
(l=0). (9) 


Equation (9) leads at once to the integral 
I, = constant = (say) , (10) 


which clearly insures the constancy of the net integrated flux and determines the effec- 
tive temperature of the star. Another integral of the equations is obtained by considering 
equation (8) for the case / = 1. For, according to this equation, for / = 1 we have 


dh, 


3 E. Hopf, Mathematical Problems of Radiative Equilibrium (Cambridge Mathematical Tract No. 31), 
Cambridge, England, 1934. 
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or, since J; is a constant (eq. [10]), 
Ion 3Fr+a, (12) 


where a is an arbitrary constant at our disposal. This integral is readily verified to be the 
same as what is more generally known as the ‘“‘K-integral.’’* 
Before we proceed to outline our method of approximation we shall first formulate the 
boundary conditions of the problem. 
For the problem of the radiative equilibrium of a stellar atmorphere the boundary 
condition at 7 = 0 is 
I(0,u) =0 for —1<y<0, (13) 


since no external radiation is assumed to be incident on the star. The implication of the 


boundary condition (13) for the 7;’s as r > 0 can be derived in the following manner. 
From equation (5) we conclude that for 7 > 0 


2 +1 ; 
=f. I(r, Pi (pu) du. (14) 


Passing now to the limit ts = 0 and remembering that, according to equation (13), 
I (0, #4) vanishes for all negative values of u, we obtain 


2 1 
=f (0, Pi(u) dy. (15) 
On the other hand, for t = 0 (cf. eq. [5]) 
1 (0, = (0) Pm (x) (16) 
m=0 


Hence, combining equations (15) and (16), we must have 


= (0) f Pru) Pm (a) da. (17) 


m=0 


Thus the boundary condition (1) is equivalent to the infinite set of linear relations (17) 
among the J;’s as r > 0. 

It is now apparent that we cannot satisfy the entire set of relations (17) when a solu- 
tion for J is sought by retaining only a finite number of terms in the expansion (5). But 
it is clear that the more the number of the relations (17) we are able to satisfy, the more, 
in general, we may expect the derived solution to be accurate. On the other hand, only 
as many of the relations (17) can be satisfied as there are disposal constants in the solu- 
tion for the /;’s. It thus appears that the order of the approximation will be indicated 
by the number of constants of integration (other than F) which the solutions for the 
I;’s contain. Thus, as we shall presently show (§§ 3, 4, and 5, below), the solutions ob- 
tained by retaining only the first three, five, or seven terms in the expansion (5) involve 
respectively one, two, or three disposable constants of integrations. These, then, will pro- 
vide the first, second, or third approximations to the solution of equation (1).5 


* E.g., see Milne, of. cit., p. 121, eq. (125), or Unséld, op. cit., p. 97, eq. (29.13). 


5 One of the oversights of Gratton (loc. cit.) consists precisely in his failure to recognize this point. 
Thus he obtains his ‘‘second’”’ approximation by retaining terms up to and including J;. However, as 
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For further reference we may note here that the relations (17) for the cases / = 0, 1, 
and 2 explicitly reduce to 


27) (0) =Io(0) +31,(0) —$753(0) +7575(0) +...., (18) 

31,(0) = 37 (0) +31, (0) +$22(0) — (0) +...., (19) 
and 

37 (0) = 47, (0) + (0) + (0) — (0) Pins ’ (20) 


where on the right-hand sides we have retained all terms up to and including Js. Remem- 
bering that 7; = constant = 3F/4, we can write more conveniently in the forms 


219 (0) +¥e13(0) (0) +.... = %F, (21) 
21 (0) +§72(0) +.... (22) 

and 
(0) (0) +.... =F. (23) 


So far we have considered only the boundary conditions at r = 0. Turning next to the 
boundary conditions at tT = ©, it is evident that® 


Iy~3Fr and I,(r) for 1>0 boundedas r> (24) 


3. The first approximation.—The first approximation is obtained by retaining the 
terms J, /;, and J: in the expansion for J and considering only the first three equations 
(namely, for / = 0, 1, and 2) which result from equation (7). The equations for / = 0 
and / = 1 lead, as we have already seen, to the integrals (10) and (12). The equation for 
l= 2is 


Sah. 


Remembering that /; is a constant and that we have put J; = 0, it follows that 
| In=0. (26) 
Accordingly, our solution is 
3Fr+a; (27) 


and involves, as we see, the one disposable constant a. To determine this constant, we 
should use one or the other of the relations (17). Now the relation for / =-1 is given pref- 
erence over all the others, since (as was first pointed out by Milne) this insures that the 


we have indicated, the genuine second approximation is obtained only when the term J, is also included. 
A consequence of this is that Gratton’s solution (after correcting some further arithmetical errors in his 
evaluation of the constants of integration) fails to satisfy certain other necessary conditions of the 
problem (see n. 7, below). 


6 Cf. Hopf, op. cit. 
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emergent flux is the same as the constant net flux in the interior. Accordingly, we use 
equation (22) to determine a. We readily find that 


a=}3F. (28) 
Thus, 
= +9), (29) 
and the corresponding law of darkening is 
I (0,4) =3F (30) 


We therefore see that the first approximation on our present method reduces simply to 
the solution in the form first derived by Milne. 

4. The second approximation.—As we shall now show, the second approximation for J 
is obtained by retaining the terms up to and including /4 in the expansion (5). This re- 
quires the consideration of the first five equations which result from equation (7). Thus, 
in addition to the integrals (10) and (12), we have now to consider the following equa- 
tions: 


ta," 
3 
4 
= | 
From the first and the last of the foregoing equations we obtain the integral 
(32)? 
Equations (31) also lead to the following differential equation for 73: 
, 3dl,_, 9 
or 
23 


The solution of this equation appropriate for our present purposes is (cf. eq. [24]) 
I;= (35) 


where A is an arbitrary constant and 
a= V$§= 1.399. (36) 
From equations (31) and (32) we now obtain 


—$Aae; (37) 


—4tAae™. 


7 The existence of this integral reveals the inconsistency of an approximation in which /; is included 
but J, is not (see n. 5). A consequence of this inconsistency is that when the solution for /, including 
only terms up to /;3, is obtained, it is found that /(0, 0) # Jo(0), contrary to what should be expected to 
be identically true (since at the limb we necessarily see only the surface layers). 
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The solutions for the J’s are thus seen to involve only two disposable constants (namely, 
A and a) and therefore lead, as already stated, to the second approximation. 

To determine the constants A and a, we use equations (21) and (22) as our boundary 
conditions. Under our present assumptions (namely, that /5, J¢, etc., can all be ignored) 
these equations reduce to 


219 (0) +1473 (0) =4%F (38) 
and 
(0) +312 (0) =3F. (39) 


Eliminating J (0) between the foregoing two equations, we obtain 
g12(0) — (0) — (0) =F. (40) 


This equation directly determines A. We find 


(41) 
1+%a 


Accordingly, our solutions for 72, 73, and J; are 
I,=0.3102Fe-7; —0.5175Fe-* ; I,=0.4136Fe-*7. (42) 


Equation (38) or (39) now determines J, (0). We find 


I, (0) =0.4397F. (43) 
On the other hand, since (cf. eq. [12]) 
I, (0) =a—312(0) , (44) 
we have 
a=0.5638F. (45) 
Thus, our solution for Jo is 
= 3Fr+0.5638F —0.1241Fe-2, (46) 
which for later comparisons (§ 6) we write in the form 
Io(r) =3F [7+ ], (47) 
where now 
= 0.7517 — (48) 


To obtain the corresponding law of darkening, we start from the exact relation 
dr 


1(0, = To(r) (49) 
0 
and substitute for J) (7) the solution (46). We thus find that 
0.1241 


5. The third approximation.—To obtain the third approximation for J, we retain all 
the terms up to and including J¢ in the expansion for J and consider the first seven equa- 
tions which result from equation (7). Accordingly, in addition to the two equations 
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which lead to the integrals (10) and (12), we have to consider the following five further 


equations: 
7 dr 
4dI,, 5 
5 , 6 dIg_ 
6 dl; 


From the first, third, and fifth of the foregoing equations we obtain the integral 
(52) 


which is our present analogue of the integral (32) we had in the second approximation. 
Again, from equations (51) we obtain 


36 6 dl, ) 


143 dr? 13 dr 


5 
ONT de? "11 
Similarly 
99 dr? 9\dr 7 dr? 
dls 
{ (54) 
9 
Equations (53) and (54) can be re-written alternatively in the forms 
59 _ _ 20 
117 dr? 63 dr?’ 
(55) 
23 dls 20 


45 dr? dr?” 
To solve the foregoing equations we make the substitutions 


I;=Ae" and I[,=Be-*, (56) 


W 
S 
A 
\ 
a 
I 


1er 
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where A, B, and a are oer for the present piste We obtain 


23 20 sla 
Hence a? must satisfy the bogus 
Solving this equation, we find that a can have either of the two values 
a,=1.9825 or a,=1.1464. (59) 
Corresponding to these two roots for a, we must have 
—1.2706A, and B,=1.2368A2. (60) 
Accordingly, the solutions for J; and 7; appropriate for our present problem are 
I3= + (61) 
and 
—1.2706A,e—%17 (62) 


Using equations (51) and (52), we can now evaluate the remaining /’s. We find 
Iz = —0.8497 — 0.4913 Ages" , 
I,g=+0.0121 — 1.2995 , (63) 
Ig — 0.7733 . 


We thus see that our solution for the /’s involves three disposable constants (namely, Aj, 
A», and a) and leads, therefore, as we have already stated, to the third approximation 


for J. 

To determine next the constants A;, As, and a we use equations (21), (22), and (23) 
as our boundary conditions. Eliminating J» (0) between equations (21) and (22), we 
have 


(0) — (0) (0) (0) (0) = (64) 


Equations (23) and (64) directly determine A; and A». We find 
A, = —0.3931F and A,=0.2384F. (65) 

Substituting these values for A; and Az in equations (61), (62), and (63), we obtain 
(+0.3340e-17 — 0.11 F, ) 
I3;= F, 
I,= (—0.0048¢-%7 — 0.3098 F, > (66) 
I5= (+0.4995 F , 
Ig= (—0.5402e-%7 — 0.1844e-%7) F. 
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The value of J» (0) can now be determined from equation (64). We find 
I,(0) =0.4440F . (67) 
Equations (12), (66), and (67) now determine a. We find 
a=0.5307F. (68) 
Thus the solution for Jo is given by 
Io(r) = 3Fr +0.5307F — 0.1336F + 0.0469F ; (69) 
or, writing 7» (7) in the form (47), we now have 
=0.7077 — 0.1781 + (70) 


Finally, to determine the law of darkening we again start with the equation (49) but 
substitute for 7) our present solution (69). In this manner we find 


1(0, ») =F (0.5307 


6. A comparison of the second and the third approximations in their predictions regarding 
q(r) and the law of darkening: the exact law of darkening.—As is well known, Hopf 
and Bronstein have succeeded in deriving certain exact results concerning the solution 
of the equation of transfer (1) satisfying the boundary condition (13). Thus it has been 
shown (seg Hopf, of. cit.) that, if the solution for 7) be written in the form 


= $F [r+ q(7)], (72) 
then q (7) is a monotonic increasing function of 7 and that, moreover, 
(0) = = 0.57735 (73) 
and 


The integral on the right-hand side of equation (74) was evaluated numerically, and it 


was found that 
q(o) =0.710447. (75) 


These exact values for g (0) and q (~ ) should be compared with those given by equations 
(48) and (70) on our second and third approximations, respectively. We have 


q(0) = 0.5862 ; =0.7517 (second approximation) (76) 
(0) = 0.5920 ; =0.7077 (third approximation) . 


It is seen that, while the third approximation effects a marked improvement over the 
second in the agreement of g() with the exact value (75), it worsens somewhat the 
agreement of g (0) with its exact value (73). However, it appears that the third approxi- 
mation actually represents a substantial improvement over the second over the entire 
range of 7 except in the very immediate vicinity of the boundary r = 0. That this is so is 
indicated, for example, by the comparison made in Table 1 of the functions g (7) given by 
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our second and third approximations with that given by Eddington’s second approxima- 
tion.’ This conclusion is further strengthened when we compare the laws of darkening 


TABLE 1 


A COMPARISON OF THE FUNCTION g(r) DERIVED ON THE BASIS OF THE SECOND AND THE 
THIRD APPROXIMATION (EQS. [48] AND [70]) WITH THAT DERIVED ON 
EDDINGTON’S SECOND APPROXIMATION 


q(r) q(r) 

Approxi- Approxi- Second Approxi- Approxi- Second 

mation mation mation mation 

mation mation 

.5974 .6053 .620 2 .7108 . 7030 .694 
.6078 .6172 .635 .7208 .7070 .697 
.6175 .6280 .645 .7283 .7091 .699 
.6266 .6375 .653 .7340 .7102 
.6429 .6537 .665 .7416 .7106 .703 
.6503 .6605 .670 7441 .7104 .704 
.6695 .6768 .679 .7098 .705 

.6895 .6912 .686 | .7492 .7092 .706 

TABLE 2 


A COMPARISON OF THE LAWS OF DARKENING GIVEN BY THE SECOND AND THE 
THIRD APPROXIMATION (EQS. [50] AND [71]) WITH THAT 
GIVEN BY THE EXACT FORMULA (EQ. [77]) 


1(0, 1(0, w)/Z(0, 1) 
a Second Third Second Third 
Approxi- Approxi- Exact Approxi- Approxi- Exact 
mation mation mation mation 
0.4397 0.4440 0.43301 0.3484 0.3530 0.34390 
0.5299 0.5363 0.54011 0.4199 0.4264 0 42896 
0.2 0.6168 0.6232 0.62802 0.4887 0.4955 0.49878 
eee, 0.7014 0.7068 0.71123 0.5557 0.5620 0.56487 
0.7842 0.7884 0.79210 0.6214 0.6268 0.62909 
0.8657 0.8684 0.87156 0.6860 0.6904 0.69220 
0.9463 0.9475 0.95009 0.7498 0.7533 0.75457 
1.0261 1.0258 1.02796 0.8130 0.8156 0.81642 
Recht, Shane 1.1052 1.1035 1.10536 0.8757 0.8774 0.87788 
ae es 1.1838 1.1808 1.18238 0.9380 0.9388 0.93906 
TAs eee 1.2620 1.2578 1.25912 1.0000 1.0000 1.00000 


given by our two approximations with the exact law of darkening. This comparison is 
made in the following paragraph. 


8 Milne, op. cit., pp. 122-124, and Unsild, of. cit., p. 102. 
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It has been shown by Hopf (op. cit.) that for the case under discussion the problem of 
darkening admits of an exact solution. Thus, expressing the emergent intensity J (0, u) in 
the form 

V3 


Hopf has shown (op. cit., p. 105) that 


m2 


Although the foregoing solution has been known for over ten years, it does not appear 
that the integral defining ® (u) has been evaluated. Accordingly, we have evaluated the 
function numerically for the values = 0, 0.1, 0.2,...., 1.0 and tabulated the re- 
sulting values of J (0, uw) in Table 2. In Table 2 we have also tabulated the functions 
I (0, 4) as given by the equations (50) and (71), i.e., by our second and third approxima- 
tions, respectively. It is seen that, except for u very close to zero, equation (71) for J (0, u) 
provides an appreciably better approximation to the exact law of darkening than does 
equation (50). And, moreover, the agreement of either (50) or (71) with the exact law of 
darkening over the whole range of y is entirely satisfactory. 


In conclusion I wish to record my indebtedness to Miss Frances Herman, who carried 
out most of the numerical work connected with the preparation of Tables 1 and 2. 
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STELLAR ROTATION AND LARGE-SCALE CURRENTS 


WASLEY KROGDAHL 
Washington, D.C., and Yerkes Observatory 
Received October 6, 1943 


ABSTRACT 


The general problem of stellar rotation is considered. It is shown that in viscous stars rotation will 
give rise to meridional currents, which, in turn, cause a differential rotation. A formal solution for the 
differential rotational velocity is given in terms of the meridional velocities. The equations to be satis- 
fied by the meridional velocities are given, and an approximate method for their solution is suggested. 


Ever since the presentations of what may be termed the “classical” investigations on 
stellar rotation—i.e., those in which the fundamental theorems of uniform rotation were 
first given—interest in the subject of rotating stars has been directed more and more to 
the particular topic of differential rotation. This course has been, to a large extent, a 
consequence of von Zeipel’s enunciation of his well-known theorem, which would entail 
the impossibility of uniform stellar rotation under certain general conditions in stars in 
radiative equilibrium.! Moreover, a complete and satisfactory theoretical analysis of 
differential rotation in general, and of the sun’s differential rotation in particular, was 
wanting. It is natural, then, that numerous investigators have sought to solve both the 
general problem of differential rotation and various special cases; however, their success 
has, in some respects, been of a restricted nature because of either great generality or 
great specificity. 

The problem is that of finding a system of axially symmetric currents which either can 
or must exist in a star having a nonzero net angular momentum. The problem can be 
considered relative to any of four fundamental sets of circumstances: the stellar mate- 
rial may be either viscous or nonviscous (perfect), and at the same time there may or 
may not exist large-scale meridional currents in the star. However, Gunnar Randers 
has shown? that, if the fluid is perfect, stable meridional currents cannot exist; and S. 
Rosseland® has found that in nonviscous stars having no meridional currents the differ- 
ential rotational velocity is indeterminate to the extent of an arbitrary function. There- 
fore, attention may be confined to viscous stars. 

Motions in a viscous fluid will, of course, be subject to the decelerating forces due to 
viscosity; and the fluid velocity will, therefore, in general be a secularly varying vector- 
field. For the problem to be considered, however, it will be postulated that the velocities 
are stationary—i.e., independent of time—for it is only the equilibrium state of a rotating 
star which at present interests us. One therefore seeks only sustained motions. 

To characterize a motion in a viscous fluid as “sustained” is to imply the existence of 
a sustaining force. In a rotating star such a sustaining force must be ultimately traceable 
to the star’s over-all rotation; it must eventually be shown that any sustaining force is 


1 However, as Dr. S. Chandrasekhar has pointed out to the writer, von Zeipel’s theorem does not ap- 
pear in a theory in which only terms of the first order in w} (where wo denotes the angular velocity) are 
retained; it is only when terms of higher order in w% are included that von Zeipel’s theorem follows (cf. M. 
Schwarzschild, A p. J., 95, 452, 1942). This should perhaps be emphasized, since in the case of the sun, for 
example, w%/4rGp, ~ 10-5, and we cannot really be certain that the fundamental equations are valid in 
the higher order approximations, i.e., to at least one part in 10.!° 


2Ap. J., 94, 123, 1941. 
3 A strophysica Norvegica, 2, 173, 1936; ibid., p. 249, 1937. 
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generated by the rotation, for presumably sustained (steady-state) motion would not 
exist in a static star (in which no preferred directions could be defined). 

With these assumptions and restrictions, consider now the problem of general large- 
scale motions in a rotating star. 


I 


The distribution and behavior of the matter in a star must satisfy six general equa- 
tions; they are (in the convenient tensor notation): (a) the hydrodynamical equations‘ 


where p is the density, p the total pressure, B the gravitational potential, v,, the (co- 
variant) velocity vector,® A the coefficient of viscosity, g*’ the fundamental tensor, and 
where 


is the divergence of the velocity (g being the determinant | g,;|); (0) Poisson’s equation 
Ve ox Vi 4nGp ; (3) 
(c) the equation of continuity 
1 
(d) an equation of state 
p=P(p,T), (5) 
(where 7 is the temperature), usually taken to be 
k 
p + (6) 
(e) the radiant flux equation 
past 


where p” is the radiation pressure, x the opacity, c the velocity of light, and F, the radia- 
tive flux vector; and (/), the sab ile of conservation of energy 


1 


87 


(8) 


where 7 is the total subatomic and thermal energy per unit mass of stellar material, F; 
the total (radiative plus conductive) flux, ¢ the rate of generation of energy per unit mass, 
and g the heat generated per unit time and volume by viscosity. Equations (1), (3), 
(4), (5), (7), and (8) are the fundamental equations of the problem. 


‘ 4 — J. McConnell, Applications of the Absolute Differential Calculus, p. 282, London: Blackie & 
on, 1936 


bd og indices preceded by a comma denote covariant differentiation with respect to the corresponding 
variable. 
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For future convenience and consistency it will be found of service to express these 
equations in dimensionless variables according to the transformations. 


) 
P=P.y = (4nGp?r**) 
p=pS, 
B= | (9) 
T=T.o, 
= 


/ 


Using these variables and taking the co-ordinate system x! = r (the radius vector), 
x? = yw = cos ¢g (where ¢ is the polar angle), x? = @ (the azimuthal angle), with the cor- 
(contravariant) velocities = (44Gp, r*)1/2 = (4aGp.r*)'/? B*, and 

= (4nGp.r*)'/? w*, the hydrodynamical equations (1) for a stationary (i.€., 0Um/dt = 0), 
axially star are, in full, 


1. a ap*] 17.,0%* da* 


dé 
an 


/ 


The condition that all velocities vanish with the rotational velocity plus the obvious 
properties | that a* and 6* must be even functions of w (the over-all angular velocity), 
whereas w* is to be odd in wo, make it possible to choose velocities of the form 


= wo (1 + + wjw’+....), 
(13) 
= (B+ 


0- 
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)) 
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where w, w’ peeeeg Gy a’,....,8,8',...., are functions of and u. Similarly, one may 
for convenience represent the ‘dimensionless quantities X,...., Corresponding to 
the pressure, density, potential, etc., by expansions in w?, namely, — 


n(&, = mo (E) (&, +...., 
x (E, = x0(E) (E, w) +...., 


(14) 


etc. By introducing expressions (13) and (14) into equations (10)-(12) and equating 
terms of like order in wo, one obtains the equations to be satisfied by the rotational per- 
turbations of the various orders. This procedure is simply an expedient by which w, 
may be eliminated from the equations to be dealt with. Specifically, the equations of 
the first-order perturbation effects are 


1 dno 1 1 07a 0a 


15 

d 


From the form in which equations (15) and (16) have been written it is clear that they 
are gradient equations. They must therefore satisfy the condition (applied to the right- 
hand sides) 


which, when applied to the left-hand sides, becomes, after collecting and rearranging 
terms, 


+ (1-2) ae (apt) [#8 + 9) 
where 
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At this point the equations may be made somewhat more manageable by a suitable 
choice for the forms of the solutions. Thus, without loss of generality,’ one can take 


k 

k 

x = x , (21) 
k 


= >) ax (£) Pe (x) , 
k 


u) = (E) Se (u) , 
k 


where P;.(u) is the Legendre polynomial of kth order and where 


dP, k (k +1) 


(22) 
du 


The substitution of (21) into (19) reduces the latter to a set of equations in & alone, 
namely, 


d 1 duo ) 
d 


x +1) (£81) [ (Ea) (23) 


—k(k+1) (88) |= | 


where 


=F | So fo nae (24) 


It is to be noted that equation (23) is an equation of condition. It merely imposes re- 
strictions on a, and §; (or on m,x% and ¢1,.) which, when satisfied, make it possible to write 
the first two hydrodynamical equations, (15) and (16), as gradient equations. With solu- 


6 It is well known that any suitably behaved function of u may be expressed uniquely in the Lapresite 
polynomials; hence it is obvious that there are no unnecessary restrictions placed on 9), ¢@), x, or a 
by assuming solutions of the form (21). In addition, 8 could likewise be expressed in terms of the P;, but 
it is a simple matter to show that the only solution in this form which satisfies the equation of continuity 
would allow of a grouping of terms which would bring it to the form given in (21). 
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tions of the form assumed and satisfying (23), equations (15) and (16) reduce simply to 
the set of equations 


= —k (k+1) 


dé 2d 
El dé g (25) 
k 
where 


Thus far, the only equations of the problem to be considered have been the hydro- 
dynamical equations, which have been brought to the form given in equations (17), (23), 
and (25). In addition, the various functions must satisfy the perturbation equations de- 
rived from (3), (4), (5), (7), and (8). From Poisson’s equation ([3]) one readily obtains 
the first-order equation 


dx; 
aig (® mrt) xy = (27) - 
while from the equation of continuity ({4]) one finds (for stationary motion) 
1 da 1 dfo 
— —k(k+1 =0. 28 
dé a): (k +1) By (28) 


Combining equations (25) and (27) and setting 


1/1 d 2 
one has (from ng [24], [25], [27], and [29]) 
d®,, N1y k 
where 


(31) 


My k 
dé 
Equations (23), (24), (28), and (30) are the principal equations pertaining to the veloc- 


ities in the meridional plane. A discussion of their solution will be considered later. 
First consider the third hydrodynamical equation. 


II 


Consider equation (17), to be satisfied by the differential rotational velocity. It can 
be written in the form 


he 


é 
A 
it 
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Let w have the form? 


=} we Me (u) (33) 
k 
where 
Mi (u) (34) 
Then, since (from [34]) 
+k (k+1)M.=0, (35) 


equation (32) reduces to 
1 


k 


But 
P;= 2741 “ 1 M;-1] 
and 
Ms 


“7 Mil, 


so that, for each k, one has from (36) 


1 Qk+1 | 


2k —1 2k+3 


= 


A formal solution of (37) can easily be obtained by the method of variation of parameters; 
it is 


1 fo 1 
= get (ws x—* g, (x) ax) + (w+ 


xf gy (x) dx), 


(38) 


7 The function M; is found as a solution of equation (35), which must necessarily be satisfied if w is to 
have a solution separable in — and yw. The set of functions M; is complete, and all M; are bounded for 
— 1 <¢ w <¢ 1, the only other properties which they are required to possess. 
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where w, and w; are constants of integration. However, if this solution is to remain 
finite at £ = 0, one must have 


1 
k+1 


so that, in final form, 


—k 
= (m S g, (x) (x) dx. (39) 


The constant w, is fixed by the boundary condition that the stress due to w vanish at 
the surface. This requires that® 
d Wk 
—=0 = 
Here two points are worthy of note.” The first is that, quite generally, if there are no 
meridional currents (i.e., if g, = 0), there can be no differential rotation, for then 


which will satisfy the condition (40) only with w, = 0. This makes it clear that a differ- 
ential rotation is a consequence of the existence of meridional velocities and does not call 
the latter into being to preserve stability. Second, consider the behavior of w; as vp > 0. 
The function g, will become infinite if a, and 6; remain finite, in which case conditions of 
stability will be violated, as Randers has shown. On the other hand, if a, and 6; are 
made to vanish, the function g; may be made to approach any function whatever. This is 
precisely Rosseland’s result. 

It should also be pointed out that, as was stated before, the differential rotation must 
be motivated and sustained by some generating force. This force is represented by the 
function g, and, as is to be expected, the vanishing of g entails the vanishing of w. 

It has been shown that the vanishing of a and B implies the vanishing of w. It is a 
simple matter to prove the converse theorem. Thus, from (32) it follows that for w = 0 
the condition 

Eup 


must be satisfied. In general, such a condition is incompatible with the equation of con- 
tinuity and the hydrodynamical equations, except for a = 0 = 8, for it imposes a special 
condition on an already determinate function a. Hence, in general, the condition w = 0 
(i.e., no differential rotation) implies the result a = 0 = 8, proving the theorem. 

Now it follows further from (23) and (24) that with no meridional velocities the func- 
tions /; all vanish identically, thereby satisfying the necessary and sufficient condition 
that 7 be a function of ¢*—that there exist a pressure-density relation. Hence it has been 
shown that the existence of a pressure-density relation is both a necessary and a sufficient 
condition for the existence of uniform rotation, and conversely. 

It is an interesting corollary of this theorem that the structure of a uniformly rotating 
star depends only on the pressure-density relation of the corresponding unperturbed con- 
figuration. For, as has been shown, if the star rotates uniformly, there must exist a 
pressure-density relation. This relation must be the same as that which governs the un- 
perturbed configuration, for the pressure and density at every point are the sums of their 
respective unperturbed values and first-order perturbations. But, since a pressure- 


8 See eq. (47). §W. S. Krogdahl, Ap. J., 96, 128, 1942. 
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density relation does exist, the first-order perturbations can be found from a knowledge 
of the pressure- and density-gradients of the unperturbed configuration.!° The per- 
turbed pressure and density thus found satisfy the original pressure-density relation. 

As a special case, one may consider two static configurations—one a polytrope of in- 
dex n = 3.25 with K an absolute constant, the other a model in radiative equilibrium in 
which the energy generation is uniform and the opacity is given by Kramers’ law; the 
latter model is also known to be built on a polytropic model with n = 3.25. The pressure- 
density relations, expressed parametrically as p = p(r), p = p(r), are identical in the 
two cases. The corollary states that the corresponding uniformly rotating configurations 
will likewise have identical pressure-density relations. 


Ill 


The problem is now reduced, in principle, to that of finding solutions to equations 
(23), (24), (28), and (30) and imposing on them the appropriate initial and boundary 
conditions. 

Consider the boundary conditions which are to be satisfied by the velocities. First is 
the obvious condition that the velocity normal to the surface must vanish. Now the sur- 
face is defined by the condition 


f= (41) 


Hence the normal to the surface will have the components 
dg ly k 


(42) 


Therefore, the condition on the normal velocity, namely, 
= 0 (r=1»), 
gives the first-order condition v! = 0 (r = 79) or 
a; = 0 (= &). (43) 


Three additional conditions on the velocities are that they be such that the viscous 
stresses shall vanish at the surface. This means that 


= 0 , (44) 


where Tmn = 3(0m,n + Un,m) is the velocity stress tensor. Using (42), the three first-order 
conditions are seen to be 


Ti =0, Tn =0, = 0 (r=1); (45) 
and, since 
dvi 1 r2 v? 1 dv? 


10 Tbid. 
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together with (43) they reduce to 


At & = 0, on the other hand, the radial velocity and its first derivative must vanish; 
i.e., 
a,=0, (€=0). (48) 


The second of conditions (47) is seen to be a condition on ®;, but this may be trans- 
lated into a condition on a; by use of the equation of continuity. For, from the equation 
of continuity, one has 


day 1 dfo 


and from the boundary conditions it is possible to express a; and 6, near the boundary 
ay = af) Eo) 2+ (E— Eg)? 

= Bi + Bi?) (E — Eo) *+..... 
Also, the function (2&!+ {7 'dfo/dé) can be put in the form 


1 


Introducing these series forms into the equation of continuity, 
(k+1) + ] CE — (3 FAK) +1 af?) 
—k (k +1) BM] (&— &)*+....=0. 


or 
Bo =0 


and, since in general 2 + I‘) # 0, 


(2) 


Therefore, the second of conditions (47) may be replaced by the condition 
da 
(£= &); (49) 
and from equations (47), (48), and (49) it is clear that any acceptable solution a, must 
be expressible as 
a, = A, (€), (50) 


where A;(£) is a bounded continuous function of £ in the interval 0 < & < &. 

Several further properties of the velocities may be known without ‘actually obtaining 
their solutions. First, it is clear that for all nonzero a, and 6;, k must be even to obtain 
solutions for a and 8 respectively even and odd in y, while for all nonzero w,, j must be 
odd to obtain a solution w which is an even function of u. This property is consistent with 
equation (37). Further, a9 and 8) must vanish, since they correspond simply to a radial 
expansion, which cannot exist in a star in equilibrium. 
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Second, consider the boundary condition dw;,/d— = 0 ( = &) on w;. From (39) one 


has 
—k—2 


1 g 


whence 


(k +2) 


At & = &, 


~k-3 
(k ged = 0: 


dw; 


dé 4 = (k—1) + 


or 


+2) 
2k+1 


x*gidx=0. 
0 


81 = (3 EB. — ay) ; (54) 


and, since in general the second of equations (53) holds only if g: = 0, one must in gen- 
eral have ag = 0 = By." 

The only other initial conditions to be specified are those which must be satisfied by 
the temperature and density perturbations. If one chooses to consider configurations for 


Eo 
(53) 


But, by (37), 


which” 
fy o1,4=0 (€=0), (55) 
then, after requiring also that 
_ 


the problem is completely determinate. 

One now has to solve the equations (23), (24), and (30), subject to ‘one (28), (7), (8), 
and the initial and boundary conditions. Clearly, any attempt to find explicit and rigor- 
ous solutions of these equations would soon involve one in an impossible maze of equa- 
tions and conditions. At best, therefore, any direct procedure must, for practical reasons, 
be at some point short-circuited. But, first, it may be asked whether any nontrivial solu- 
tions (i.e., not identically zero) for a;, Bx, 71,x, etc., exist. .To answer this question, note 


1! The possibility 
a2 = 


is excluded, since, by combining it with the equation of continuity, one finds that 
d 2 
di log 0, 


2 The particular initial values of {1, , and o1, , are immaterial and are taken as in (55) only for conven- 
ience. Therefore, a specification of these initial values does not amount to fixing two disposable constants, 
but rather to the adoption of a particular convention. Equations (56) do, however, correspond to the 
fixing of two constants of integration, namely, those which arise from the radiant flux equation and the 
equation of conservation of energy—equations (7) and (8), respectively. 


a= 


= 
| 
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that by combining equations (23), (28), (29), and (30) one can obtain an equation in 


a, of the form 
bs « (424). (57) 


Since this is an ordinary linear differential equation, a nontrivial solution can, in general, 
be expected to exist, provided the number of conditions imposed upon it does not exceed 
the number of available adjustable parameters; such a solution will be unique if the 
number of conditions equals the number of adjustable parameters. From this result there 
follows another, of considerable importance for general purposes. Thus, it is,clear that 
to require that 


d d 1 dv ) 


—k(k+1) =0 | 


is equivalent to imposing on a, an infinite number of conditions in addition to the two 
initial conditions and three boundary conditions. Hence, this relation can be expected 
to be true only if a, = 0 (= 6, by the equation of continuity). That is, f, = 0 in general 
implies a, = 0, and conversely. 

The significance and interpretation of this result can be seen in the following argu- 
ment. Ideally, a star in static gravitational equilibrium will be a radially symmetric con- 
figuration. Mathematically this means that any function describing some property of 
the star will be a function of a single variable, such as the distance from the center. 
Thus, if two or more stellar functions are specified as functions of the radius vector, this 
specification may be regarded as a statement in parametric form of the dependence of 
these functions on each other. This fact may be looked upon as a sort of mathematical 
degeneracy implied by the radial symmetry. If the star is made to rotate, however, one 
degree of the degeneracy is removed, for now the specification of any function through- 
out the star must be made in terms of two independent variables (such as é and uy, for ex- 
ample) and any two stellar functions will, in general, be independent of each other. In 
particular, the pressure and density will im general be independent functions. (As will 
be seen, an important exception to this general rule is that where the star rotates with 
a constant angular velocity [see next paragraph]). Thus, if the pressure and density are 
independent stellar functions, then in any region within a rotating star any other func- 
tion is expressible as a function of both p and p. Specifically, if the gravitational potential 
is ¥ = B(p, p), the isobaric and equipotential surfaces will not coincide. Therefore, the 
pressure- and gravitational-gradients cannot be parallel; their difference is an accelera- 
tion (vector) which must be compensated in some way if the star’s structure is to be sta- 
tionary. In a viscous star the compensation will be afforded by the viscous and inertial 
forces acting in opposition to the accelerating forces. This is merely a verbal statement 
of the hydrodynamical equations and, when put thus, seems obvious. 

This argument does not apply if, contrary to assumption, the potential is either inde- 
pendent of p or if p is a function of p. But, if ¥ is to satisfy Poisson’s equation, it cannot 
be independent of p. On the other hand, the necessary and sufficient condition that p be 
a function of p is that!® 


13 See“above, n. 9. 
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which has already been seen to imply, in general, the vanishing of a and 8. Clearly, then, 
the function f is the generating function which gives rise to the meridional currents. It 
is a result of the over-all rotation in the sense that rotation removes the degeneracy 
which implies the vanishing of f. 

It is clear, now, how the differential rotation is brought about. The over-all rotation 
creates a force, in the manner just described, which generates large-scale currents in the 
meridional plane. These, in turn, provide a rotational acceleration which gives rise to a 
differential rotation; the meridional currents are second-order effects, the differential 
rotation a third-order effect. The cause-and-effect relation of the two phenomena is 
sometimes reversed, and for this reason it is well to specifically point it out. 


IV 


It is clear that to obtain rigorous solutions for all the fundamental equations simul- 
taneously is an almost impossible task. A convenient method of approximation would be 
desirable. One method is suggested by equation (50). This method makes use of a well- 
known theorem that any bounded continuous function may be approximated in a par- 
ticular interval to any arbitrary degree of exactness by a suitable polynomial. Therefore, 
let A; in (50) be given as a polynomial in — with landletesmined coefficients. To deter- 
mine the coefficients, first determine the solutions aD, oi, ...., ete. , corresponding to 
the particular polynomial A{™ = &. In so doing, use will be made of equations (23), 
(24), (28), (29), (30), (31), and (5), but not of (7) and (8). Then the coefficients may be 
so chosen as to satisfy equations (7) and (8) as nearly as possible.'* Even this procedure, 
while representing a considerable simplification, is at best laborious, and it will not be 
undertaken here. However, it is quite simple and may prove instructive to derive the 
meridional streamlines for the various “harmonic components” of the meridional veloc- 
ities. 

The streamlines for the kth harmonic component of the meridional velocities are de- 
fined as curves along which 


dé du 
= 59 
or 
d dP, 
ak 2 dP, 
(1 
But, by the equation of continuity, 
Bi _ 1 
k(k+1)— Soar) log (£?Soax) , 
so that, from (60) 
P, 
dlog (1 — 42) “tt 0 
and the streamlines are the curves 
2 2 
[ (l—p =¢,=constant (k= 4, 6,8,....). (61) 


Figure 1 shows the shape of the streamlines for the case Ay = &, {) = (43 26)* 25. where 
93.25 is the Emden function of index 3.25. Note that the streamline c, = 0 is formed by 


14 That is, either a best fit, in the sense of least squares, or exactly at a prescribed number of points 
equal in number to the number of adjustable coefficients in the polynomial A x. 
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the outer boundary ({) = 0), the axis (u = 1), the equatorial plane (u = 0), the cones 
over which dP;/du = 0, and the spheres over which a, = 0. 

Finally, it may be found of interest to fit a solution of the form (33) for to the rota- 
tion observed in the sun. No pretense is made of having thus obtained a rigorous solu- 
tion; and the results are in no sense intended as data against, or in support of, the theory 


Fic. 1.—Meridional streamlines of the first non-vanishing harmonic component for the special case 
Ag = & £0 = 


TABLE 1 
Latitude red." 


given in Section II, for any smooth function of u could be well fitted by a polynomial in 
the dP;/du. By using the observed solar rotational velocity at latitudes 60°, 30°, and 0°, 
the angular velocities at other latitudes may be calculated from a three-term series. 
Table 1 shows the observed" and “predicted” values for four latitudes. 


It is a pleasure to acknowledge my indebtedness to Dr. S. Chandrasekhar for his un- 
failing encouragement and helpful advice and to Mr. Robert M. Snow for numerous 
fruitful discussions. 


6 Observed data are taken from the Handb. d. Ap., 4, 167, Table VII. 
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NOTES ON STELLAR SPECTRA* 


Otto STRUVE 
McDonald and Yerkes Observatories 
Received September 15, 1943 


ABSTRACT 


The P Cygni star HD 169226 = MWC 293 shows radial velocities ranging from — 143 km/sec for 
H to — 50 km/sec for Mg u. The emission lines of H give + 89 km/sec, and the emission lines of He 1 
give + 141 km/sec. 

The star HD 217050 = MWC 394 has a peculiarly poor spectrum on the violet side of \ 3647, with 
only Nim very strong. The other lines, Feu, Ti 11, Cr 11, etc., must be suppressed by continuous H 


absorption within the shell. 
The color of HD 160529 = MWC 266 has been found to be similar to that of a late F or early K star. 


The spectral type is A4sea. 

Pleione has shown, in July, 1943, an increase in intensity of Sc 1 and Fe 1, which now definitely es- 
tablishes the evolution of its shell toward lower ionization. 

AG Pegasi showed, in the summer of 1943, P Cygni absorption lines of increased strength, while Z 
Andromedae had stronger emission lines of Fe 1 than in the preceding fall and winter. 


HD 169226 


This star is MWC 293! and is described by Merrill as showing an H@ line of the P 
Cygni type on a spectrogram taken at Mount Wilson on July 20, 1924. Two spectro- 
grams were taken at the McDonald Observatory, on July 7 and July 8, 1943, with the 


TABLE 1 
RADIAL VELOCITIES OF HD 169226 
Atom No. of Lines Velocity Atom No. of Lines Velocity 
H absorption...... 7 —142.8 km/sec || Cau absorption... 1 — 10.8 km/sec 
H emission........ 3 + 88.6 N wu absorption... . 1 => 78.9 
Het absorption... . 7 — 84.0 O 1 absorption... . 7 — 61.0 
Het emission...... 1 +141.4 Sim absorption... 3 — 95.9 
Mg 1 absorption... 1 — 49.8 C 11 absorption... . 1 — 80.6 


Cassegrain quartz spectrograph, giving a dispersion of 40 A/mm at \ 3933. The star 
is very red and therefore represents a difficult object, even for a powerful instrument. 
The P Cygni character is very pronounced. In fact, among the P Cygni type stars which 
I have observed, it resembles P Cygni more than any of the other members of the group. 
But the emission lines are conspicuous only for H and weak for Het. The other atoms 
have only absorption lines. Ca 11 is probably interstellar and is very strong. There is a 
conspicuous broad absorption feature at \ 4430, agreeing with Beal’s unidentified inter- 
stellar band in this position. The mean radial velocities from the two plates, in Table 1, 
show a large range, but there is no obvious relation with ionization potential, as in the 
case of P Cygni. The large positive velocities obtained from the emission lines are of 
interest in connection with the analogous result recently obtained by Beals? for MWC 
374. 


*Contributions from the McDonald Observatory, University of Texas, No. 81. 
tq = 18518™6; 6 = — 12°15’ (1900), mag. 9.1. 2 J.R.A.S., Canada, 37, 245, 1943. 
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HD 217050 


The spectrum of this shell in the ordinary photographic region has been recently de- 
scribed by R. B. Baldwin.* The star is MWC 394 or HR 8731.4 Several early Yerkes 
spectrograms are valuable in establishing the variations of this spectrum. These plates 
have already been used by E. B. Frost in establishing the fact that the emission lines of 
H are variable in this star.* The shell features were, however, not recognized as such 
at that time and remained unnoticed. On July 28, August 4, and October 9, 1913, the 
H lines showed strong sharp cores, which are typical of shell absorption. On August 5 
and 30, 1918, these shell absorptions were completely missing: they had disappeared, 
together with the emission lines. The next Yerkes plate is one taken June 4, 1929. It 
shows the sharp cores of H, together with double emission lines and probably a few 
faint, fairly narrow lines of Fe 11. A longer series of plates was started on September 15, 
1941, and shows the shell spectrum very much as described by Baldwin. I have the im- 
pression that there are fairly rapid changes in the intensities of the a Cygni lines, but 
they never disappear as in AX Monocerotis or HD 218393 and only seem to be weaker 
at certain times. I have not been able to find a period, but the suspected changes take 
place within a few days. The radial velocities of the Yerkes plates are given in Table 
2.6 They are based upon the lines He, H6, Hy, HB, Catt K and half a dozenof the 
stronger Fe 1 lines. 

The ultraviolet spectrum of HD 217050 is of considerable interest. I obtained a 
strongly exposed spectrogram with the McDonald Cassegrain quartz spectrograph on 
Eastman Process emulsion and have measured on it the lines listed in Table 3. The 
date of the plate is June 22, 1943, 10"28™ U.T. Considering the fact that the ordinary 
photographic region is rich in lines which, though fainter than in Pleione or 48 Librae, 
crowd our plate on the red side of the region which I have measured, the poverty of 
lines on the violet side of the Balmer limit is an interesting example of continuous H 
absorption within the shell—a phenomenon to which I have already called attention 
in the case of 48 Librae.’ There is another interesting similarity between HD 217050 and 
48 Librae: both stars show the ultraviolet lines of Ni 1 with great strength and appar- 
ently with little effect from the continuous H absorption. At least, there seems to be no 
noticeable weakening of these lines as compared to Ni 11 4067. It is probable, though 
perhaps not certain, that Ni 11 occupies a relatively high level within the shell. 


HD 160529 


This star is MWC 266° and was last discussed by Swings and Struve® and by Merrill 
and Burwell.'? The spectrum shows very strong a Cygni lines, whose radial velocities 
suggest expansion in varying degrees. Merrill remarks that the interstellar line \ 6284 
is very strong. Swings and I pointed out earlier that the H lines are weak but that the 
star is so red that it was difficult to observe it in the ultraviolet region. Repeated efforts 
to get the ultraviolet region with our normal dispersion of 40 A/mm at A 3933 having 
failed, I obtained in the summer of 1943 several plates with the dispersion of 100 A/mm 
at \ 3933. These plates were obtained in part to get the ultraviolet spectrum and in part 
to determine the reddening. For the latter purpose the slit was opened quite wide and 


3 Ap. J., 97, 388, 1943. 

4 Also BD+47°3985; a = 22552™7; 5 = +48°9’ (1900), mag. 5.2. 

5 Ap. J., 49, 61, 1919. 

6 The radial velocities after August 24, 1943, were added in the proof. 
7Ap. J., 98, 98, 1943. 

8q@ = 17535™3, § = —33°27’, mag. 6.7, Sp. A4sea. 

9 Ap. J., 91, 592, 1940. 10 Ap. J., 98, 356, 1943. 
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TABLE 2 
RADIAL VELOCITIES OF HD 217050 
Date U.T. Quality Date U.T. Quality 
1929 June 4...) 8:48 |—12.1 km/sec. p 1943 July 22...) 7:08 |—25.9 km/sec. g 
1941 Sept.15...) 7:19 | 20.1 g July 28 4:10; 18.0 p 
Sept. 26...) 7:04 | 13.0 g July 28 5:20 | 10.8 p 
Sept. 27.../ 1:45 | 15.4 g Aug. 10 3:48 | 21.1 g 
Oct. 126 g Aug. 10 5:03 | 29.9 g 
Oct. 16...| 1:16 | 21.5 g Aug. 17 5:50 | 20.8 g 
Nov. 18...| 2:12 | 18.6 g Aug. 18 4:42 | 23.4 g 
Nov. 19...| 2:02 | 27.2 g Aug. 24 3:22 | 29.9 g 
Dec. 20...| 0:07 | 22.4 g Aug. 24 4:52 | 20.8 g 
1492 July 3...) 6:18 | 24.6 g Sept. 21 1:25 | 23.8 g 
Sept. 1.../ 5:06 19.5 p Sept. 23 p 
Sept. 16...| 2:57 | 13.1 g Sept. 24 is] iA g 
Sept g Sept. 28 2:34 | 23.4 g 
Sept. 25...| 6:02 | 19.0 g Sept. 29 1:38) 23:2 p 
Sept. 27...| 5:09 | 21.3 g Oct. 6 1:22 | 19.7 g 
Sept. 28...| 2:48 | 19.9 g Oct. 7 1:28; 27:3 g 
Oct. 6...| 2:07 | 30.2 g Oct. 8 1:06 | 27.4 g 
Oct. 7...11:31 g Oct. 15 1:39 | 19.2 g 
Oct. | g Oct. 20 0:50} 11.1 g 
Oct. 19...) 1:48 | 20.7 g Oct. 21 1:50 | 19.0 p 
Oct. 28...| 0:24 | 24.9 g Oct. 27 0:41 | 25.3 g 
Nov. 13...| 1:00 | 19.7 g Oct. 29 0:12 | 29.7 g 
Nov. 27...| 0:00 | 20.0 g Nov. 28 1:12] 394 g 
Dec. 4...| 0:42 | 19.4 g Dec. 5 3:05 | 24.1 g 
Dec. 5-..| 1:03 | 17.6 g Dec. 24 2:49 | 16.0 g 
Dec. 12...| 0:52 | 23.9 g Dec. 25...|23:47 | 23.1 g 
1943 Apr. 8...|10:10 | 17.8 g 1944 Jan. 3 2:10} 27.4 p 
May 28...| 8:08 | 26.7 g Jan. 10 1:00 | 26.3 g 
June 11...) 7:26 | 21.2 g Jan. 17...|23:56 |—28.0 p 
July 22...| 6:07 |—25.9 g 
TABLE 3 
ULTRAVIOLET ABSORPTION LINES OF HD 217050 
Int. Identification Int. Identification 
0 Cr 11 39 .80(50) 1 Fe 21.27(6) 
On Ti 11 72 .80(400) Inn | He134.24(2); 
1 Nit 73 .98(4) (Cr 31.72(40)) 
In (Fe 11 87 .84(2)) 1 Ti 59 .76(150) 
1 (Cr 11 03 .32(100)) 1 H31 61.22 
1 Ni 07 .32(8) 1 H30 62.26 
5432.98... In Cr 11 33 .30(75) ¢ 2 H29 63.40 
2 Nit 54.16(5) 3 H28 64.68 
3 65 .62(1) 3 H27 66.10 
1 Fe 11 68 .68(8) 4 H26 67 .68 
2 Niu 71.35(2) H25 69.47 
5 Nim 13.94(8) 6 H24 71.48 
Onn | Tit 20.26(18) 6 H23 73.76 
3 Nit 76.76(3) 7 H22 76.36 
9B. 1 Fe 81.19(1000) 7 H21 79.36 
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was oriented north and south, so that all the light of HD 160529 and the neighboring 


comparison stars was admitted, in spite of the large amount of atmospheric dispersion _ 


at 6 = — 33°. The comparison stars were taken near the meridian and immediately 
after HD 160529; and, since they are located within a fraction of 1° from this star, differ- 
ential extinction can be neglected. The exposures were 5 minutes each; the night was 
very clear and the seeing steady. The results are summarized in Table 4. We estimate 
that the energy distribution of HD 160529 resembles that of a late F or early K star. A 
long exposure with narrow slit, taken on July 7, 1943, shows that the H lines are very 
narrow and weak, though they can be seen to about H24. They are thus much weaker 
than in normal supergiants, like a Cygni, or in shells, like 48 Librae, but they are 
stronger than in v Sagittarii, where Greenstein could measure #24 on plates of more than 


TABLE 4 


COMPARISON OF ENERGY DISTRIBUTION IN HD 160529 AND NEIGHBORING 
STARS ON JULY 7, 1943 


Mac. 
STAR ‘ines HD REMARKS 
(U.T.), SPECTRUM 
Ptm. Ptg 
HD 160529....| 5:47 (Oe5) 6.7 7.6 
160748....| 5:58 Ma 6.5 8.0 Redder than HD 160529 
160461....} 6:50 AO 7.8 ie Bluer than HD 160529 
160410....| 6:59 FO 8.4 8.2 Bluer than HD 160529 


twice the dispersion and much more contrasty emulsion.!! The spectrogram of HD 
160529 shows no noticeable Balmer discontinuity—which is in accord with the intensities 
of the H lines. The metallic lines in the ultraviolet region appear to be normal in in- 
tensity. 

PLEIONE 


Spectrograms taken by Messrs. C. A. Bauer and G. Miinch on Eastman Process 
emulsion show that the evolution of this remarkable shell spectrum since the winter of 
1942-43 has proceeded in the direction of increased intensity for Sc 11, Fe 1, and probably 
Tit. This change shows almost conclusively that the general evolution” since 1940 
has been in the direction of lowered ionization. There have been strange departures 
from normal ionization conditions, but the general trend has been from Ni 0, I.P. 18.2 
volts, in 1940, to Sc u, I.P. 12.8 volts, in 1943. This result will help in establishing the 
ionization sequence of other shells. The recent plates were taken on July 19 and July 
29, 1943. 

HD 193182!8 


In 1923 Merrill announced that HD 193182, or MWC 632, has a bright Ha line and 
that H@ is an absorption line with bright borders. The star was classified by him as F, 
although the Henry Draper classification was given as AO. In 1939 Morgan" pointed 
out that this star has all the characteristics of high luminosity usually associated with the 
designation c, and hence he assigned it to class cAQ, with an absolute magnitude of 
— 5.0. The star is located in the vicinity of y Cygni, where several other supergiants of 
early type were noted. But, while several of these supergiants are very greatly reddened 


1 4p. J., 91, 438, 1940. 
2 4p. J., 97, 428, 1943. 
13.q = 20138, 5 = +39°16’, mag. 7.0. 14 Ap. J., 90, 632, 1939. 
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by space-absorption, Morgan stated that “HD 193182 is probably not reddened more 
than one- or two-tenths of a magnitude.” Quite recently Merrill’® suggested that this 
spectrum may be a “‘false c-star,”’ like 48 Librae or Pleione. 

In order to test whether the spectrum is that of a shell, several plates were obtained 
in June and July, 1943, at the McDonald Observatory with the Cassegrain quartz 
spectrograph giving a dispersion of 40 A/mm at \ 3933. The spectrum is exceedingly 
rich in lines and leaves no question as to its origin: the hydrogen lines consist of very 
narrow and deep cores, which can be seen to H32. These cores are symmetrically super- 
posed over wide wings suggestive of Stark effect; the lines of Mg 11 4481 and of Si 11 are 
very weak and show that dilution of radiation is conspicuous. There are probably present 
broad and shallow lines of He 1 4026 and 4472, so that the underlying stellar spectrum 
may be that of a fairly rapidly rotating B star. The shell spectrum is of particular 
interest, because it is rich in lines of relatively low ionization potential; for example, 
Fet is strong. This spectrum is, therefore, the latest shell spectrum thus far observed, 
except that of 14 Comae Berenices."* In this latter star there are large differences in 
rotational broadening observed in lines of different atoms. No such effect is seen in HD 
193182 with the dispersion used in this work. On the McDonald spectrograms Hf is a 
sharp narrow absorption line, probably without bright edges. Ha is not shown. 


AG PEGASI 


A spectrogram of this star’? taken on July 5, 1943, is interesting because of the great 
strength with which the P Cygni absorption lines are now visible in the higher members 
of the H and He1 (28P — n®D) series. Plate VII shows the progressive change in this 
spectrum since 1939. 

Z ANDROMEDAE 


The principal change shown by a spectrogram taken on June 28, 1943, from plates 
taken in January, 1943, is a strengthening of the bright lines of Fe1. The decline in 
excitation recorded some months ago'® has thus continued. The spectrogram is re- 
produced in Plate VIII. The line Fe 1 4045 is now quite easily visible, although it had not 
been measured in 1940, 1941, and 1942.'° 


Iam indebted to Miss J. Ringstad and Miss G. Peterson for assistance in the 
measurements of spectrograms of HD 217050. 


18 Ap. J., 98, 364, 1943. 

16 Ap. J., 94, 305, 1941. 

17 BD + 11°4673 = HD 207757; see also Merrill’s detailed investigation, Ap. J., 95, 386, 1942. 
18 Swings and Struve, Ap. J., 98, 93, 1943. 


19 Swings and Struve, Ap. J., 98, 200, 1943. Note added in proof: On a spectrogram of November 14, 
1943, the excitation is again much greater. The lines of Fe 1 are much weaker, and those of [Fe vi] and 
{Ne v] much stronger than in June, 1943. The line [Fe vit] 3586 is very diffuse. 
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ABSTRACT 


HD 163181.—A new orbit has been derived for this eclipsing variable from 55 McDonald Observatory 
plates. The elements agree with those determined in 1928 at Mount Wilson, except that e is 0.10 instead 
of 0.065. The velocity-curve shows a conspicuous rotation effect at phases on both sides of 6 days—which 
is Gaposchkin’s secondary minimum. The system consists of a very luminous B star with strong absorp- 
tion lines and an equally luminous star which is associated with an emitting nebulosity of H and Het. 

HD 78316 (x Cancri).—A new orbit has been derived from 84 Yerkes single-prism spectrograms and 
3 McDonald Cassegrain plates. The elements agree well with those recently determined at Victoria 
except for an unexplained difference in y and a small difference in w. The spectrum contains unusually 


strong lines of Mn 11. 
HD 163181! 


This star is of special interest because it has bright lines of variable intensity, which 
have sometimes been likened to the bright lines in the spectrum of 8 Lyrae. A compari- 
son of the two stars is, however, at first quite disappointing. On the first few spectro- 
grams taken at the McDonald Observatory the bright lines were very weak; and, al- 
though on later dates they became somewhat stronger, they never even remotely re- 
sembled the conspicuous emission lines of 8 Lyrae. However, a more critical examination 
of the 55 spectrograms obtained between July 9 and September 15, 1943, shows a number 
of very interesting features which may help in elucidating some of the puzzles which 
abound in our discussions of Be stars. 

The binary nature of HD 163181 was discovered at the Mount Wilson Observatory, 
and an orbit was published in 1928 by Humason and Nicholson.? The period is almost ex- 
actly 12 days, and the range in velocity is 385 km/sec. The eccentricity is 0.065, and the 
longitude of periastron is 23.2°. Humason and Nicholson discovered the variability of 
the emission lines and also of the absorption lines and established the principal character- 
istics of these variations. Through their kindness I was able to examine several of their 
spectrograms. There seems to have been no appreciable change in the intensity or ap- 
pearance of the bright and dark lines. 

The star took an added importance when S. Gaposchkin,’ in 1938, announced that it is 
an eclipsing variable with two almost equal minima, whose depths are 0.43 and 0.42 mag., 
respectively. Since then the variable has also been observed by Hertzsprung? and his as- 
sociates, who find a 8 Lyrae-type light-curve and closely confirm the results of Gaposch- 
kin. I have adopted Gaposchkin’s epoch of principal minimum, JD 2428429.79, and his 
period of 12.004211 days and have computed my phases with the help of these data. The 
radial velocities are listed in Table 1. Separate columns give the mean results from all 
lines except those of interstellar calcium, from the calcium lines, and from the H lines. 


* Contributions from the McDonald Observatory, University of Texas, No. 86. 
1 V453 Sco; a = 17%49™7; 6 = —32°27’ (1900), mag. 6.6, Sp. BO. 

2 Ap. J., 67, 341, 1928. 

3 Harvard Bull., No. 909, p. 20, 1938; Ap. J., 89, 125 and 322, 1939. 


* Quoted by Mrs. C. P. Gaposchkin from a letter by Professor Hertzsprung to Dr. Harlow Shapley 
dated December 9, 1940. 
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RADIAL VELOCITIES OF HD 163181 


PHASE 
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| No. No. 
| of Cau of H 
| Lines Lines 
July | 0 2 
| 0 1 3 9 5 
| 2 3} 2 6 = 1 
2 2 6 0 8 
me 2 4 8 | 3 
| 3 2 3 1 7 : 
| 4 Me 2 4 0 | 3 
| 4 2 6 5 | 2 
5 sone 1 7 8 
5 1 3 7| 4 
| 8 py} 1 3 0 | 0 
1 8 | 0 
| 1 1 | 2 
| 2 1 3 7 | 3 
7 1 0 | 9 
7 2 2 | 9 
8 1 3 2 | 3 
Aug. 3 2 3 7 | 0 
4 2 5 2 
5 2 2 4 | 3 
5 2 2 2 | 9 
6 | 1 2 | 7 
6 1 0 | 9 
10 i 2 
11 3 | 5 
0 aes 2 1 | ) 
5 2 2 1 7? 
10 1 | 
10 fa: 2 5 
0 5 | 
0 ee | 2 2 | 7 
1 4 3 
2 1 | 2 8 
3 2 3 8 | ° 
Sept. 8 | 1 | 5 0 
10. 1 | 3 ‘3 
10| he: 9 
10| 1 ex. 2 
10. 1 | 3 8 
11] 2 a 5 
11| 2 2 8 
0. 1 2 2 
0. 1 3 
0. 1 9 
0. 1 2 A 
4. 1 2 5 
4. 1 | 3 “2 
5 1 | 3 3 
5. 1 3 4 
5. 1 he 8 
6. 1 A 
6 1 3 4 
6. 1 3 5 
7 1 3 
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The plates were grouped in 13 normal places as is shown in Table 2. The values of O—C 
are those which result from a least-squares solution. Table 3 shows a comparison of the 
velocities from the H lines with those from all stellar lines. The velocity curve is given in 
Figure 1, where the observations refer to the normal points. Because the period is so 
nearly 12 days, the observations are grouped at intervals of about 1 day, and there are 
no intermediate observations. The orbital elements were derived by means of Schlesing- 
er’s least-squares method and are shown in Table 4. The probable error listed in the line 


TABLE 2 
NORMAL PLACES 
No. | Limits of Phase | o-c 

0.08- 0.16 0.120 2 + 0.65 

0.90- 0.96 0.931 7 — 0.16 

1.12- 2.15 1.753 6 + 0.48 

2.92- 3.10 3.030 3 — 5.80 

3.91- 4.09 4.002 4 — 3.82 

4.90- 5.10 4.986 5 +10.92 

5.89- 5.96 5.935 6 +21.03\ Rotation 

6.89— 7.03 6.953 7 —11.17{ effect 

7.91-— 8.07 7.978 4 — 6.82 

10.00-10.01 10.005 2 + 6.05 

10.90-10.95 10.924 5 —10.26 

See 11.91-11.99 11.947 3 — 3.37 

TABLE 3 
COMPARISON OF H AND OTHER LINES 
Hi+H 
Normal Mean H —(All Normal Mean H All — 
Place Phase Plates Lines) —(All Lines) Place Phase Plates Lines) —(All Lines) 

0.120 2 —11.0 + 2.2 ol! 6.953 7 + 2.9 +25.9 
0.931 7 .—26.4 —26.4 7.978 4 + 7.0 — 0. 
1.753 6 —15.7 —18.6 8.900 1 + 6.6 +22.2 
3.030 3 — 2.1 — 8.3 10.005 2 —14.7 + 9.6 
eae 4.002 4 + 6.5 +13.3 re 10.924 5 + 4.4 +32.3 
4.986 5 — 9.5 —15.4 11.947 3 — 6.6 — 0.8 
5.935 6 — 5.6 + 3.4 


for y is that for Schlesinger’s unknown I’. The value of T is given in days counted from 
zero phase at JD2428429.79. The normal points were all given equal weight, in order to 
reduce the influence of physical factors, such as rotation, which tend to distort the veloc- 
ity-curve. The value of Av? was reduced from 1793 to 966. The probable error of a 
single normal point is +7 km/sec; that of a single plate is +8.7 km/sec. These values 
are large; Humason’s and Nicholson’s observations give an even larger value, namely, 
+9.2 km/sec for one plate. 

The agreement of my elements with those of Humason and Nicholson is very satis- 
factory. There has been no change in w, and the only significant difference is in e. The 
Mount Wilson value was 0.065. The observations show at once that there is a pro- 


NEW ORBITS FOR SPECTROSCOPIC BINARIES 213 


nounced rotation effect. The residuals are positive at phases near 5 and 6 days and nega- 
tive at phases near 7 days. The range in these residuals is of the order of +25 km/sec, 
so that the rotational velocity must be considerable. It would be premature to compute 
it without an accurate velocity-curve over the entire range of phases concerned. The 
effect is also unmistakably present in the Mount Wilson velocity-curve, though it was 


+200 
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Fic. 1.—Velocity-curve of HD 163181 


TABLE 4 
ORBITAL ELEMENTS OF HD 163181 


Preliminary Final Probable Errors 
y=— 42.8 km/sec y=— 39.8 km/sec + 3.5 km/sec 
K= 189.5 km/sec K= 192.6 km/sec + 2.9 km/sec 
w=+46.2° w=+23.7° +22.2° 
e=0.04 e=0.10 + 0.014 
T=4.78 days T=4.04 days + 0.76 day 


not originally interpreted in this way. With my elements the time of mid-eclipse should 
occur when v + w = 90°, or when the radial velocity is 


y + Ke cos w = —22.2 km/sec. 


The Mount Wilson elements give 
y+ Ke cos w = —30.3 km/sec. 


Whichever value we use, there is a strong indication that the spectroscopic rotation effect 
is slightly unsymmetrical with respect to this point. Both velocity-curves give positive 
residuals of considerable size at that phase at which the velocity is approximately — 26 
km/sec. The observations also suggest that the rotation effect begins at least 1 full day 
before mid-eclipse and lasts at least 1 full day after it. The limits may well be as great 
as +1.5 days. This does not accord well with Gaposchkin’s photometric result, D = 
0.09. P = 1.33 days, for the entire duration of eclipse. If the observations affected by the 
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rotation effect are excluded, the probable error of one plate is reduced to a value which 
is consistent with the internal precision of the measurements. 

Probably the most remarkable result of Gaposchkin’s work is the fact that his princi- 
pal eclipse falls 0.5P earlier than the eclipse of the B star which we have just discussed. 
Yet there is no indication in the spectrum of a second set of absorption lines which would 
correspond to that body which is eclipsed at phase 0, if we depend upon an extrapolation 
of Gaposchkin’s period, or at phase 0.6 day, if we again use the spectrographic criterion 
v + w = 270°. My observations in this respect agree with those of Humason and Nichol- 
son, who searched for the other component but could not find its absorption spectrum. 
This apparent paradox is of the utmost interest, and its solution will undoubtedly lead 
to important results. 

I believe that a satisfactory solution of the paradox is obtained from the following con- 
siderations. In the vicinity of my phase 0 or Humason and Nicholson’s phase 9.5 days, 
the absorption lines are very strong. This is precisely the phase of Gaposchkin’s principal 
mid-eclipse. The B star stands in front, the unknown star is eclipsed. The absorption 
lines, which are weak at all other phases, are greatly enhanced, and the amount of this 
enhancement appears to be very closely the same for all lines. The critical test comes 
from the absence of any dilution effect in the lines of Het or in the intensity of Mg 1 
4481. These lines are normal. The H lines are strong and narrow, but they lack the typi- 
cal sharpness of H lines produced in a tenuous shell. We conclude that we are observing 
the absorption spectrum of a star, not that of a shell. It is clear that this star must be 
the normal B-type component of the binary system whose absorption lines, at all other 
phases, are weakened by the superposition of the strong continuous spectrum of the un- 
known component. This latter, then, does display a continuous spectrum, even though it 
does not show absorption lines at the times of elongation. At our phase 0, when the ab- 
sorption lines of the B star are strong, the emission lines are very weak or absent. We 
conclude that the unknown star has emission lines in its spectrum and that the gases giv- 
ing rise to them are eclipsed at phase 0. 

At phase 6 days, when Gaposchkin observes his secondary eclipse, the unknown star is 
in front, and the emission lines are little affected. The absorption lines of the B star are 
weak and show the rotation effect, but they never quite disappear. Since no absorption 
lines belonging to the unknown star are observed at the elongations—in spite of the fact 
that Gaposchkin finds for its total light Z, = 0.46 as against L; = 0.54 for the B star— 
we must conclude that the unknown star is not responsible for the absorption lines at 
phase 6 days. The eclipse is probably not total, though we cannot be quite certain of this 
because of the gaps in our observations. It may be annular, because at phases near 0 days 
the emission lines do seem to disappear completely. 

The general run of the values O—C in Table 2 shows that there is no rotation effect on 
the two sides of phase 0. The B star is in front, and its lines should, of course, not show 
such an effect. But the bright H lines disappear at phase 0. Hence we might expect a ro- 
tation effect in the measures of the H absorption lines, if the latter are to any extent pro- 
duced by absorption of the continuous light of the unknown star in the emitting shell. 
Table 3 might suggest that such an effect does exist: for Hé and Hy we have pgsitive 
residuals at phases 8.9, 10.0, and 10.9 days, and negative residuals at phases 0.9, 1.8, and 
3.0 days. During the larger part of this interval the emission lines are very weak or ab- 
sent, so that there can be no serious trouble from blending. But the H lines which we 
measure are those of the B star, not of the unknown star. Hence the phenomenon, if it is 
real, is not a simple rotation effect. It is more probable that the H lines of the B star are 
blended with those of the unknown star and that the residuals of the H velocities reflect, 
very slightly, the motion of the unknown star. The question is whether or not we can be 
certain of any absorption lines of the unknown star. At phase 8.0 days the emission lines 
are strong and are greatly displaced toward the red. They are single, without central re- 
versal, so that, if there is any absorption line due to the unknown star’s shell, it must 
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blend with the H lines of the B star and shift it toward the red. But at phase 8.0 days the 
residual is nearly zero: there is no absorption from the shell. At phases near 6 days the 
emission lines are double. But the eclipse of the B star is presumably annular, or partial, 
so that the H absorption lines may belong to it. The emission lines are too weak to show 
whether the absorptions are rendered shallow by the superposition of emission or are cut 
below the continuous spectrum by high-level gases. We must, therefore, leave undecided 

the origin of the H residuals in Table 3. 

Our picture is complicated by the fact that the emission lines are quite dissimilar at 
the two elongations. At phases around 3 and 4 days they are very weak and shifted to- 
ward the violet of the absorption lines. At phases around 9 days they are strong and shift- 
ed conspicuously toward the red. There can be no question that the velocity of the emis- 
sion line taken as a whole shifts, relative to the absorption line, with a phase difference of 
(.5P. Whether the emission line also shifts with respect to the center of mass of the sys- 
tem is not quite certain because it is exceedingly difficult to measure the bright lines with- 
out being influenced by the disturbing presence of the absorption lines. It is my impres- 
sion that the emission lines do shift but that their range is considerably smaller than that 
of the absorption lines. This does not mean, however, that the unknown star is the more 
massive one, because it is possible that the velocity of the emission lines does not repre- 
sent the velocity of the unknown star. 

Why, then, is the appearance of the emission lines so different at the two maxima? 
I suggest that the unknown star has a continuous spectrum of approximately the same 
intensity as the B star and that it has very weak absorption lines which are never seen at 
the elongations or at secondary minimum. It does, however, have associated with it a 
set of emission lines whose velocities are systematically shifted toward the red by perhaps 
1.5 A, as are the emission lines of 29 Canis Majoris or of the Wolf-Rayet component of 
HD 193576.5 In 29 Canis Majoris the red shift occurs in the brighter component of the 
binary, and the displaced emission lines oscillate in phase with the absorption lines. In 
HD193576 and HD 163181 the displaced emission lines belong to one of two approxi- 
mately equally luminous stars, and they oscillate with a phase difference of 0.5P from the 
absorption lines. It is quite probable that the red shift of the emission lines is a fairly 
general phenomenon among spectroscopic binaries of very early type. We cannot ex- 
plain it, but we are fully aware that it constitutes an important problem in astrophysics. 
It is almost certainly related to the red shift observed by C. S. Beals® in the emission 
lines of the P Cygni-type star, MWC 374. In this remarkable spectrum the shift is prob- 
ably caused by self-absorption in a decelerated expanding shell. There is, at present, no 
evidence of expansion in the gaseous shell of the unknown component of HD 163181. 
Even if such expansion does exist, it would be difficult to avoid the complication of a vio- 
let shift, which should probably be present in an expanding shell that is small enough 
to fit within the rather small amount of space allowed by Gaposchkin’s revised photo- 
metric elements. These considerations, together with O. C. Wilson’s results for HD 
193576, leave the problem of the red shift unsolved, though I should not want to go so 
far as to question, with Wilson, the expanding-shell hypothesis in the case of a normal 
Wolf-Rayet star. Self-absorption in the manner considered by Beals for MWC 374 un- 
doubtedly is important in some stars. Unfortunately, neither HD 163181 nor HD 193576 
has forbidden lines. If such lines were present, their velocities should not be affected by 
self-absorption, and the test would be easy. 

The H emission lines of HD 163181 are strong at phases 7, 8, and probably 9 days. 
But at phase 10 days they are much weakened. They remain about constant in intensity 
until phase 11 days, being then clearly double, with the stronger component on the red 
side. At phase 0 they are double and very weak. But their presence—enhanced, perhaps, 
by the longer exposures necessary at minimum light—suggests that the H shell is never 


50. C. Wilson, Ap. J., 95, 402, 1942. 6 J.R.A.S. Canada, 37, 241, 1943. 
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totally eclipsed but remains as a luminous fringe around the totally eclipsed unknown 
star. At phases 1, 2, 3, and 4 days the emission lines are very weak, presumably because 
of the general red shift, which would bring their strongest portion roughly i in coincidence 
with the absorption lines of the B star. It is difficult to decide whether the intensities of 
the latter are consistent with this interpretation, because of the weakness of the emission 
lines. But the weakening of the emission lines at phase 10 days, only 1 day after elonga- 
tion, would then indicate that the eclipse of the emitting shell of the unknown star is al- 
ready well under way 2 full days before mid-eclipse. It is entirely possible that the emit- 
ting gases are not symmetrically distributed around the unknown star, asl it is not possi- 
ble to separate such an effect from the proposed red shift. 


HD 78316 (76 x CANCRI)? 


The binary nature of x Cancri was discovered in 1904 by Frost and Adams from 5 
spectrograms taken at the Yerkes Observatory. An orbit, based upon 25 Yerkes plates, 
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Fic. 2.—Velocity-curve of 76 « Cancri 


was obtained by N. Ichinohe® in 1907. The eccentricity of this orbit was 0.15, and the 
velocity-curve was sufficiently unsymmetrical to place the star among those objects for 
which a new determination is desirable in order to test whether there has been a change 
in the line of apsides during the last 36 years. A new series of single-prism spectrograms 
was started at the Yerkes Observatory in 1929 and was continued, with interruptions, 
until May, 1943. In the meantime J. A. Pearce and Phoebe Riddle® have published the 
results of their measurements of 32 high-dispersion spectrograms of « Cancri taken be- 
tween March and May, 1938, with the Victoria 72-inch reflector. The dispersion used for 
this work at Victoria was much greater than that used at Yerkes, and the precision of the 
Victoria orbit is greater than that obtained by me. Nevertheless, there are certain fea- 
tures in my velocity-curve which appear to justify its publication. In the first place, the 
Yerkes observations have a gap between 1935 and 1941 and therefore do not duplicate in 
time the observations made at Victoria. In the second place, my velocity-curve differs 
appreciably in shape from that computed with the Victoria elements, and a small differ- 
ence in w is indicated. Since this difference is not consistent with the fact that the value 


Tq =. 9h2m3; § = +11°4’ (1900), mag. 5.1, Sp. B8. 
8 Ap. J., 25, 315, 1907. ® Pub. A.A.S., 10, 65, 1940. 
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obtained by Pearce and Riddle is almost identical with the value obtained by Ichinohe, 
there is no evidence of a rotation of the line of apsides. But Ichinohe’s orbit is based upon 
a relatively small number of plates, and, while his mean error per plate is only +4.2 
km/sec, the minimum of the velocity-curve was not observed sufficiently to allow a very 
critical determination of w. 

The 84 new Yerkes observations and 3 McDonald observations are listed in Table 5. 
The phases were computed with the period obtained at Victoria, and zero phase is the 
Victoria value of the time of periastron passage. The observations were arranged in 15 


j TABLE 5 
RADIAL VELOCITIES OF 76 « CANCRI 


Plate Date U-T. | Phase |(Kv%eq)|| Plate Date U.T. | Phase ec) 

1R9180.....| 1929 Apr. 14 | 350m] 24450 | 451.1 || 1R13040....| 1942 Dec. 12 | 11»36m| 64108 | —23.6 
9188... Apr. 15 | 244 | 3.404| 478.2 || 13049...) Dec. 16 | 1003 | 3.651 | +80.4 
9197... Apr 16 | 446 | 4.488| 457.6 || 13050... Dec. 16 | 1105 | 3.694| +79.2 
9995... Dec. 19 | 820 | 3.483| +80.9 || 1943 Jan. 15 | 706 | 1.562 | +413.9 
10005....| 1932 Jan. 19 | 635 | 2.444| +57.1 || 13057 an. 15 | 810 | 1.606] + 4.4 
10884...) 1933 Dec. 6 | 1135 | 5.585 | —0.3|| 13078 eb. 16 | 710 | 1.598] +10.1 
11663....| 1935 Mar.11 | 5 36 | 5.028] +27.0|| 13079 Feb. 16 | 814 | 1.643] 413.2 
11668 Mar.13 | 411 | 0.575 | —55.7 || 13089 Feb. 18 | 525 | 3.526] +82.1 
11672 Mar. 15 | 443 | 2.597] +52.6 || 13090 Feb. 18 | 685 | 3.574] +78.3 
11678 Mar.16 | 549 | 3.643/ 472.4 || 13091 Feb. 18 | 752 | 3.628] +87.8 
11687 Mar.24 | 350 | 5.168| +25.6 || 13092 Feb. 18 | 851 | 3.669] 485.3 
11694 Mar.27 | 359 | 1.781] 421.4 || 13097 Feb. 19 | 538 | 4.535] +61.7 
11702 Mar.29 | 347 | 3.773| +474/5 || 13098 Feb. 19 | 651 | 4.586] 452.1 
11705 Apr. 7 | 350 | 6.382 | —59.9 || Mar. 2 | 759 | 2.847| +74.3 
11707... Apr. 13 | 3.35 | 5.978 | —22.7 || Mar 4 | 805 | 4.851] +465 
1941 Nov 18 | 1134 | 0.696 | —52.3 || Mar.13 | 811 | 1.068] —25.4 
12654... Nov. 20 | 1146 | 2.704 | || Mar.21 | 316 | 2.470] +59.9 
12659 Dec. 6 | 956 | 5.842] — 4.8 || 13109... Mar.21 | 425 | 2.518] +53.8 
12669... Dec. 11 | 1010 | 4.459] || 13110. Mar.21 | 540 | 2.570] +56.7 
12676... Dec. 14 | 1118 | 1.113| —23.2 |] Mar.21 | 643 | 2.614] 467.7 
12700... 1942 Jan. 6 | 755 | 4.792] 441.0 13117222) Mar 22 | 448 | 3.534| +84.6 
12707... Jan. 13 | 913 | 5.453] + 4.7 || 13118... Mar.22 | 554 | 3.580] +78.3 
12717 Jan. 23 | 800 | 2.617] +4001 13119... Mar.22 | 657 | 3.624] +82.1 
12724 Jan. 24 7 45 3.606 +79.0 13126.... Mar. 30 4 08 5.113 +24.1 
Feb. 12 | 820 | 3-450) 72.7) Mar.30 | 545 | 5.181| 
12736 Feb. 13 | 604 | 4.356] 458.1 || 13130... Mar. 31 | 334 | 6.089| —33.2 
Feb. 20 | 6:35 | $30.6 |] Mar.31 | 445 | 6.139] —41.5 
12747... Feb. 21 | 558 | 5.958) —41.1 || 13132... Mar.31 | 549 | 6.183] —43.2 
12749....| Mar. 5 | 552 | 5.168] 421.7 || 13134... Apr. 1 | 413 | 0.724] 
12752... Mar.19 | 405 | 6.307| —49.8 || 13135... Apr. 1 | 515 | 0.767 | —42.2 
12753...  Mar.19 | 458 | 6.344] —44.2 || 13142... Apr. 8 | 237 | 1.264] — 8.9 
12757....,  Mar.30 | 402 | 4.519] |] 13143... Apr. 8 | 339 | 1.307] 15.3 
12759....| Apr. 1. | 355 | O.121| |] 13144... Apr. 8 | 453 | 1.358] — 3.4 
12766 | Apr. 12 2 55 4.686 +52.4 13154.... Apr. 21 1 43 1.439 ~ 1.4 
1272 Apr. 13 | 405 | 5.735] —13.6 || 13155... Apr. 21 | 252 | 1.488] + 3.2 
Apr. 15 | 3.25 | 1.314] — 2:5 || 13156... Apr. 21 | 356 | 1.532] +4140 
12794 | Apr. 29 2 00 2.468 +53.5 13162.... Apr. 22 1 40 2.438 +47.1 
12976... Nov. 6 | 1147 | 2.082] 443.1 || Apr. 22. | 243 | 2.481] +57.2 
12991 Nov.17 | 1046 | 0.252] —51.3 || Apr. 22 | 3.46 | 2.525] +54/9 
Nov.17 | 1148 | 0.205| —52.8|| ay 1 | 145 | 5.048] 437.1 
Nov. 27 | 1134 | 3.893) 477-3 |] 1 | 248 | 5.092) 432.9 
13019 Dec. 4 | 1043 | 4.464} 466.3 |] 131795.) May 3 | 218 | 0.678| —S0.1 
13020 Dec. 4 | 1150 | 4.511 | 02.2 || May 3 | 345 | 0.738| —47.6 
13039 Dec. 12 | 1031 | 6.063} —39.1 | 


normal places (Table 6), and the values of O—C were computed separately for each plate 


with the help of the preliminary elements and were then averaged for each normal place. 
A least-squares solution was made following the method of Schlesinger and adopting 
equal weights for all normal points. This was done in order to minimize the possible ef- 
fect of physical causes which are sometimes known to distort a velocity-curve. The value 
of 2Av* was reduced from 123.8 to 63.6, and the probable error of one normal velocity is 
+1.6 km/sec; that of one plate is of the order of +3.5 km/sec. The elements are shown 
in Table 7. The time of periastron passage is expressed in days counted from phase 0, or 
from JD 2429012.081. The probable error listed in the second line is that obtained for 
Schlesinger’s variable I’. The principal difference from the elements of Pearce and Riddle 
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is in the value of y, for which they find +25.3 km/sec. It is possible that this difference 
of 4.2 km/sec is due to a systematic difference between the two instruments, though the 
amount is alarmingly large. The published systematic differences for stars of classes B8- 


A3 are: 
Yerkes — Lick = +0.95 km/sec., 


Victoria — Lick = —0.95 km/sec .!° 


However, these determinations (which incidentally are in the wrong sense) probably do 
not apply in the present case. In September, 1941, a change was made in the Bruce spec- 


TABLE 6 
NORMAL PLACES 


Limits of Mean | No. of Limits of Mean | No. of 
No Phase Phase Plates vated, No. Phase Phase | Plates 
0.12-0.30 0.22 3 —53.8 +2.3 3.65-3.89 3.73 5 +77.5 +2.1 
» 0.58-0.77 0.70 6 —44.7 -—2.3 ) 4.36-4.59 4.49 8 +60.4 +0.3 
1.07-1.36| 1.24] 6 —13.61 40.5 || 11...... 4.09-4.85 | 4.78] 3 +47.8| —1.3 
1.44-1.64 1.55 7 + 7.3 +0.7 4.98-5.18 5.10 8 +30.0 —2.6 
1.78-2.08 | 1.93| 2 431.0] +1.5 13...... 5.45-5.84| 5.65| 4 +3.6 
2.44-2.52 2.47 8 +57.2 —3.1 5.96-6.18 6.07 7 —35.0 +0.1 
2.57-2.85 | 2.66| 6 464.1] —1.4]] 15...... 6.31-6.38 | 6.34) 3 —47.8| —3.4 
Sor 3.40-3. 64 3.55 11 +78.2 +1.0 
TABLE 7 
ORBITAL ELEMENTS OF 76 « CANCRI 
Preliminary Final Probable Errors 
y=+21.1 km/sec y=+21.1 km/sec +0.8 km/sec 
K= 66.9 km/sec K= 66.1km/sec +0.6 km/sec 
w=161.0 w=174.6° +2.4° 
e= 0.14 e= 0.137 +0.086 
T= 0.00 T=+0.21 day +0.065 day 


trograph which resulted, according to Hynek’s measures of ¢ Tauri, in the following dif- 


ference: 
Yerkes new — Yerkes old = —2.0 km/sec. 


Measures of several plates of standard-velocity stars have given a similar result. Since 
formerly the systematic error of the Yerkes plates was +1.0 km/sec, it should now be 
about —1.0 km/sec. It is difficult to account for the rest of the discrepancy. The plot of 
the velocities in Figure 1 certainly does not warrant the idea that there has existed since 
1941 a systematic difference of —4 km/sec, which did not exist in 1929-1935. If there is 
a difference, it is in the opposite sense. It is not impossible that the velocity of the center 
of mass of x Cancri is slightly variable. This question must be postponed until more ma- 
terial on the systematic errors is available. 

In spite of this uncertainty, the other elements are known with considerable precision. 
There is no spectrum of the secondary component, and it is probable that what Ichinohe 


10 J. H. Moore, Pub. Lick Obs., 18, xi, 1932. 
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TABLE 8 
LIST OF LINES IN THE SPECTRUM OF x CANCRI 
Wave Length Int. Identification Wave Length Int. Identification 

3289.40...... 0 Fe 9.35 3 7.53 
in Cru 5.42, Fei 5.81 255 1 Fe 1 9.88 
Fe 11 7.89 3507 .93 0 Fett 7.39 
$305.19 ..... 1 Fe 2.86, Fe 11 3.47 1 Ti 11 0.85 
1 Fe 11 14.00 1 Fem 5.82 
1 Tit 2.93, Fe 11 3.07 1 Tim 5.41 
1 Tit 6.76, Cr tr 8.35 

3585.90... 0 Crit 5.31 
$335.12...... O T6145 Cr15.28,|| 3627.55...... 1 | Few7A7 

Cru 5.46 1 Cr 1.49, Crit 1.72 

2 Cr 11 6.32 Onn | H 202.81 
1 Tim 7.85 Ti 5.20 
0 Cru 9.80, Ti 0.34 1 H 19 6.83 
GOAN 0 Titi 1.84 inn | A 181.56 
3349.02...... 9.00, Tiu 9.41 || 3697.11...... 2nn | 177.15 
0 Fe 1 6.26 3704.03 3nn | 16 3.86 
Cr 8.50 In Cat 6.04, Ti 11 6.23 
3360.05 1 Cr 11 0.30 3708 .40...... 1 Mn tt 8.06 
1 Ti 6.18 4nn | A 15 1.97 
3368.03...... In Cr 11 8.05 1 Crit 5.19, Cru 5.45 
0 4.35 5nn | H 14 1.94 
| Cr 11 9.82 7nn | A 13 4.37 
$367 .67...... 2 Ti 7.84 3 Ca tr 6.92 
In Ti u 4.55, Crit 4.32 2 Fet 3.37 
... 1 Cru 9.54 10nn | 120.15 
S808. 1 Crit 3.32 2n Ti tr 9.30 
10% 1 Cru 8.77 In Ti 1.32 
3419.47...... 1 (Fer 8.51) 4+ Fet 3.80, Fem 4.11 
1 Crit 1.20 3165.84... ..... 0 Crit 5.62 
3422 .83 Crit 2.74 0 Fei 7.21 
In Cr 11 8.94 |) 10nn | 110.63 
0 Cr 11 3.30 0 Y 4.34, Tim 4.65 
0 Fem 6.11 2 V 8.36 
5 Mn 8.96 0 Fem 1.51 
3442.04...... 7 1.98 10nn | H 107.90 
1 Co 11 6.40 5012. 2 Fei 2.96 
3451.69...... 1 Fei 1.23, Fem 1.32 1 Feu 4.12 
18% In Cr 7.62 1n Feu 1.97 
7 0.33 1 Fem 4.91 
3464.45 Onn | Ferm 4.50 3827.69...... 1 Feu 7.08, Fe 1 7.69, 
1 Fe 11 8.68 | Fe17.83 
0 Niu 1.35 3832.49...... 0 Mg12.31 
7 Mn u 4.06, Mn 114.15 || 3835.52...... 10nn | 95,39 
0 Crit 5.13, Fem 5.25, || 3844.35...... + Mn it 3.98 

Feu 5.74 0 V 17.32, Mg 8.24 
3 2.91 5 Sit 6.03, Mn 6.53 
3488.79...... 4 Mn tt 8.68 3862.84...... 5 Sit 2.60 
3490.42...... 0 Fe1 0.60, Ti 1 1.06 3878.08...... 1 Fet 8.02 
3492.48...... 3679: 3 (Fe 8.58) 
3494.41...... 0 Feu 4.67 <> ee | 10nn | A8 9.05 
3495.80...... 3 5.84 3896.00...... | Onn | V1 6.16 


219 
_| 
ce 
be 
of 
ce 
is 
er 
a- 
n. 
he 


220 OTTO STRUVE 


TABLE 8—Continued 


Wave Length Int. Identification Wave Length Int. Identification 
3896 06. ..... 2n Fe1 7.89, 8.01 || 4192.51...... 0 Niu 2.02 
3900.68...... 2n Tin 0.54, Al 0.68 0 Sim 8.25, Fer 8.31 
1n Fe1 2.95 AZOD 49... in Tiu 0.40, Mn 1 0.25 
3905.56...... in Crit 5.64, Sir 5.53 4202.65...... 1n Fer 2.03, V 11 2.34 
3906.27. <2... 0 Feu 6.04 4905.49... | 2 V 11 5.09, Fe 1 5.48 
0 Tit 3.46 420663. 3 Mn 6.43 
0 Feu 8.51 in? S52 
2 Fez 0.26 0 Fet 9.36 
1 Mnu 6.47 1 2.22 
¢ 7 Cau 3.67 3 Feu 3.17 
1 Feu 5.94 A285. 22... 0 5.94 
1 Fe10.82 2 Mn 9.72 
0 Fe1 2.44 1 Alu 0.75 
3 Mnu 3.82 4242 .48......| 3 Cru 2.38, Mn 1 2.38 
3958 .58...... 0 Cr 8.07 4244 .36...... 4.25 
396838... ... 3 Ca 8.47 1 Sc 11 6.83 
3969.98...... 10nn | H €0.08 4248.17...... 1 Mn11 7.95, 8.23 
0 Fe1 1.77, Fem 1.61 | 2n 1.77 
3983 .86...... 5 Fet 3.96, Cr1 3.91 | 3 Crit 2.62, Mn 3.02 
£002.30. ...<. 0 Fem 2.55 42590.38...... 3 9.26 
4013.82... .... 0 Fet 3.82, Fe1 3.79 07... 0 Feu 3.89 
A03426...... 2n Feu 4.55 226 3n Car 
3nn | Her6.19 Fe1 6.97 
415...... 0 Tiu 8.34 4269.26...... 0 Cru 9.28 
1 Fe1 0.49, Mn 1 0.75, 1 Feu 3.32 

Cru 0.28 0 Cru 5.57 
2 Fe 1 2.95 1 Fe 8.13, Crm 8.10 
4034.38...... 0 (Mn 4.49) AS... 3 Mn 11 2.50, Fer 2.41 
4038 .98...... 1 Mn 8.73 4284.26...... 1 Cru 4.21 
4044 .32...... 1 Feu 4.01 .. 0 Tiu 7.88 
4045.80...... 0 Fer 5.81 4202 3 Mn 112.25 
0 Fe 3.60 3 Tit 0.05, Mn 11 0.22 
4066.93...... 0 Nim 7.05 4303':43...... 2 Feu 3.17 
0 Cr 11 0.90 4305.89...... 0 Fe1 5.45, 5.71 
0 Sim 5.45, Crm 5.63 2n Ti 7.90, Fer 7.91 
10nn | 61.74 | 2 5.76 
4104.75...... 0 Mnu 5.01 3 Mn 11 6.76 
..... 2 Cru 1.01 4340.58...... 10nn | H y 0.46 
74... 1n Het 0.81 4344.17....... 4 Ti 4.29, Mn 4.03 
6 Si18.05 1 Mnu 8.49 
4130.89...... 6 Si 11 0.88 0 Ti 110.83 
4136.94...... 4 Mnu 6.91 4351.94. .. 1 Feu 1.76 
0 Mn11 0.19 in Fem 7.57 
0 Het 3.76, Fe1 3.87 4363..39...... 0 Crit 2.93 
4163.89...... On Ti 3.65 S77 1n (Mo 7.76) 
1n Tiu 1.90 1 Mn 9.74 
2n Feu 3.45 1 Feu 5.38 
1 Fe1 7.60, Y 7.55, 0 7.93, Fe 7.90 
Fe 7.70 4409 .46...... 1 Ti 9.22, Tim 9.52 
4978.78... .. 3 Fem 8.85 4414.66...... 0 Fer 5.12 
On V 3.43 1 Fe 11 6.82 
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TABLE 8—Continued 


Wave Length 


Wave Length Int. Identification Int. Identification 
0 Feu 1.63 0 (Cr 111.82) 
4434.10...... 0 Mg itt 3.99 1 Feit 5.34 
4451.66..... 0 Few 1.54 4519.28...... 1 (Fe 1 0.22) 
4406.11...... On Fe1 6.55 Comes 
4468.20...... In Ti 11 8.50 4528.33. 0 8.62, V 11 8.51 
4471.68...... 3n Het 1.48, Het 1.68 4549 .45...... 3 Feu 9.47, Ti it 9.63, 
447865. 2 Mn uw 8.74 4555.68... . =... 1 5.89 
4481.48...... 6 Mg 1.33 1 Cr 8.66 
1 Fe 1.40 4566.42... 1 Cru 8.22 
0 Feu 9.71 4590.24...... 1 Tim 9.95, Cr 9.89 


thought was the spectrum of the secondary was really due to faint lines of Mn m1, belong- 
ing to the primary, which were not identified in 1907. There is, in particular, a line of 
Mn wat \ 4478.74, which looks deceptively like a component of Mg 11 4481.33. I have 
measured a considerable number of lines on two McDonald Observatory spectrograms 
on Eastman Process emulsion. The list of these lines is given in Table 8. The appearance 
of the spectrum is that of a typical Mn 1 star with sharp lines. 
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THE SPECTROGRAPHIC PROBLEM OF U CEPHEI* 
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ABSTRACT 


The velocity-curve of this eclipsing binary was redetermined from measures of 115 McDonald Ob- 
servatory spectrograms. If this curve is interpreted in the usual way, the elements are: P = 2.493 days, 
y = —5 km/sec, K = 120 km/sec, e = 0.20, w = 40°, and T = phase 2.16 days after minimum light. 
Although the eccentricity is much smaller than that obtained by Carpenter from plates taken nearly 20 
years ago, the asymmetry of my velocity-curve does resemble his and confirms the discrepancy between 
spectrographic and photometric observations. In the partial phases of the eclipse the lines of the B8 
star show a large rotation effect which suggests that the equatorial velocity of rotation of the brighter B8 
star is 200 km/sec. That of the fainter, G2, star is less than 50 km/sec, in spite of the greater size of this 
star. The rapid angular velocity of the B8 star, which is far in excess of the angular velocity in the orbit, 
constitutes a particularly interesting cosmogonical problem. During the partial phases there appear 
sharp lines of H and Ca 1 which are displaced from those of the limbs of the B8 star and which give ve- 
locities of about +100 km/sec before totality and —150 km/sec after totality. These lines are not pro- 
duced by the G2 star but must originate in a tenuous stream of gas which may circle around the G2 star 
in a manner similar to that of the streams in 8 Lyrae and SX Cassiopeiae. The spectrum of the B8 star 
shows systematically weaker lines in the first partial phase than in the second. Similarly, the lines of the 
receding stream (before totality) are weaker than the lines of the approaching stream (after totality). 
The lines of Het confirm Redman’s conclusion that lines broadened by Stark effect must be weaker at 
the limb of the B8 star than in the center. On the other hand, Ca 11 does not appreciably change in 
equivalent width and therefore becomes more conspicuous (narrower and deeper) as totality is ap- 
proached. This constitutes an argument against the hypothesis that “pure absorption” is predomi- 
nantly responsible for the formation of lines in the B8 star. It is suggested that the asymmetry of the 
velocity-curve may be caused by the same stream which gives rise to sharp H lines in the partial phases. 


1. The problem.—The eclipsing binary U Cephei' presents one of the most intriguing 
puzzles in stellar spectroscopy. Its light-curve unmistakably points toward an orbit in 
which the value of e cos w is close to zero, because its shallow secondary minimum is al- 
most exactly halfway between its deep primary minima. On the other hand, Carpenter? 
found from 34 Lick Observatory spectrograms that the velocity ‘curve is very unsym- 
metrical, giving e = 0.474 and w = 25.0°, so that 


ecosw=0.43. 


If this conclusion from the spectrographic results were correct, the secondary eclipse 
would occur 0.51 day or 12.2 hours later than has actually been observed with the pho- 
tometer. The uncertainty with which the secondary eclipse has been placed at phase 
0.5P probably does not greatly exceed 0.1 day. We are, therefore, confronted with a 
serious contradiction between the photometric and the spectrographic results. Most 
photometric observers have, perhaps somewhat arbitrarily, assumed that e = 0. The 
depth of the secondary eclipse in visual light is only of the order of 0.1 mag., and in violet 
light it is even less. Hence, the shape of this shallow minimum cannot be accurately de- 
termined; and there is, at present, no reliable photometric method to determine for this 
star e sin w, from which e and w could be separately obtained. Nevertheless, the period 
is so short (2.493 days) that the assumption e = 0 must certainly be regarded as a close 
approximation to the truth. Among the photometric investigators, only two have at- 


* Contributions from the McDonald Observatory, University of Texas, No. 87. 
1 HD 5679; a — 0553™4, 6 = +81°20’ (1900), mag. 6.85-9.19, Sp. HD AO. 


2 Ap. J., 72, 205, 1930. 
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tempted to determine e and w, although others* had recognized a displacement of the sec- 
ondary minimum by about 0.1 day. Completely at variance from all other light-curves 
is the result of J. Ellsworth,* who obtained a displacement of secondary minimum by 
+0.3 day from phase 0.5P and, with the help of Carpenter’s radial velocities, computed 
e = 0.221, w = 15°. It is doubtful whether this result is correct, because Walter® and 
Carrasco,’ who observed U Cephei at nearly the same time as Ellsworth, failed to get any 
noticeable displacement of secondary minimum. More recently the entire problem has 
been rediscussed by N. Tchudowitchev.’ Using his own photographic observations and 
those of all preceding workers, he represented the observed epochs of primary minimum 
by means of an expression of the form 


O-C=AT, +aP-E—*~ cos (1) 


where A7> is the error of the initial epoch (determined by Wendell), AP is the error of the 
period (also determined by Wendell), ¢ is the eccentricity of the orbit, w is the longitude 
of the periastron, and Aw is the change ir. during one revolution. The last, periodic, 
term results from an assumed slow rotatien of the line of apsides in a slightly eccentric 
orbit.’ The observations gave 


O— C= —0.0050 + 0.0000165E — 0.063 cos (0.028°E + 235°) . (2) 


This led to the following values: 
€=0.080; w=99°(forE= 8000); Aw=0.028°. 


The determination rests upon several unreliable old observations and may be seriously in 
error. The most recent observations, upon which the present investigation is based, give 
an appreciable departure in O — C from formula (2). Nevertheless, the general idea is 
probably correct. Tchudowitchev shows that qualitatively, at least, the advance of the 
line of apsides by 0.028°£ per revolution, or by 27 in 99 years, is consistent with the fact 
that the earlier light-curves by Wendell (1895-1912), Dugan® (1914-1917), and Baker ° 
(1914-1916) show a positive shift of secondary minimum of the order of +0.1 day, while 
the most recent curves by Walter (1931-1933) and Tchudowitchev (1934-1935) show a 
small negative shift of the order of —0.03 day. This could be interpreted as an advance 
in w from about 0° to 90° in an interval of a quarter of a century. The photometric data 
require an eccentricity of e < 0.1, which is entirely inconsistent with Carpenter’s veloc- 
ity-curve. 

The question arises whether Carpenter’s velocity-curve is in error or whether it is dis- 
torted by other effects. My discussion of the radial velocities of SX Cassiopeiae'® has 
shown that such distortion can result from currents in the outer atmosphere of a binary 
system which resemble the currents first found in the system of 8 Lyrae. The only other 
investigation of the radial velocity of U Cephei known to me is a study by Z. Kopal" 
of a series of 34 slit spectrograms obtained at Victoria by R. O. Redman for an investiga- 


8 R. H. Baker, Bull. Laws Obs., 2, 257, 1921. 

4 Pub. Lyon, 2, fasc. 1, 82, 1936. 

5 A.N., 250, 333, 1933. 

6 Boll. Obs. Madrid, 1, Nos. 12-13, 1933. 

7 Bull. Engelhardt Obs., Kasan U., No. 17, 1939. 


8 This term was first used by Tisserand to explain a periodic term in the distribution of minima of 
Algol, C.R., 120, 125, 1895. 


® Contr. Princeton U. Obs., No. 5, 1920. 
10 Ap. J., 99, 89, 1944, Harvard Bull., No. 914, 1940. 
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tion of the center-to-limb variation of the absorption lines during the partial eclipse. 
Kopal measured 5 spectrograms taken in full light between phases 0.921 and 1.057 days 
(counted from minimum light), 15 spectrograms taken during the partial phases, — 0.234 
to +0.137 day, and 10 spectrograms taken during the time of totality when the spec- 
trum was that of the faint secondary. From these latter 10 observations, with the help 
of an additional velocity obtained by Joy, Kopal derived the values 


y=+7.341.0 km/sec and K,=200+11 km/sec, 


under the assumption that e = 0. Using the value of y thus obtained, and again taking 
= 0, he next derived from the 5 plates taken at almost the same phase 


K,=62+6 km/sec. 
This value is completely different from Carpenter’s result, 
K,= 109.9 + 3.1 km/sec . 


Kopal’s value was determined by means of the formula 
Vi= K, sin (M +7) 


where ¥ is taken from the observations of the secondary and where M (the mean anom- 
aly) has practically the same value for all 5 observations. The mean phase of these ob- 
servations was 1.0 day, which is only 0.225 day removed from inferior conjunction. 
Hence, a small systematic error in the mean velocity measured on the 5 plates will be 
magnified by a factor of 2 in the determination of Ay. Similarly, an error in y, which in- 
cidentally is very sensitive to the precision with which the phases are known in Kopal’s 
work, will have a large effect upon K;. Carpenter found y = —6.0 + 3.7 km/sec. The 
difference of 13 km/sec between the two results is somewhat too large to be ignored. 
The probable error found by Kopal for y is probably too small—perhaps partly because 
of an uncertainty in establishing the exact time of central primary eclipse. This latter 
epoch was computed by Kopal with an average period of 2.492934 days, which he had 
found to represent satisfactorily 48 observations of minimum light covering an interval 
of 20 years. The period of U Cephei is believed to show appreciable fluctuations in a 
rather complicated manner, so that, according to Carrasco,” two sine terms are needed, 
in addition to a term containing £ and another containing E*. Experience has shown that 
errors amounting to 1 hour, or 0.042 day, may occur when a constant period is used over 
an interval of 20 years. A shift in the epoch of minimum by 0.04 day would change vy by 
almost 20 km/sec. However, Dr. Kopal has advised me that for the Victoria observa- 
tions there could have been no appreciable error in the time of minimum light. Hence, it 
would seem that the discrepancy between his and Carpenter’s results cannot be explained 
in this manner. Whatever the cause may be, it would seem advisable, in view of the un- 
favorable distribution of Kopal’s observations in full light, not to make use of his deter- 
mination of K,. His actual measurements, however, are in very good accord with my new 
velocity-curve. His mean velocity of — 31 km/sec at phase 1.004 days falls well within the 
permissible range of my observations, even if allowance is made for the uncertainty of his 
and my phases. 

The problem of this investigation is to derive a new velocity-curve in order to ascer- 
tain whether Carpenter’s results were correct or whether Kopal’s procedure has given a 
more reliable value of Ki, and consequently of the absolute dimensions of the system. 

2. The observations.—During the long nights of November it is possible to observe 
U Cephei, with its high northern declination, for 12} hours each night. Since the period is 
almost exactly 2.5 days, it is possible to cover the entire cycle in five consecutive nights. 


12 Bol. Obs. Madrid, 1, No. 16, 1933. 


| 
| 
| 
( 
I 
~. 


U CEPHEI 225 


Whenever it is desired to make certain that secular effects do not influence the results, 
such a quick coverage is especially desirable. The observations described in this paper 
were made at the McDonald Observatory, using the quartz Cassegrain spectrograph at- 
tached to the 82-inch reflector. The linear dispersion of all plates was 40 A/mm at A 3933, 
The observations were made on nine consecutive nights—November 5—13—and cover 
parts of four consecutive cycles. The best exposure time outside of eclipse was found to 
be 10 minutes, for fair seeing. In poor seeing the exposure times were increased to 15 
minutes or even longer. During the principal eclipse the star becomes so faint and so red 
that exposures of 1 hour are required to obtain satisfactory densities. The total number 
of spectrograms measured is 115. Table 1 gives a list of the stellar wave lengths used for 
the determination of the radial velocities. 


TABLE 1 
LIST OF STAR LINES USED FOR RADIAL-VELOCITY DETERMINATION 
B8 STAR 
G2 STAR 


The phases of all exposures were computed by means of the formula given in the 1940 
edition of Schneller’s catalogue of variable stars:'* 


Principal minimum = 2424253.660 + 2.492934E. (3) 


The value of the period is identical with the one used by Kopal. The epoch is almost ex- 
actly 1049 cycles earlier than the value adopted by Kopal. Hence, the phases computed 
by me are on precisely the same system as thoseused by him. 

Because of the great uncertainty that is attached to extrapolations of linear formulae 
for the epoch of central primary eclipse, I have determined, as accurately as possible, the 
time of mid-eclipse, using the finder of the 82-inch telescope for making visual compari- 
sons of U Cephei with two neighboring stars. The mean result from estimates on Novem- 
ber 8 and 13 gives 


Mid-eclipse = 1943, Nov. 8, 9°20" U.T. = JD 2431036.889 . 


This value may be in error by +10 minutes or +0.007 day. I have neglected the small 
correction from the earth to the sun. The observed minimum occurs too early, if com- 


pared to formula (3), by 
(O—C) =—0.044 day, 


so that on the system of phases here used mid-eclipse actually took place between phases 
2.442 and 2.456 days. This departure of a little more than 1 hour is unquestionably real. 


13 Kl. Veroff. Berlin-Babelsberg, No. 21, p. 195, 1939. 
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The same departure is found when we use Kopal’s formula. It is customary to determine 
the departure (O — C) from Wendell’s elements: 


Mid-eclipse = 2407890.301 + 2.4928840E . (4) 


The result is 

(O — C) Wende = +0.160 (E= 9285). 
Miss Edith Jones'* found in 1937: O — C (Wenden) = +0.1515 for E = 8262, and 
Tecza! obtained almost the same result.’® My determination of (O — C) suggests that 
after a prolonged interval, during which the period of U Cephei increased, it has now 
started to decrease very slightly. If the observations by Miss Jones, Tecza, and myself 
are plotted on Tchudowitchev’s curve relating (O — C)wenden and E, we find that the 
last observation falls considerably below the extrapolated curve and suggests a bending 
which may be due either to an error in Tchudowitchev’s period for the apsidal motion or 
to the remaining inequalities which were considered by Carrasco but which have not, as 
yet, been accounted for physically. 

Table 2 contains the radial velocities as measured on my plates. The lines are excep- 
tionally poor, except during totality; and the settings are very-uncertain. Hence, a large 
amount of scatter is to be expected. 

3. Description of the system of U Cephei.—The photometric and spectrographic in- 
vestigations are in full accord reg: “ting the general characteristics of the system. The 
brighter star is of type B8 or B9 (not AO, as was heretofore believed, because the He 1 
lines are visible and He 14472 = Mg 11 4481) and has a slightly elliptical shape, with 
?max = 0.18a and rmin = 0.17a. The light of this star is 0.84 of the total luminosity. The 
larger star is of type gG2 (not G8 or KO, as was heretofore believed).!’ Its larger semi- 
axis iS fmax = 0.42a, and its smaller semiaxis is rmin = 0.32a. According to Tchudowit- 
chev, this star is much more elliptical than the B8 star. It contributes only 0.16 to the 
luminosity of the system (in visual light) and can be observed only during the principal 
eclipse. This eclipse is total, and the lines of the B8 star disappear completely during the 

total phase. The system has an appreciable amount of darkening at the limb (the coeffi- 
| cient of darkening u = 3 for visual light, according to Dugan) and varies with the wave 
| length, according to Rosenberg.'* It has, however, been shown by Krat'® that, although 
Rosenberg’s discussion furnishes a qualitative indication of the change of darkening 
with A, his numerical values are dependent upon the elements of the light-curve which 
he used. Dugan found an asymmetry in the light-curve, in the sense that the total light 
prior to eclipse is smaller than after eclipse and that the descent toward total phase is less 
steep than the ascent after totality. This asymmetry has been confirmed by most other 
observers, except by Rosenberg, who, instead, found a considerable decline in brightness 
during totality. Dugan suggested that the asymmetry of the light-curve can be explained 
by assuming that the B8 star rotates in a period shorter than 2.5 days and that this rapid 
rotation carries the tidal bulge forward by 24°. Tchudowitchev finds a tidal lag of 26° 
and a libration at principal minimum of 8°. He attributes this libration to the G2 com- 
ponent, which presumably rotates with a period of 2.5 days—equal to that of orbital 


revolution. 


14 Bull. Harvard, No. 905, 1937. 18 Acta Astr., C, 3, 131, 1937. 

16 T am indebted to Professor R. Prager, of the Harvard Observatory, for a list of recent references 
of photometric observations of U Cephei. 

17T am indebted to Professor W. W. Morgan for helping me with the classification of this spectrum 
and for furnishing a series of standard spectra having approximately the same dispersion as my spectro- 
grams of U Cephei. 

18 Ap. J., 83, 67, 1936. 

19 Verdnder. Sterne, Nishni Novgorod-Gorki, U.S.S.R., 5, No. 7, 214, 1938. 


PLATE X 


Phases 1.939,2.318,2.469,0.038,0.041,0.102, 


SPECTRA OF U CEPHEI BETWEEN H¢ AND H6 


The first spectrogram was taken near maximum radial velocity; the last, near minimum radia 
velocity. The second plate, phase 2.378 days, shows narrow lines of Hf and Hé unsymmetrically sup 
posed over the broad lines from the limb of the B8 star. The third plate shows the spectrum of the 
star during totality. The plates at phases 0.038 and 0.041 show the doubling of Ca m and H after totality. 
The plate at phase 0.102 shows unsymmetrical lines of H and Ca 11 after totality. 


if 


i 
H 
i 

bd t 

~ 


TABLE 2 
RADIAL VELOCITIES OF U CEPHEI 


Date Velocity Re- Date Velocity Re- 
No. 1943 GE: Phase (Km/Sec) | marks* No. 1943 Ut. Phase (Km/Sec) | marks* 

Nov. Nov. 
5 | 4>48™) 1.752 | + 80.2 2505 9 | 8>48™ 0.934 | — 62.1 
2526. ss 5 18 | 1.773 | + 81.6 2596... 9 00 | 0.942 | — 43.6 
co) 7 11 | 1.851 | + 93.3 2591.:: 9 12 | 0.950 | — 68.7 
2529... . 7 23 | 1.860 | +107.8 2599... 10 45 | 1.015 | — 57.6 
25350. «5: 10 48 | 2.002 | +148.0 2600. . 10 59 | 1.025 | — 63.2 
Co) 11 01 | 2.011 | 4154.3 2602. . 12 27 | 1.086 | — 45.2 
2535.4 6| 8 02 | 0.395 | — 92.8 2603... 12 43 | 1.097 | — 53.2 
2536... 8 15 | 0.404 | — 85.5 2604. . 10 | 42 | 1.596} + 71.0 
0. | 9 52 | 0.471 | —115.0 2005... 0 55 | 1.605 | + 69.3 
y's ee 10 07 | 0.482 | — 97.7 2606. . 1 09 | 1.615 | + 68.9 
yoo) oe 10 26 |.0.495 | —106.2 2607... 1 22 | 1.624 | + 49.3 
Ao) 10 47 | 0.509 | —104.9 2609. . 3 25 | 1.709 | + 79.8 
tr 11 08 | 0.524 | —104.0 2610... 3 38 | 1.718 | + 59.4 
pe 11 33 | 0.541 | —116.3 2682. ..: 5 03 | 1.777 | + 58.8 
2544.... 11 56 | 0.557 | — 98.0 2613... 5 16 | 1.786 | + 67.5 
TS... 12 15 | 0.570 | —105.8 2615... 6 57 | 1.857 | +108.4 
2546.... 12 33 | 0.583 | —111.8 2616. . 7 08 | 1.864 | + 90.2 
2547... . 12 46 | 0.592 | — 95.0 2618... 8 44 | 1.931 | +117.5 
2545... 7 0 40 | 1.088 | — 77.3 2619: ... 8 55 | 1.939 | +140.2 
2599. x. 0 52 | 1.096} — 53.8 262). ... 10 45 | 2.015 | +151.4 
1 04 | 1.104 | — 51.1 202%... ... 11 18 | 2.038 | +136.7 
(55 in 118 | 1.114] — 43.6 2622; .. 11 59 | 2.067 | +127.8 
2503: ...: 3 14] 1.195 | — 53.6 2623... 11 0 40 | 0.102 | — 63.8 (1) 
2558: 3 32 | 1.207); + 2.4 2624... 1 01 | 0.116 | — 69.2 (1) 
3 58 | 1.225 | — 37.8 1 20 | 0 130 | — 49.0} (1) 
i 8 0 40 | 2.088 | +130.9 2626... 1 36 | 0.141 | — 31.8 (1) 
0 53 | 2.097 | +135.0 2627 =... 2 01 | 0.158 | — 56.3} (1) 
2559... 1 22 | 2.117 | +117.0 2628: 217 | 0.169 | — 55.1} (1) 
2500. ... 1 42 | 2.131 | +113.4 2630... 4 18 | 0.253 | — 81.5 
2561... 2 09 | 2.150 | +127.6 2631... 4 32 | 0.263 | — 80.7 
2563; .. 3 44 | 2.216 | + 86.1 2632. 4 47 | 0.273 | — 83.1 
2564... 3.59 | 2.226 | + 69.9 2634... 6 26 | 0.342 | — 84.6 
4 13 | 2.236 | + 82.0 263506. 6 42 | 0.353 | — 86.2 
2566)... 4 30 | 2.248 | + 54.6 2636... 7 00 | 0.366 | — 95.2 
2567... .. 4 51 | 2.262 | + 50.8 2649... 12 5.02 | 1.284 | + 28.0 
2568... . 5 14 | 2.278 | + 33.6 2650... 5 14} 1.292 | + 40.8 
5 41 | 2.297 | + 26.3 5 1.299} + 9.6 
v7) ae 6 12 | 2.318 | +117.4 2054... 6 47 | 1.357 | + 21.9 
yi ee 6 49 | 2.344 | +150.6 2655... 6 58 | 1.364 | + 30.4 
7 32 | 2.374 | +192.6}| (1) 2056... 7 09 | 1.372 | + 22.7 
8 15 | 2.404 | — 39.0 (2) 2660. . 9 22 | 1.464} + 38.2 
2574... 9 18 | 2.448 | — 26.0] (2) || 2661... 9 33 | 1.472| + 36.4 
251. . : 10 22 | 2.492; + 5.8 (2) 2662. . 10 02 | 1.492 | + 46.5 
2516... 11 18 | 0.038 | + 21.6} (2) 2663. . 10 13 | 1.500 | + 27.9 
2577 ... 12 04 | 0.070 | —159.2| (1) 2664. . 10 25 | 1.508 | + 32.2 
Ze... 12 34 | 0.091 | —101.5| (1) 2665... 10 36 | 1.516} + 61.9 
2579. . 9! 0 54] 0.605 | —114.3 2672... 13 | 3 43 | 2.229] + 71.6 
2580... 1 06 | 0.613 | —104.3 2673... 3 54 | 2.236 | + 66.2 (1) 
2581. . 1 18 | 0.621 | —111.4 2641. 6 57 | 2.364 | +163.6} (1) 
1 32 | 0.631 | —107.9 2678... 718 | 2.378 | +195.2| (1) 
2585... 4 03 | 0.736 | — 95.1 2679... 7 45 | 2.397 | — 18.1] (2) 
2586. . 4 31 | 0.755 | — 83.5 2680. . 8 19 | 2.421} — 21.2] (2) 
4 43 | 0.764 | — 94.1 2681. . 9 29 | 2.469} — 2.4] (2) 
4 54 | 0.771 | — 59.1 2682... 10 29| 0.018 | + 8.0} (2) 
2590... 6 21 | 0.832 | — 47.2 2683... 11 02 | 0.041 | —244.4} (1) 
fs.) | 6 38 | 0.843 | — 70.2 2684... 11 35 | 0.064 | —156.9| (1) 
Yt 6 56 | 0.856 | — 78.5 2685... 12 00 | 0.081 | —135.2}| (1) 
1 8 33 | 0.923 | — 57.1 


* (1) Velocity of B8 star from Cau, Het, Mg ut, and Si lines; (2) velocity of G2 star. 
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The spectrum outside of total eclipse consists of a strong, normal set of H lines and of 
very weak and shallow lines of Ca 11 3933, He 1, Mg 1 4481, and Siu. These lines are all 
exceedingly broad; and even the H/ lines have broad centers, suggestive of rapid axial 
rotation. The lines other than those of H are exceedingly weak and are difficult to meas- 
ure.?° Carpenter did not mention them at all; and Kopal did not measure them, although 
Redman had mentioned them in his paper.”! The lines are very weak on my plates, but I 
have found them very useful for radial-velocity measurements. An estimate from the 
width of Ca 11 3933 gives for the equatorial velocity of rotation 


V (rot B) = 180 km/sec . 


This may be in error by as much as +50 km/sec. It is not possible to improve the deter- 
mination appreciably by measuring actual contours, because the total intensities of these 
lines are very low. The strong H lines are, of course, not suitable, because they are in- 
fluenced by Stark effect. It is, however, certain that the rotational velocity is very much 
larger than the value of 50 km/sec obtained by Kopal from his measurements of the 
Victoria plates. He had used this value in support of the hypothesis that there exists 
synchronism of rotation and revolution in the early-type star, whose radius he had found 
to be 1,800,000 km. A rotational velocity of 50 km/sec is ruled out by two additional 
considerations: (1) so small a velocity produces a hardly noticeable broadening of the 
lines in full light, with a dispersion of 40 A/mm (as used in my observations, at \ 3933); 
and (2) the lines of the secondary are actually sharp on my plates, yet, if synchronism ex- 
isted for both stars, they should be at least 50 per cent broader than the lines of the B8 
star. I estimate from my plates that the rotational velocity of the G2 star is 


V (rot G) < 50 km/sec . 


We have here a result of exceptional interest, which in a most unexpected manner con- 
firms Dugan’s hypothesis for the explanation of the tidal lag of the early-type compo- 
nent. The B8 star rotates three or four times faster than the G2 star, despite the fact 
that its radius is 50 per cent smaller. . 

4. The velocity-curve-—The results of my measures are plotted in Figure 1. The mini- 
mum is at phase 0.56 day and gives V = —105 km/sec. The maximum is at phase 2.00 
days and gives V = +135 km/sec. The range in velocity is therefore 240 km/sec and 


K,=120 km/sec. 


This is close to the value obtained by Carpenter (Ki; = 110 km/sec) and is almost twice 
the value derived by Kopal. Although the scatter is quite large, K, cannot be greatly in 
error. The curve is unsymmetrical, the interval from minimum velocity to maximum 
being 1.44 days and the interval from maximum to minimum being 1.05 days. This 
asymmetty is very pronounced and is, without question, real. The maximum is narrow, 
while the minimum is fairly broad. If we wish to interpret these results in terms of orbital 


elements, we should find 
e=0.20, w=+40°. 


Unfortunately, it is futile to try to improve these elements. The velocity-curve is so seri- 
ously distorted on its important descending branch by the rotation effect, and it also 
shows such a peculiar deformation at maximum, that no amount of refinement can give 


2°On many plates outside of eclipse the lines of Catt K and Het look double or unsymmetrical. 
Sometimes they are single but fairly narrow, while at other times they are very broad and diffuse. There 
appears to be no regularity in these observations, and I believe they are not physically real. Spectral 
lines which are always near the limit of visibility are greatly influenced by accidental grain patterns in 
the emulsion. A somewhat similar conclusion has been reached by Redman. 


21 96, 488, 1936. 
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elements of greater precision. The curve departs very markedly from Carpenter’s veloc- 
ity-curve. In terms of orbital elements this difference is principally illustrated by the two 
sets of values for w and e: 


Carpenter: w= +25° ; e= 0.47, 

Struve: w==+40°; e=0.20. 
Carpenter’s curve shows a much broader minimum than my curve, and I am unable to 
account for this difference. Moreover, my curve descends to V = —105 km/sec, while 


Carpenter’s minimum falls near V = —80 km/sec. The two maxima are also somewhat 
different. It is impossible to say whether there has been a real change in the velocity- 


km/sec 
+200 T T 


+ 3 


+100}— 


—150} 


-200 L | L 1 4 | 4 L 


2 elipse 3 days 


Fic. 1.—Velocity-curve of U Cephei. Solid circles represent the velocities of the B8 star. Crosses 
designate the velocities of the G2 star. The rotational disturbance during the partial phases of the 
eclipse is particularly striking. The rest of the curve is noticeably unsymmetrical, suggesting an eccen- 
tricity which is incompatible with the photometric investigations. 


curve since Carpenter’s observations were made, in 1923-1925. It is entirely possible 
that the great uncertainty of his and of my measures may account for the difference. Al- 
though the difference in e¢ is large, it is significant that the asymmetry of my velocity- 
curve does resemble that found by Carpenter; it should be noted that the inequality of 
the “crest-intervals,””? which is the most conspicuous feature of my curve, depends upon 
w as well as upon e. The same inequality could be obtained with e = 0.3, w = 60° and 
with e = 0.40, w = 70°. The reason for preferring e = 0.20, w = 49° is the almost straight 
ascending branch of my velocity-curve. I believe we can be certain that the shape of the 


2 See E. S. King, Harvard Ann., 81, 231, 1920. 
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velocity-curve is not in harmony with the photometric results. The value of y which best 


fits my observations is 
y=—5 km/sec, 


in excellent agreement with Carpenter. The measures of the G2 star at totality give the 
slope of the secondary’s velocity-curve between phases 2.37 and 0.05 days. Although the 
scatter of my points is quite large—partly because of underexposure of several plates and 
partly because of the greater slit width used for these exposures—the slope is determined 
with a fair degree of precision and is almost identical with the slope of the primary curve 
undisturbed by rotation. The mass ratio is thus close to 


If we assume that the shape of the velocity-curve of the G star is similar to that of the 
B star, its semiamplitude would be 


K,= 120 km/sec. 


If, on the other hand, we assume that the velocity-curve of the B star is distorted by un- 
known causes, while that of the G star is genuine, we can compute Ke by means of 


Kopal’s formula: 


2 
K:tt+y. 


We find from my velocity-curve that for 
Phase 2.60 days: V2= +50 km/sec , 


Phase 2.30 days: V2= —60 km/sec . 


Hence, 
P 


Ky 


This is much less than the value found by Kopal. I am unable to account for the fact 
that the slope of my observations is different from his. Possibly the best that can be done 
at the present time is to use the mean slope of Kopal’s measures and mine, which would 
be consistent with Ke = 173 km/sec and mg/mg = 0.7 or 0.8. 

At the time of partial eclipse the velocity-curve is distorted by the rotation effect. 
This phenomenon is discussed in the following section. The best possible elements from 
my velocity-curve are 


P= 2.493 days , w= 40°, 
y= —5 km/sec, T = phase 2.16 days , 
e=0.20, a,sini=4X10%km, 
i _ 
K,=120 km/sec , 0.420. 


The dotted line which in Figure 1 represents the radial velocities of the G2 component 
intersects the interpolated velocity-curve of the B8 star at y = —5 km/sec. This is a 
confirmation of our elements. Although my value of y agrees closely with that deter- 
mined by Carpenter, the minima of the two curves differ appreciably, and the agreement 
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in y may be accidental. Kopal*’ has remarked that in U Cephei “‘the presence of a third 
body seems to be ruled out by the constancy of the radial velocity of the center of mass,” 
but it is not clear upon what evidence this statement is based. Kopal’s own value of y 
differs considerably from Carpenter’s and from mine. There is probably not enough in- 
formation available at the present time to ascertain whether y remains constant. 

If the spectrographic elements were consistent with the photometric observations, 
mid-eclipse should take place when v + w = 90° or 270°, or when the radial velocity is 


V.=7y+Ke cosw=+13.5 km/sec. 


This predicted mid-eclipse would fall at phase 2.400. The observed value was 2.449. 
The difference of O — C = +0.049 is too large to be caused by errors in the photometric 
observations. 

5. Spectroscopic phenomena during partial eclipse-—The durations of the entire eclipse 
and of its total phase are 10"5"6 = 0.421 day and 1554™8 = 0.080 day, respectively.”4 
Since on the system of phases used in this paper mid-eclipse occurs at phase 2.449 = 
—().044 day, the phases of the four contacts should be 2.239, 2.409, 2.489, and 0.166 days. 
The uncertainty of each of these phases is estimated to be +0.007 day. If we follow the 
velocity-curve forward from the maximum, we find good agreement until phase 2.300. 
Thereafter the velocity-eurve rises very suddenly to a maximum of +195 km/sec at 
phase 2.38 and then drops to a minimum at —244 km/sec at phase 2.53 = 0.04 day. 
This is the rotation effect of the B8 star. At the two critical phases we observe the veloc- 
ities of the receding and of the approaching limbs of the B8 star, and we may therefore 
take for the rotational velocity at the equator (since i ~ 90°) 


195+244 
2 


: V (rot) = = 220 km/sec. 


Our direct estimate from the widths of the undisturbed lines was 180 km/sec. The agree- 
ment is satisfactory and confirms the unusual absence of synchronization of the periods of 
rotation and revolution. The rotational velocity of the B8 star is twice its orbital veloc- 
ity, and the rotational effect observed during partial eclipse is the largest found for any 
eclipsing star. The lines used for the determination of the radial velocities in the partial 
phases include all lines of He 1, Mg u, Ca 01, and Si 1 but do not include the lines of H. 
The reason for omitting H will be explained at once. 

Table 3 lists the details of my measures during the primary eclipse. The material is 
arranged in such a way as to show at one glance the division into three groups: first par- 
tial phase, totality, and second partial phase. The velocities of the H lines have been seg- 
regated from those of the faint lines of the B8 star: He1, Mg u, Cau, and Si 1. It will 
be seen that the H lines show a smaller rise than do the lines of He 1, Mg 11, and Ca 11 at 
phase 2.38 days. There may also be a somewhat smaller drop of the H lines at phase 
2.53 = 0.04 day, although the departure is not as striking as it was before totality. The 
question arises whether the H lines or the group of He 1, Mg u, Ca 11, and Si 1 show the 
true rotational effect of the B8 star. 

In order to answer this question, we turn to a description of the B8 spectrum during 
the partial phases. Nothing abnormal can be seen until phase 2.278 days. At phase 2.297 
the H lines are distinctly unsymmetrical, with a strong core on the violet side of the dif- 
fuse, Stark-broadened line. At phase 2.318 these cores have become somewhat stronger, 
and they are very sharp. These fine lines are superposed over the fairly broad and un- 
symmetrical feature produced by the partly occulted rotating disk of the B8 star. They 


23 Ann. New York Acad. Sci., 41, 46, 1941. 
*4 Dugan, Contr. Princeton U. Obs., No. 5, p. 34, 1920. 26 Formula (3). 
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do not look like the unsymmetrical minima of the rotational contours. The sharp cores 
are still more conspicuous at phases 2.344, 2.364, 2.374, and 2.378 days, although they 
always remain narrow and sharp. It is clear that I have measured these cores, and not the 
underlying broad lines, which at all these phases are displaced toward the red with respect 
to the cores. The lines of He 1, Mg 11, and Ca 11 are too weak to show the narrow cores, 
even if they existed. Hence, these lines show the undisturbed rotational shift. 


TABLE 3 
RADIAL VELOCITIES OF U CEPHEI DURING ECLIPSE 
RapiAL VELOCITY 
DATE Exp. RE- 
No. 1943 U.T. H MARKS 
ot, Com, H G Star 
Mgu, Siu 
8 4 30 2.248 20 + 54.6 
8 4 51 2.262 22 + 50.8 
8 5 14 2.278 22 + 33.6 
8 5 41 2.297 30 + 26.3 
8 6 49 2.344 43 +150.6 
13 6 57 2.364 20 +163 .6 
ee 8 7 32 2.374 35 +192.6 + 90.4 —46.0 
_ 13 7 18 2.378 20 +195.2 + 37.2 —16.2 (1) 
8 10 22 2.492 + 5.8 
: ae 13 11 02 0.041 30 —244.4 —194.4 +25.4 (2) 
8 12 04 0.070 30 —159.2 (4) 
8 12 34 0.091 28 —101.5 
11 0 40 0.102 21 — 63.8 
11 1 O1 0.116 20 — 69.2 
11 1 20 0.130 — 49.0 
11 2 01 0.158 15 — 56.3 
11 217 0.169 15 — 55.1 


i) haunt component of Ca 1 K gives +49.0; (2) of Ca 11, H6 and Hy, gives —53.3; (3) of Hé gives —42.8; (4) of CaurK 
gives —43.8, 


In order to be quite certain that the sharp H cores represent superposed lines and are 
not to be identified with the unsymmetrical minima of the rotational contours from the 
partly eclipsed disk of the B8 star, it is necessary to investigate these contours. Thus far, 
the only investigation of such contour (as distinct from studies of the shift of the center 
of gravity of unsymmetrical lines) was made by Struve and Elvey”* in the case of Algol. 
The eclipse of this star is partial; and the resulting contours become very unsymmetrical, 
with the steeper edge on that side toward which the line, as a whole, is shifted. The ob- 


26 91, 667, 1931. 
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servations very clearly demonstrated that the lines were really unsymmetrical in the 
sense predicted. But it should be realized that, while in a binary system in which the 
direction of axial rotation of each component coincides with the direction of orbital mo- 
tion the shift of the line is always in the positive sense before totality and in the negative 
sense after totality, there is no such uniformity with regard to the asymmetry of the 
line. Under certain conditions the red wing may be the steeper, while under other condi- 
tions the violet wing may be the steeper. 

We shall assume that the B8 star has a circular disk during the eclipse and that its 
equatorial velocity is vo. Then 


If we consider the center of the star’s disk as coinciding with the center of co-ordinates, 
the x-axis falling along the equator, and the y-axis along the axis of rotation, then, if the 
radius of the B8 star is 1, 


v= Vox 

and 
Ado" 


The contour of an infinitely narrow line will be converted by rotation, in the absence of 
an eclipse and without limb darkening, to 


If the intrinsic line contour is given by an expression J(«), then the actual broadened 
contour is 


+o 
I'(x) = fT (y) (x-y) dy. (6) 


Equation (5) suggests a very simple graphical procedure for drawing the contours of lines 


distorted by the eclipse. If the radius of the B8 star is taken equal to 1, then V1 — x rep- 
resents the equation of a half-circle whose diameter is given by x = +1. For any value 
of «, the square root is the half-length of the chord drawn on the disk of the B8 star at 
right angles to the abscissa. If there is an eclipse, some of the chords are cut short by the 
disk of the eclipsing star. If we draw to scale the relative positions of the disks of the 
eclipsing and eclipsed stars and measure the uneclipsed sections of vertical chords on the 
B8 disk, we obtain at once a representation of the corresponding ‘‘pure”’ rotational con- 
tour. The actual contour can then be easily obtained by carrying out integration (6). 
In Figures 2 and 3 I have shown several typical cases. Figure 2, a, corresponds to the 
system of U Cephei. The inclination is 90°, and the radius of the eclipsing star is 1.8 
times that of the bright star. At first contact the contour /(x) is undisturbed. At second 
contact, /(x) = 0, since totality sets in at that moment. The intermediate contours are 
for }, 3, and 3 of the partial phase, in time. The five contours correspond to phases 2.239, 
2.281, 2.324, 2.367, and 2.409. The theoretical contour for “edge 3,” or phase 2.367, cor- 
responds to the time when we observe the narrow H core, giving a velocity of the order of 
+100 km/sec, as compared to a velocity of nearly +200 km/sec for the lines of He 1, 
Mg u, and Cau. It is quite certain that the theoretical contour does not give rise to a 
sharp minimum at one-half of the displacement of the center of gravity of the line. There 
can be no doubt that broadening by Stark effect will produce for the case “edge ?”’ an al- 
most symmetrical, although greatly shifted, contour. Incidentally, the contours at 
“edge }” and “edge 4” are interesting because of the sharp cutoff on the violet side which 
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they predict. We do not observe this cutoff, which is another argument in favor of the 
suggestion that Stark effect obliterates the details in the H lines, while the other lines are 
too weak to show the contours clearly. Figure 2, b and c, and Figure 3, a, b, and c, corre- 
spond to other conditions and show the variety of contours that may be observed in 
eclipsing variables. The contours, as drawn, have not been normalized to retain their 
areas constant. To do this, we should have to take 


I" (x) = (x) 


But it is not necessary to carry out this normalization if it is only desired to compare the 
shapes of the resulting curves. In Figure 2, 6, the diameters of the two stars are the same. 


Red Violet Red Violet Red Violet 


ist 


2nd contact 


Fic. 2.—Typical contours of intrinsically narrow spectral lines of rapidly rotating stars in different 
stages of eclipse. The first case corresponds to U Cephei, where the eclipsing star is 1.8 times larger than 
the bright eclipsed star. In the second case the diameters of the two stars are the same, and in the third 
case the eclipsing star is the smaller. In all three cases the inclination is 90°. The contours have not been 


normalized to preserve constancy of the equivalent widths. 


In Figure 2, c, the eclipsing star is the smaller. In Figure 3, a, 6, and c, the inclination dif- 
fers from 90°. 

I believe it is safe to conclude that, while the asymmetry of the contours in U Cephei is 
such as to make the violet edge of the line slightly steeper than the red, the observed 
sharp core cannot be due to this asymmetry and must be produced in a stream of gas 
whose motion is in the same direction as the direction of rotation of the B8 star but whose 
velocity is only about one-half of the equatorial velocity of this star. The sharpness of the 


U CEPHEI 235 


superposed absorption lines shows that the gas is at very low pressure, because there are 
no wings produced by Stark effect. We may tentatively attribute these lines to a stream 
which starts from the B8 star and passes along the following side of the G2 star, to return 
toward the B8 star along the preceding side of the G2 star. This interpretation would be 
consistent with that advanced in the case of SX Cassiopeiae.?’ 

The lines of the G2 star are not at all visible at phases 2.318 and 2.344. They are very 
weakly, but unmistakably, present at phase 2.364 and later. But the radial velocity of 
the G star at these phases is of the order of — 40 km/sec. The H cores give +100 km/sec. 
Hence, they are not to be identified with the H lines of the fainter component. Further- 


Fic. 3.—Contours of intrinsically narrow spectral lines of rapidly rotating stars in different stages of 
eclipse. The ratios of the diameters of the two stars are the same as in Fig. 2, but the inclination is not 90°. 


more, the sharp H cores are present at H¢ and at other extreme violet members of the 
Balmer series, which are not strongly enough exposed to show the G star at all. At phase 
2.378 there is a clear indication that Ca 1 K is double. But it is impossible to separate the 
effect of the G star completely. The exposure lasted 20 minutes—from 2.371 until 2.385. 
The continuous spectrum of the B8 star is quite strong, as is shown by the broad ultra- 
violet Balmer lines. The lines of the G star are quite weak, and it is almost certain that a 
20-minute exposure would fail to register its continuous spectrum at \ 3933. The first 
G-type line on the violet side is Fe 1 4005, which appears very weakly. I am almost cer- 


27 Struve, op. cit., p. 94. 
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tain that the duplicity of Ca 1 K has nothing to do with the G-type spectrum. Either it 
may be due to an effective separation of a Ca 1 line in the stream from one produced in 
the rapidly rotating limb of the B8 star, or it may be due to an emission line in a “‘flash- 
spectrum” just preceding total eclipse. This feature is not brought out in the notes to the 
radial velocities of Table 3. 

In connection with my observation of the ultraviolet lines of the B8 star at phase 
2.378, it is of interest to recall Redman’s** measures of the equivalent widths of HB, Hy, 
and #6 in the partial phases of U Cephei. He found that these lines become gradually 
weaker as totality is approached and that they may disappear completely at the extreme 
limb of the B8 star. Redman*® has remarked that this behavior is to be expected, not 
only in the case of lines produced by “‘pure absorption,” but also in the case of lines 
broadened by Stark effect. The evidence from my plates seems to confirm the weakening 
of the H lines toward the limb, accompanied by an appreciable reduction in width. But 
at all phases at which the B8 star was not totally eclipsed, its Balmer lines are clearly 
visible. Even at phase 2.397 (from 2.386 until 2.408) the higher Balmer lines of the B8 
star are present. If it is recalled that second contact came at 2.409, the curve given by 
Redman for the equivalent widths predicts an intensity of the H lines which is roughly 
one-quarter of that observed outside of eclipse. Although my plates are not suitable for 
accurate spectrophotometry, I should have estimated the ultraviolet H lines to be con- 
siderably more intense—perhaps of the order of 0.5 their normal equivalent widths. 
These ultraviolet lines are completely free of the disturbing effect of the G star, for which 
Redman made suitable corrections in the case of HB, Hy, and H6. 

During the totality, only the G-type spectrum is seen. On one plate, phase 2.448, I 
suspect the presence of central emission lines at Ca 11 K and H. There is no change in the 
spectrum during the total phase. At phase 0.018 the specirum is still entirely that of the 
G star, in spite of the fact that third contact came at 2.489 = —0.004 day. Perhaps my 
principal minimum is slightly in error. But, it is also possible that the asymmetry in the 
light-curve observed by Rosenberg near third contact reduces the brightness of the B8 
star to such an extent that it does not show until phase 0.038. At this phase the lines of 
the G star are still weakly present. In fact, they are seen as late as phase 0.041. The ultra- 
violet Balmer lines of the B8 star are strong at phase 0.038; but they look much narrower, 
and deeper, than outside of eclipse. At phases 0.038, 0.041, 0.064, and perhaps 0.070 the 
lines Ca 11 K, H6, and Hy are conspicuously double; probably HB, He, and H¢ are also 
double on one or two plates. This doubling is much more conspicuous than that observed 
on one plate before totality. The violet components are in all lines fainter and more dif- 
fuse than the red components. This is especially well shown in the case of Ca 1 3933. The 
violet components of Ca 11 and H agree with He 1 and Mg 11 and undoubtedly represent 
the narrowed unsymmetrical contour of the approaching limb of the B8 star. Their veloc- 
ity is of the order of —150 to —200 km/sec. The red absorption components of H and 
Cait K give about —45 km/sec and must be caused by the approaching side of the 
stream of gas whose presence we already detected in the pretetality observations. Again 
these lines cannot be identified with the H/ lines of the G star because the velocity of the 
latter is now of the order of +20 km/sec. They are much stronger than before totality 
and more diffuse, but the H lines have the appearance of lines broadened by turbulence 
and not by Stark effect. It is interesting that the stream on this side shows a fairly strong 
line of Ca 11 K, while on the other side we only suspected a trace on one plate. 

I am not yet entirely certain whether the narrow space between the two absorption 
components of Ca u K, Hé and Hy, is continuous spectrum or true atomic emission. The 
rotational contour becomes quite narrow for phases which are close to totality, but I am 
a little surprised to find it sufficiently narrow to give such a clear separation of the ab- 
sorption lines. On the other hand, Ca u K shows the doubling more distinctly than H6 


28 Op. cit., p. 504. 29 Observatory, 59, 287, 1936. 
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or Hy, and this is what should be expected from the absence of Stark effect in Cam. At 
phase 0.081 the two absorption components are blended, producing a single unsymmetri- 
cal line, with a strong and rather broad core on the red side. This core corresponds to the 
narrow sharp core observed on the violet side of the unsymmetrical line before totality. 
This appearance of the H lines continues, with gradually diminishing asymmetry, until 
phase 0.141, which is well past the } stage of the partial eclipse. It is about at this same 
phase that the measured radial velocities, which are now indistinguishable for the B8 
star and the stream, again agree with the velocity-curve computed without regard for the 
rotation effect. 

The question arises whether the stream of gas which produces the extra “satellite” 
lines in U Cephei (to use an expression coined by Baxandall to describe a similar phenom- 
enon in 6 Lyrae) may not, in reality, be the rotating atmosphere of the G2 star. There is 
a distinct possibility that we observe the light of the B8 star shining through the outer 
atmosphere of the G2 star. The velocity, +108 km/sec before totality and — 194 km/sec 
after totality, is, however, much too large to be consistent with this hypothesis. We had 
found that the equatorial velocity of rotation of the G2 star, V., is less than 50 km/sec. 
Yet, the light of the B8 star would give us an integrated line for regions of the G2 star 
from latitude 0 to some higher latitude, where the rotational velocity is less than V.. 
Although the size of the G2 star is such that this effect is small, it is difficult to accept a 
value for V, greater than 50 km/sec. Another difficulty is the fact that only H and Ca 
show the satellite components. In the normal spectrum of the G2 star Ca 11 is enormously 
stronger than H, and Fei and Ca1 are quite similar to H. The stream would therefore 
have a different stage of ionization and excitation than the normal atmosphere of the 
G2 star. Finally, geometrical conditions are such that a hypothetical ¢ Aurigae effect 
would be greatest at or before the halfway point of each partial eclipse. In reality the 
satellite lines are conspicuous between the three-quarter point and totality. I am in- 
clined, therefore, to favor the stream hypothesis, although the three arguments ad- 
vanced do not entirely rule out the ¢ Aurigae effect for the extra components observed at 
phases 2.378, 0.041, 0.064, and 0.070 days. 

Redman has pointed out that it would be of great interest to follow lines other than 
those of H during the partial eclipse of the B8 star. This is especially true because we have 
shown that the H lines themselves are not simple lines, originating in the atmosphere of 
the B8 star, but are complicated by absorption lines produced in a detached stream of 
gas. The lines of He1 are probably least affected by disturbing factors. If we consider 
first the second half of the eclipse and follow their appearance from phase 0.169 back- 
ward, we observe that they, together with Ca 1 K, Mg 11 4481, and probably the lines of 
Si 1, become narrower and deeper—precisely as we should expect from the rotational 
effect. At phases 0.116, 0.102, and 0.091 all these lines are fairly conspicuous, although 
they are still less intense than we would expect to find them in a nonrotating B8 star. 
These phases are close to the halfway mark of the partial eclipse (phase 0.081 day). From 
there on, and especially at phases 0.070, 0.064, 0.041, and 0.038, the He 1 and Mg 11 lines 
not only cease to increase in depth but become appreciably weaker. At phase 0.038 they 
are present and can be measured, but they are much less conspicuous than at phase 
0.116. The Mg m1 line behaves like He 1, but I am not entirely certain that this effect is 
free of influence from the continuous spectrum of the G star. There is no such difficulty in 
the case of He 1. With the short exposure times during partial eclipse the line \ 4026 is 
quite free of any continuous spectrum of the G star. The Ca 1 K line continues to in- 
crease in intensity as we approach totality: the violet component at phases 0.038 and 
().041 is much more conspicuous than the undisturbed line outside of eclipse: it is narrow- 
er and very much deeper, approaching the kind of line we should expect to observe in a 
slowly rotating star. The same is certainly true of Si 11 3862, which is quite conspicuous 
at these two phases, while it is exceedingly difficult to see outside of eclipse. I believe that 
the equivalent widths of Si 11 and Ca 11 remain essentially constant throughout the par- 
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tial eclipse, while He 1 remains constant until about phase 0.116 and then decreases to- 
ward totality. 

The difference in the behavior of He 1 (and perhaps Mg 11), on one side, and of Ca 1 
and Si 11, on the other, is interesting. Those He 1 lines which could be used in these com- 
parisons belong to the series 2?P — n*D and 2'P — n'D. Both are strongly influenced by 
Stark effect. Hence, it is probable that we must regard the weakening of He 1 toward the 
limb (and, by analogy, the phenomenon observed by Redman for #7) as in the nature of 
an absolute-magnitude effect. The lines appear weaker at the limb because the pressures 
are lower, and not because the mechanism of line formation is predominantly one of 

“pure absorption.” 

In the first half of the eclipse the lines of Het are systematically weaker than in the 
second half. There is a marked increase in the conspicuousness of He 4471 from phase 
2.236 to phase 2.364, the line becoming narrower and deeper. The same may be true of 
He 1 4026 and one or two other He 1 lines. But at this phase all He 1 lines are weaker than 
at the corresponding phase after totality. Certainly, the He rlines do not continue to be- 
come more noticeable afier phase 2.364. The line Ca 11 K does, however, become more 
noticeable as totality is approached. At phase 2.374 it is fairly conspicuous, though the 
Het lines are barely visible. But this line, too, is systematically weaker at phases 2.374 
and 2.378 than at phases 0.038 and 0.041. There appears to be a real difference in the in- 
tensities of lines produced by the receding and approaching limbs of the B8 star, the 
former having the weaker lines. We have already seen that the H and Ca 11 components 
which we have attributed to the stream are also stronger after totality than before. Both 
phenomena are related, but it is impossible to determine their causes. The well-known 
asymmetry of the light-curve in the partial phases may be related to this spectroscopic 
difference in spectral line intensities. According to Dugan, the tidal bulge is carried for- 
ward by the rapid rotation of the B8 star. This bulge makes the star slightly fainter just 
before eclipse than after its end. We must, then, conclude that the lines originating from 
the tidal bulge are fainter than those which come from other parts of the B8 star’s sur- 
face. 
6. Discussion.—The velocity-curve contradicts the photometric observations. The 
former gives e = 0.20 and w = 40°. The latter, according to Tchudowitchev, are con- 
sistent with e = 0.08 and w = 99°. No amount of refinement in the solution can remove 
this fundamental difficulty: in both Carpenter’s and my velocity-curves the rise from 
minimum to maximum takes approximately 1.4 days, while the descent from maximum 
to minimum requires only 1.0 day. The observations do not reveal a definite clue to the 
solution of this problem, as they do in the case of SX Cassiopeiae. There is no dilution 
effect in the absorption lines, outside of eclipse, and hence no direct proof of the existence 
of an extended shell. But the existence at partial phases of superposed lines, originating 
in a stream of gas, links the problem with that of SX Cassiopeiae and suggests that the 
asymmetry of the velocity-curve may be produced by the same stream. Kopal remarked 
in his paper ‘that the profile of Hy (as well as of other members of the Balmer series) is 
distinctly asymmetrical during full light and in the early partial stages, when the rota- 
tional effect is still insignificant.’’ My plates show no asymmetry at phases near 1.0 day, 
where Kopal’s plates, in full light, were taken. I have suspected some asymmetry in the 
contours of the hydrogen lines near maximum radial velocity (phase 2.0 days) in the sense 
that the deepest part of the line is displaced toward the red from the geometrical center 
of the line. But I am not certain that this effect is real. 

The absence of strong emission lines may not affect the stream hypothesis. But, in the 
case of SX Cassiopeiae, the eclipse of the violet and red emission components clinched the 
argument. We conclude that, while there is a possibility that the velocity-curve of U 
Cephei is distorted by absorption i in gaseous streams, there is no such compelling evi- 
dence as in the case of SX Cassiopeiae. 
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ABSTRACT 


Attention is called to the fact that several stars are now known to exhibit appreciable spectroscopic 
eccentricities while their photometric orbits are practically circular. It is pointed out that previous at- 
tempts to explain spurious spectroscopic eccentricities have not been successful; and an attempt is made 
to attribute the spectroscopic eccentricities to distortion of the velocity-curve in a mass of gas in the 
outer parts of the stars whose density is so low that the distribution of mass in the system remains sensibly 
unaffected. The asymmetry of the light-curve of U Cephei may possibly be explained in terms of a 
theory proposed by Cowling. 


Recent researches by O. Struve! have brought to the focus of attention striking dis- 
crepancies between spectroscopic and photometric eccentricities (as well as spectro- 
scopically and photometrically determined times of conjunctions) of certain eclipsing 


ELEMENTS 
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Spectroscopic Photometric 
U Cephei.. e=0.47+0 02 (p.e.) e=0.0 
w= 25° + 3° (p.e.) Many investigators 
Carpenter, Ap.J., 72, 205, 1930 
w=40° 
Struve, Ap.J., 99, 228, 1944 
MR Cygni........ e=0.12 e<0.01 
w= 103° Swope, Harvard Bull., No. 
Pearce, J.R.A.S. Canada, 29, 914, 1940 
413, 1935 
u Herculis........| e=0.05+0,01 (p.e.) e<0.01 
w=93° + 10° (p.e.) Martin, B.A.N., 8, 265, 1938 
Sterne, Proc. Nat. Acad., 27, 
175, 1941 


variables. In systems like SX Cassiopeiae this may perhaps be explained by a special 
mechanism suggested by Struve on page 96 of the first paper cited; but Struve himself 


. pointed out that such a mechanism seems of little avail in a more “‘normal’’ system like 


U Cephei, showing few if any spectral peculiarities characteristic of SX Cassiopeiae. The 
aim of this note is to point out that such discrepancies between photometric and spec- 
troscopic elements may quite possibly be a common, rather than an exceptional, feature 
of close binary systems and call for an explanation on a more general basis. 

The fact that the orbital eccentricities, as determined spectroscopically, usually come 
out too large in comparison with the photometric eccentricities has from time to time 
been commented upon by several writers;? and the accompanying tabulation of a few such 


_ 1 Ap. J., 99, 89 and 222, 1944. The writer is greatly indebted to Professor Struve for having enabled 
him to see the latter paper in advance of publication. 


2 Cf., e.g., Luyten, Minnesota Pub., Vol. 2, No. 2, 1935. 
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representative systems shows that these discrepancies may be of varying order of 
magnitude. 

The photometric evidence that the orbits of all these systems are practically circular 
is so sound and definite that there is but little room for doubt that, in the measured 
Doppler displacements of lines of these stars, shifts of undisclosed origin must be super- 
posed upon purely orbital velocities to yield a resultant which may sometimes bear only 
a distant relation to the actual orbit. Hopes were entertained for some time that dis- 
crepancies as small as that found for u Herculis could perhaps be explained by effects, 
upon radial velocity, of limb- and gravity-darkening of distorted rotating components,’ 
or by the reflection effect‘—both of which were found, indeed, to simulate an orbital ec- 
centricity with the apsidal line parallel to the line of sight. Upon closer scrutiny these 
effects proved, however, to be too small to account for the discrepancy between spectro- 
scopic and photometric eccentricity for u Herculis° and still less so for the other stars 
listed above, not to mention the fact that the apsidal line of the spectroscopic orbit of 
U Cephei is manifestly not parallel to the line of sight. The discrepancy, therefore, re- 
mains unexplained, and we may ask with Struve, “‘. . . . if we cannot trust the velocity- 
curve .... for which the photometric observations provide the crucial test, how can we 
trust the velocity-curve of any other spectroscopic binary to give us an orbit which has 
a dynamical meaning?” This is a far-reaching question which deserves close attention. 
Since the theory has thus far been unable to provide much guidance to its approach, we 
shall attempt to survey the problem in the light of available statistical data. 

It is to be stressed at the outset that eccentric eclipsing binaries, for which photo- 
metric observations may provide the crucial test, are, unfortunately, rare. The displace- 
ment of minima specifies, as a rule, only one component of the eccentricity, namely, 
é cos w;° while é sin w can be obtained only from the unequal durations of both minima 
with usually much less accuracy. In systems with stationary or slow-moving apsides 
a separation of e and w is therefore an inherently difficult task; in consequence of which 
for only very few variables are genuine photometric determinations of e,and w available. 
From these we may conclude that, in all systems for which spectroscopic observations 
have indicated a circular orbit, photometric observations are consistent with it. This 
seems to be true for very close binaries (8 Lyrae) as well as fairly well separated ones 
(Algol); no exception from this rule is known thus far. On the other hand, spectroscopi- 
cally eccentric orbits may or may not be also photometrically eccentric. For Y Cygni, 
for instance, the spectroscopic and photometric values of e and w are in almost perfect 
agreement.’ For the majority of systems a consistency of spectroscopic and photometric 
elements cannot be readily tested because the photometric observations yield only 
e cos w;® while for certain others (such as those listed above) the spectroscopic eccen- 
tricities are definitely at variance with the photometric evidence; and the significant fact 
is that such eccentricities seem to be wholly, not partly, spurious. 

Of the three systems listed above, U Cephei appears to be not only the most striking 
representative of its group but is also observationally best known. It is, therefore, tempt- 
ing to try to relate the occurrence of anomalous spectroscopic eccentricity with certain 
other peculiarities shown by U Cephei, namely, the asymmetry of its light-curve and 


3 Sterne, Proc. Nat. Acad., 27, 168, 1941. 

4 Kopal, Proc. Amer. Phil. Soc., 86, 351, 1943. 

5 Tbid., pp. 361-362. 

®In more precise terms, the observed displacement of minima due to central eclipses specifies 
cos w(1 — é)-3. 

7 Dugan, Princeton Contr., No. 12, 1931. 

8 A combination of spectroscopic e and photometric e cos w need, therefore, not necessarily lead to a 
true longitude of periastron—a point to be borne in mind in any future discussion of apsidal motions in 
close binary systems. 
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suspected fluctuations of orbital period. It has been known now for more than thirty 
years that the light-curve of U Cephei exhibits asymmetry throughout the whole cycle. 
The principal asymmetry*generating term varies between minima as sin 26 (@ being the 
phase-angle), and its coefficient increases with effective wave length of observation 
exactly as the fractional luminosity of the secondary component.® This proves that the 
asymmetry of the light-curve between minima is due to some peculiar action of the 
secondary, not the primary, star. After the light-curve has been rectified for the variation 
of light between eclipses, the primary minimum still remains asymmetric; nearly all ob- 
servers have found that the brightness of U Cephei within minima diminishes less rapidly 
than it increases '° in visual, as well as photographic, light. 

Various suggestions have been offered to account for this asymmetry. Dugan"! 
thought it to be a consequence of “‘tidal lag” invoked by fast rotation of the primary 
component, while Walter® attributed it to bodily librations of the secondary star. 
Taylor” invoked certain resonance phenomena for its explanation, while Cowling! sug- 
gested that it may be due to a lag between variation in light and radius of a distorted 
configuration exhibiting forced nonradial oscillations. Of these the hypotheses by Dugan 
and Walter, although interesting innovations at the time of their publication, can no 
longer be maintained. It has been recognized for some time that tides cannot lag by any 
appreciable amount in systems consisting of fluid components—and there is very little 
room for doubt that the stars behave almost like perfect fluids; neither can there be any 
libration. The secondary component of U Cephei could librate in the way envisaged by 
Walter only if it were solid from the core to the surface; not otherwise. If synchronism 
exists between rotation and revolution in systems consisting of fluid components, the 
tidal bulge will form a standing wave (or a system of waves of the same phase) on the 
surface of the distorted star. If the latter rotates faster (or slower) than its mate re- 
volves, the tidal wave will merely sweep along its equator so as to remain always pointed 
to the center of mass of the disturbing body. The viscosity of stellar matter is too small 
to allow the tides to lag by any observable amount." 

On the other hand, the theoretical basis of Taylor’s explanation of asymmetries of 
eclipsing binary systems was largely invalidated by Cowling, whose own hypothesis may 
be summarized as follows: It is known from observations of Cepheid variables that the 
variation in light emitted by a radially oscillating star lags behind the variation in radius. 
Now, if the rotational and orbital velocities differ, the tidal distortion of the components 
in eccentric binary systems becomes a forced nonradial oscillation. We may then expect 
that the variation in light, which may accompany such nonradial oscillations, will also 
lag behind the vertical displacement; and this may produce the observed asymmetry. 
To this one might remark that it is not yet certain whether nonradial oscillations of com- 
ponents in close eccentric systems are accompanied by any light-variation. Such varia- 
tions in radially oscillating stars originate presumably in the deep interiors of the stars, 
while nonradial oscillations of tidal origin are a relatively superficial phenomenon. At 


Walter, A.V., 250, 342, 1933. 


10 A suggestion by Tchudowitchev (Engelhardt Obs. Bull., No. 17, 1939) that the asymmetry of minima 
oscillates in sign remains still to be confirmed. His geometrical explanation of this phenomenon as a con- 
sequence of the revolution of apsides is, however, certain to fail. Apart from the fact that any evidence 
for apsidal motion in U Cephei is still slight, an orbital eccentricity considered by Tchudowitchev is en- 
tirely insufficient to produce any observable difference in duration of the ascending and descending 
branches of the primary minimum. 


1 Princeton Contr., No. 5, 1920. 
2 Ap. J., 94, 46, 1941. 
18 M.N., 101, 367, 1942. 


4 Struve’s conclusion that the observed fast rotation of the B8 component . . . . “in a most unexpect- 
ed manner confirms Dugan’s hypothesis of the tidal lag of the early-type component” seems, therefore, 
to have been premature. The observed phenomena admit also of alternative explanations. 
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any rate, Cowling’s suggestion still remains to be tested by detailed comparison with 
the observations, and the origin of asymmetry of the light-curve of U Cephei and other 
eclipsing variables appears still uncertain. 

We may pursue the problem one step further and inquire as to whether or not the spu- 
rious spectroscopic eccentricity or the asymmetry of light-changes of U Cephei are more 
likely to be due to a mere peculiar distribution of /ight on the surfaces of its components 
or to a peculiar distribution of mass in the system itself; for, unlike the former alterna- 
tive, the latter should give rise to dynamical perturbations of simple orbital motion, 
some of which—namely, those of orbital period—could possibly be detected. It has, in- 
deed, been suspected for a long time that the orbital period of U Cephei is not entirely 
constant; at least in recent decades it has been appreciably longer than at the time of its 
discovery seventy years ago. On the other hand, the existence of complicated periodic 
variations, postulated from time to time by several writers," should still be considered as 
doubtful. The evidence for such variations was based almost exclusively on old series of 
observations from 1880 to 1910; and several considerations indicate that, with a few ex- 
ceptions," these should be given much less weight than has been done hitherto. The light- 
curve of U Cephei is now known to be asymmetric; the light within minima diminishes 
less rapidly than it increases; and even during totality it fails to remain constant. The 
bearing of these facts upon the determination of the times of minima is obvious; yet they 
were only gradually realized after decades of observations. Whenever old observations of 
the minima were published in extenso, we may reconstruct their form and often find the 
moments of central eclipse to differ appreciably from those given by the observer; but 
minima based on unpublished data (especially if the observations were few or on one 
branch only) should not be included in any precise analysis. With the advent of modern 
photometric methods the accuracy of known minima has increased considerably. The 
writer rediscussed in detail 48 minima observed between 1923 and 1940" but found the 
period to be constant within a second" and no indication of any oscillation with an am- 
plitude in excess of 5 minutes or a period of less than a century.'® Thus it appears that, 
except for a possibly abrupt lengthening of the orbital period in the years 1900-1905, 
there is no sufficient evidence of any period variation in U Cephei which could be traced 
back to dynamical perturbations. In consequence we may conclude, with some hesitation, 
that the effects generating, the asymmetry of light and (or) velocity-curve have their ori- 
gin in peculiar distribution and motions of matter in the outer parts of the stars, whose 
density is so low that the distribution of mass in the system remains sensibly unaffected. 
The detailed mechanism generating these effects is still a mystery. So long as it remains 
so, the conventional parameters defining the amplitudes, phase, or skewness of the anom- 
alous velocity-curves should be regarded as mere constants which do not readily admit 
of conventional dynamical interpretation. The really disturbing fact is that so far we 
have no idea as to how many ordinary (i.e., noneclipsing) spectroscopic binaries may 
have their elements vitiated in the same manner as U Cephei or MR Cygnior u Herculis 
—a fact which necessarily detracts from the significance of certain correlations obtained 
by statistical analyses of elements of spectroscopic binaries. 

There exists, possibly, one hopeful avenue of approach to the solution of our problem, 
and that is a detailed study of line profiles in the spectrum of U Cephei and other similar 


16 Chandler, A.J., 9, 49, 1889; Yendell, A./., 23, 213, 1903; Carrasco, Madrid Bol., 1, No. 16, 1933. 

16 A long series of photometric measurements by Wendell, rediscussed by M. B. Shapley (A p.J/., 44, 
51, 1916), being a notable exception. 

17 Unpublished. 


18 Harvard Bull., No. 914, 1940. 
19U Cephei is, incidentally, by no means the only variable for which most of the anomalous period 
oscillations seem to have just about ceased with the advent of photographic and photoelectric pho- 


tometry; Algol is another conspicuous example. 
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systems. The’ contours of the hydrogen lines of the B8 spectrum of U Cephei appeared 
distinctly unsymmetrical on the Victoria spectrograms taken one day after mid-primary 
minimum." This was questioned by Struve; but both he and Joy”® confirmed the unsym- 
metrical appearance of the lines at phases near 2 days. The settings, for any measurements 
of radial velocity, on extremely broad and asymmetric hydrogen lines of the principal spec- 
trum of U Cephei then represent a special problem; and the writer is inclined to believe 
that conspicuous discrepancies between spectrographic orbits derived by Carpenter,” 
Struve, and the writer” are largely due to this cause. 

Two final remarks remain to be added. First, it is very unlikely that, in U Cephei at 
least, any of its peculiarities could be ascribed to the presence of a third body. Struve is 
quite right in pointing out that the spectroscopic evidence on this point is still unsatis- 
factory. A stronger argument against the presence of any third body in U Cephei fol- 
lows, however, from the photometric evidence. If any third star were present, the orbital 
period of the eclipsing pair would have to exhibit regular oscillations due to the light- 
equation in the third orbit. From the writer’s analysis of minima in 1923-40" it appears 
reasonably certain that any such oscillation, if it exists, must have an amplitude of less 
than 5 minutes, or a period in excess of 100 years, or both. Either alternative would, 
however, imply that the mass of the hypothetical third body should be less than 0.10 
and its luminosity probably too small to exhibit any noticeable effect. 

The second remark concerns the speed.of axial rotation of the B8 component of U Ce- 
phei. This star, according to Struve, rotates much faster than it revolves, and he con- 
siders this to be the spectroscopically most firmly established fact of his study.” To this 
one might add that the existing observational evidence seems consistent in showing that 
the components of spectroscopic binaries revolving in eccentric orbits rotate as a rule 
faster than they revolve*4—usually the faster, the greater the spectroscopic eccentricity. 
U Cephei thus seems to present no exception. The fact that its velocity-curve between 
eclipses is thoroughly deformed by effects of unknown origin may, however, cast some 
doubt also on the interpretation of the rotational effect within minima. A development 
of the suggestion that abnormally large rotational effects observed in eccentric spectro- 
scopic binaries may also be partly spurious, as a result of the same causes which deform 
the velocity-curves between minima, must, however, be left to a more courageous in- 
dividual. 


20 In a letter to Dr. Shapley, which the present writer had the privilege of seeing. 
Ap.J., 72, 205, 1930. 


°2 Harvard Bull., No. 914, 1940. The orbital elements in this paper were derived under the assumption 
of a circular orbit, since no material was on hand to check it. 


23 Personal communication. 
24 A table of such systems was compiled by the writer in Harvard Circ., No. 443, 1940 (Part I, sec. 3). 
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NOTES 


COLOR INDICES OF PROPER-MOTION STARS 


In the course of the determination of the color indices of proper-motion stars in the 
Northern Hemisphere, six stars were found whose color indices appear to be so much 
less than expected that it seems highly probable that they are either white dwarfs or 


TABLE i 

Name a (1900) 5 (1900) Pg. Mag 7 p ne. 
L 1011-71... ..... 0530™5 + 1°21’ 16.0 0"41 205° 0.3 
4 41.6 32 59 .36 - 103 
6 20.4 66 29 12.0 .40 188 
2261-24. 11 51.6 18 56 15.8 .32 274 
19 07.8 6 34 12.8 220 
19 35.5 + 8 40 14.1 0.21 209 0.1 


“intermediates.” If they were stars of an absolute magnitude normal for their color, 
their linear tangential velocities would be in excess of 500 km/sec. Data for these stars 


are given in Table 1. 
W. J. LuyTEn 


P. D. Jose 


J. F. Foster 
STEWARD OBSERVATORY 
UNIVERSITY OF ARIZONA 
December 22, 1943 


ERRATA 


In the list of dwarf M stars in Volume 97, page 383, Table 1, the right ascension of star 
No. 42 should be 13528"0 and not 125280. 
A. N. VyssoTsky 


On page 71, Volume 99, line 35 from the top, read “VV Cephei” instead of “VV Cygni.” 


On page 115 of Volume 99, in line 17 of the article “Astrophysical Research in France in 
1940-1942,” read ‘‘120-cm reflector”’ instead of ‘‘120-inch reflector.” 
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